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Finding 
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ence sponsored by the University of California at Los Angeles in cooperation with 
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Propagation Studies Using Direction-Finding Techniques' 
Edgar C. Hayden 


(August 26, 1960) 


The most persistent and difficult problem in radio direction finding has been measure- 
ment of the direction of arriva] of an incident signal field under circumstances where multi- 


path propagation is possible. 


In the HF band, the ionosphere plays a predominant role in 


the propagation of radio waves, and in this region several mechanisms exist which promote 


splitting of a radio signal into numerous components. 
It would be of value to have a more thorough knowledge of the characteristics of the 


individual signal components. 
multicomponent signals. 


In this paper two techniques are described for the study of 
One involves use of pulse transmissions to effect “time-of-arrival” 


resolution; the other involves use of a highly directive antenna system to effect “‘direction- 


of-arrival” resolution. 


Results of the application of these two techniques in specific instances are presented. 


1. Introduction 


In the HF band (3 to 30 Me/s) long-distance 
transmission is possible with moderate power levels 
because of the existence of the ionosphere. How- 
ever, while the ionosphere makes such transmission 
possible, it at the same time has characteristics which 
make difficult the measurement of signal direction 
of arrival and the extrapolation of such measurements 
for position fixing or propagation studies. The 
chief of these difficulties is that frequently there are 
in fact several possible ray paths between transmitter 
and receiver, and these paths do not lie exactly in 
the plane of the great circle passing through the 
transmitter and receiver locations. The former phe- 
nomenon complicates the direction-of-arrival meas- 
urement problem, and the latter the problem of 
extrapolation of the direction-of-arrival measure- 
ments to the point of origin of the signal. 

The mechanisms giving rise to multiple paths 
through the ionosphere; i.e., the stratification of the 
region, magneto-ionic splitting, and irregularities of 
structure, are well known in principle, if not in details 
of behavior. 

As aconsequence of the multiple signal rays arriving 
at a receiving point, the field in the neighborhood of the 
receiving antenna exhibits “interference” or “fringe” 
effects. Most direction-measuring techniques give 
an indication which approximates the normal 
to the equiphase surface in the immediate neighbor- 
hood of the antenna system. In such an interference 
field, the equiphase surfaces are corrugated rather 
than plane, so that the indicated direc tion of arrival 
depends on the location of the instrument in the 
field. Since the field pattern moves about as the 
relative phases of the various rays change, the indi- 


1 Contribution from Radio Direction-Finding Research Group, Department 
of Electrical ne pg University of Illinois, Urbana, III. 
Paper presented at the Conference on Transmission Problems Related to 
High-Frequency Direction Finding, at UCLA, June 21-24, 1960, 








cation will fluctuate. Thus, for multiray fields, the 
normal to the equiphase surface is not a ‘satisfactory 
definition of “direction of propagation”’ for purposes 
of ray retracing. 

The nature of the problem is illustrated in figure 1. 
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Ficure 1, Wave interference pattern for ray separation of 30°. 
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When two waves, each of which is assumed to be a 
plane wave with its wave normal lying in the plane 
of the paper, arrive at a receiving site from different 
directions, a standing wave field is produced. If 
(as on a transmission line) the two are going in op- 
posite directions, the maxima and minima are each 
spaced one-half wavelength, and lie in planes per- 
pendicular to the directions of propagation. 

The equiphase surfaces are also planes perpendicu- 
lar to the direction of propagation. If the angle be- 
tween the directions of propagation of the waves is 
reduced from 180° toward zero (it is 30° in figure 1), 
the maxima and minima remain plane surfaces 
which are parallel to the bisector of the angle between 
the wave normals or rays. Their spacing, however, 
increases. It is 


d 


Su—-mM=Sm-m=5 sin e 


Sy—m is the distance between maxima 
Sm—m 18 the distance between minima 
\ is the wavelength 


2a is the angle between the wave normals 
orrays (1) 


Unlike the maxima and minima, the equiphase 
surfaces do not remain plane. They become corru- 
gated (or, if the two waves are of exactly equal 
strength, broken). “On the average,” they are 
perpendicular to the wave normal of the stronger 
ray. If the waves are assumed to be in phase at the 
origin, the intersections of the equiphase surfaces 
with the X Y plane (the plane of the paper) are given 
by the following expression: 


(2x X cos 2) mr 
n° Jo 
r sIn a 


nN 1+h 
ak OO arctan E = ta 


ee ee ee ee (2) 


It can be seen in figure 1 that the normal to the 
phase front varies from a direction between the 
wave normals (but closer to the normal to the 
stronger wave), to a direction outside the angle 
between the wave normals (but away from the 
normal to the weaker wave). Thus, as the relative 
phase between the waves changes, the normal to the 
equiphase surface at a fixed point will swing back 
and forth, across the direction of the normal to the 
stronger wave. As the two signals approach equality 
in strength, the swing approaches 90°. 

When the angle between the wave normals (2a) 
is fairly small, the distance between the corrugations 
of the phase surfaces is several times as great as the 
dimensions of asmall aperture system. For example, 
if @ is 5°, the distance between adjacent corrugations 
is over 5.54. This is about 22 times the maximum 








useful size * of a 4-element Adcock system, and about 
5.5 times the maximum useful size * of a regularly 
spaced 8-element Adcock system. All present small 
aperture direction-finding systems give an indicated 
bearing which is equal to, or very nearly equal to, 
the normal to the equiphase surface (this is the 
commonly used definition of the direction of propa- 
gation). Thus these small systems can give direc- 
tional indications which are considerably (up to 90° 
in the extreme case) in error. The indicated bearing 
can show considerable fluctuation even though the 
actual ray arrival angles are steady. 

In figure 2 the error in indicated bearing is shown 
graphically for the particular case in which 2a is 
30°. Curves are given for several values of the 
relative amplitude of the waves. The error curves 
for other values of « are similar in character. Note 
that as the relative phase varies, the deviation to 
one side of the stronger ray is less than that to the 
other side. 

Similar effects occur in systems using arrays of 
large aperture, though the error fluctuations can in 
general be confined to a smaller range. They may, 
however, be considerably greater than the true 
fluctuations in direction of arrival of the rays. 

Thus, under ordinary circumstances, direction- 
finding systems of both large and small aperture are 
subject to errors in bearing indication when the 
incident signal is multirayed. If, then, the char- 
acteristics of an individual ray are to be studied, 
some means must be provided for ray separation; 
i.e., for taking the signal apart into its components. 
Three possibilities are immediately apparent: 


(1) Wait until only one ray is present before ac- 
quiring data. 

(2) Take advantage of differences in ray path 
length to effect a resolution on the basis of 
time of arrival. 

(3) Take advantage of differences in ray path direc- 
tion to effect a separation on the basis of direc- 
tion of arrival. 


The first of these very severely restricts the situations 
under which data can be acquired. In practice it is 
of limited utility. The second requires the use of 
signals consisting of shurt pulses, the resolution 
capability being determined by the pulse duration. 
The third requires large antenna systems of high 
directivity. The resolution capability is limited by 
the beamwidth of the array radiation pattern. 

The latter two of the techniques listed are both 
useful. Each covers some situations in which the 
other fails. In particular, the time-of-arrival tech- 
nique is often satisfactory for resolution of the normal 
modes of propagation, but fails under anomalous 
conditions because of the great multiplicity of raylets 
and the extremely short pulse duration which would 
be required for resolution. On the other hand, the 


3 The maximum useful size of these systems is determined by the maximum 
permissible spacing error. The figures above are for systems whose maximum 
spacing error is limited to about 2°, 
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Ficure 2. Wave interference bearing error for two-ray signal with ray separation of 30°. 


direction-of-arrival technique can often give useful 
information under anomalous conditions because the 
angular spectrum of such signals is frequently of 
large extent, but may fail under normal conditions 
because the angular separation of the rays is small. 
Some of the advantages of each technique have been 
exploited. The results of two specific experiments 
are reported in the following sections. 


2. Ray Resolution by Time-of-Arrival 


The technique used in this experiment for taking 
apart a multiray signal involves exploitation of the 
differences in time of arrival of the various rays. 
Such a technique, while not suited to the general run 
of signals, is well-suited for use under controlled 
circumstances. With a cooperative transmitter, the 
signal modulation characteristic can be chosen to 
facilitate the receiving problem. 

For resolution of multiray signals a short pulse is 
required. Since direction-finding techniques must 
often be used with transmissions which are modulated 
by waveforms other than pulses, it would be useful 
to provide a signal to simulate other sorts of trans- 
missions as well. Comparisons would thus be facili- 
tated between measurements of the directional 
characteristics of the component modes, and the 
directional indications given on the composite signal 
by various instruments. Such a signal can be pro- 
vided by transmitting alternately a pulse short 





enough to resolve the component modes and a pulse 
long enough that the received pulses for the various 
modes overlap an appropriate amount. 

Choice of a length for the short pulse is a 
compromise between the resolution required on one 
hand, and signal-to-noise ratio, available receiver 
bandwidth, Federal Communications Commission 
regulations, and the overcrowded state of the HF 
band on the other hand. The latter two considera- 
tions seemed to be the deciding factors. In order to 
avoid interference from other transmitters on ad- 
jacent channels, and to avoid interference with other 
users of the spectrum, the bandwidth was limited to 
25 ke/s. For double-sideband transmission the 
minimum usable pulse width under these conditions 
isabout 80 usec. Fortunately this was short enough, 
for the transmitter location chosen, to resolve the 
grossly different modes, though it was only occasion- 
ally that the twomagneto-ionic components of a given 
mode could be resolved completely. 

The greatest expected difference in time of arrival 
between modes of any consequence was of the order 
of 4 msec. A long pulse of 12 msec duration was 
considered adequate to insure overlap of all interesting 
modes, and to provide a composite signal of sufficient 
duration. 

The modulation waveform chosen for the trans- 
mitted signal is shown in figure 3a. If a signal of 
this configuration were to travel to a receiving loca- 
tion via two paths, the structure of the received signal 
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Ficure 3. Signal envelope waveforms. 
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a. Transmitted waveform. 


b. Example of a received waveform; two modes, amplitude of first 3/2 that of 
second, relative phase 180°, difference in arrival time 1 msec. 


ce. Expanded sweep, short pulses only. 


might, for example, be as shown in figure 3b. The 
long pulse shows three regions, a leading edge in 
which only the first-to-arrive ray is present, a trailing 
edge in which only the last-to-arrive ray is present, 
and a central region in which a composite of the rays 
is present. Thus the long pulse permits some decom- 
position of the signal, but also permits study of the 
composite signal. If the difference in time of flight 
for the two rays is more than the duration of the 
short pulse, then two distinct short pulses will be 
observed (fig. 3c). In some cases each of these short 
pulses can be further split into magneto-ionic com- 
ponents, thus allowing an additional step in the 
decomposition of the signal. 

Because of its ability to work on pulse signals, the 
directional measurement system chosen was a con- 
ventional Adcock antenna system followed by a 
twin-channel receiver and a cathode-ray-tube 
goniometer. The elements of the system are shown 
in figure 4 
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Ficure 4. Direction-finding system, block diagram. 
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In addition to the bearing, information about the 
state of polarization was needed. A horizontal loop 
antenna was installed in the center of the antenna 
array. It can be seen in figure 4. Its output was 
received by a third channel of the direction-finding 
receiver. Since the direction-finding antenna was 
vertically polarized, and since signals were expected 
to come from an easterly direction, the output of the 
E-W channel of the direction-finding receiver was 
used, with the output from the loop receiver, to dis- 
play the state of polarization on a cathode ray tube. 
Quantitative measurement of the state of polariza- 
tion was not needed, so this simple system sufficed in 
spite of the different directional characteristics of the 
horizontally and vertically polarized parts of the 
system. 

Three cathoderay tubes, asshown diagrammatically 
in figure 5, were used in the display s system, one to dis- 
play bearing, a second polarization, and the third 
signal structure. The third tube was a two-gun tube, 
one gun being used to display the vertically polar ized 
signal, and the other the horizontally polarized signal. 
In addition, each gun was double-traced, alternate 
traces serving to display signal structure and timing 
scale. A mode selection gating pulse, variable in 
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Ficure 5. Key to interpretation of signal display. 
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duration and in position along the structure trace, 
was provided to select portions of the signal for dis- 
play on the bearing and polarization indicators. The 
duration and position of this gate are shown as a 
“step” in the timing scale base line. In some 
cases an expanded edie sweep was used, showing 
only the short pulses as illustrated in figure 3c. The 
expansion factor is roughly 10 X. 

Transmitter pulse repetition rate and receiver dis- 
play rate were synchronized by deriving each of them 
from a stable 100 ke/s oscillator. The basic pulse 
repetition rate was 25/sec. 

The displayed information was recorded with a 
specially constructed pulse-driven 16-mm camera, 
whose operation was controlled by the timing pulse 
train. The next few figures show examples of typical 
signal displays. The figures were made from indi- 
vidual frames of the actual 16-mm film recording. 
Reference to figures 3 and 5 will aid in understanding 
figures 6, 7, 8, and 9. In the figures showing the 
complete signal waveform, the camera shutter was 
open for one data cycle (2/25 see because of the 
double-tracing of the structure display to present 
both timing and data traces). The shutter was 
opened and closed in the sweep flyback periods. 
Thus one timing and signal structure trace, 
but two successive bearing and polarization traces, 
were recorded. In the figures showing the short 
pulses only, the shutter was open for five complete 
data cycles (10/25 sec), causing five timing and 
signal structure traces, but ten bearing and _ polar- 
ization traces, to be recorded. The multiple exposure 
of the short pulse signals was of material help in 
reducing the effects of noise on the bearing display. 

Figure 6 shows a two-mode signal (one-hop E 
and one-hop F’). The two short pulses are visible 
near the left end of the structure display traces. The 
long pulse shows three portions, a portion where only 
the first mode is present, a portion of composite 
signal, and a portion where only the second mode is 
present. The mode selection pulse is set to select 
the first mode and a portion of the composite signal. 
Bearings and polarization ellipses for these two por- 
tions of the signal appear superposed on the appro- 
priate indicators. The two superposed images on 
each of the indicators may be differentiated by com- 
parison of their relative dimensions with the relative 
amplitudes of the signal components as shown on the 
structure display. Note that only the vertically 
polarized (upper) component of the structure display 
is of use in interpreting the bearing display, since the 
direction-finding antenna was vertically polarized. 
Both components of the structure display are useful 
in interpreting the polarization display. This series 
of five pictures shows one fading cycle, as the relative 
phase of the components takes on successive values 
of approximately 0°, 90°, 180°, 270°, and 360°. The 
steadiness of the bearing and polarization indications 
on the single mode should be contrasted with the 
variability of these indications on the composite por- 
tion of the signal. The behavior of the bearing 
indication on the composite portion of the signal 








should be compared with the wave interference 
bearing error curves in figure 2. 

Figure 7 shows a similar fading cycle, but the mode 
selection pulse is set to pick out the second mode, 
‘ather than the first, and a portion of the composite 
signal. 

Figures 8 and 9 show a one-hop F-layer signal. 
The structure trace is in the “fast sweep”’ position so 
that only the short pulses are displayed. 

In figure 8a only the ordinary magneto-ionic com- 
ponent is present, while in 8b only the extraordinary 
component is present. These two can be differ- 
entiated by their characteristic polarizations, which 
are quite stable from day to day. Note the differ- 
ence in indicated bearing. In figure 8c both com- 
ponents are present, and are nearly resolved. Both 
characteristic polarizations are evident. The bear- 
ing display is not easily interpretable because the 
bearings of the two components are so nearly alike 
that the indications are not distinct. 

Figure 9a shows a case in which both magneto- 
ionic components are present, but are incompletely 
resolved. In the polarization display the figure 
starts to build up in one of the characteristic polari- 
zations, then switches to the other. During the 
transition, or overlap, the vertically polarized com- 
ponents are in phase, but the horizontally polarized 
components are out of phase. This causes the 
vertical hourglass shape of the polarization display 
and propeller shape of the bearing display. This 
also shows up in the splitting of ‘the horizontally 
polarized pulse, but not of the vertically polarized 
pulse. The bearing trace is in effect, swinging from 
the bearing of one component to that of the other 
through a small angle. In figure 9b, the situation 
is similar, except that the vertical components are 
out of phase and the horizontal components are in 
phase. In this case the hourglass figure on the 
polarization display is horizontal. The vertically 
polarized pulse is split, but the horizontally polarized 
one is not. The bearing trace swings from the bear- 
ing of one component to ) that of the other through an 
angle of nearly 180°, giving the effect of a split 
bearing. The wave interference bear ing error curves 
(fig. 2) will assist in interpreting the bearing displays 
in “figures 9a and 9b. Note in figure 2 that if two 
components are in phase, and their relative amplitude 
swings from zero to infinity, the bearing will swing 
from that for one signal to that for the other through 
a small angle, while if the components are out of 
phase, the bearing will swing through the supplement 
of this angle. It is this action that is taking place 
in figures 9a and 9b. Figure 9c shows a case where 
the overlap is so nearly complete that neither the 
horizontally nor the vertically polarized pulses, are 
split. However, the presence of the two components 
can still be detected by the X-shaped figure in the 
center of the polarization display (compare with 
figs. 8a and 8b in which only one component was 
present). 

The results reported here are for a particular 
single path, and a particular frequency. Though 
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they are for this reason somewhat specialized, they 
do show some of the phenomena that are observable, 
and do demonstrate the power of the technique. 
The transmitter was located at Columbus, Ohio, and 
the receiver at Urbana, Illinois. Both stations lie 
within a few minutes of the 40th parallel of latitude. 
They are separated by a distance of 450 km. The 


os 


FiGureE 6. 








frequency used was 5155 ke/s. The transmitter 
delivered a peak power of 600 watts to a horizontal 
dipole whose axis had a north-south orientation. 
The great circle bearing of the transmitter at the 


receiver was 88% degrees. Data were taken only 


during the hours between evening twilight and mid- 
night. 


Typical composite signal with one-hop E and one-hop F modes, one-hop E and composite selected for presentation on 
bearing and polarization displays; successive frames illustrate slip in phase (o) between rays over one fading cycle. 


a, =0°; b, ¢=90°; c, ¢=18)°. 
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Ficure 6. Typical composite signal with one-hop E and one-hop F modes, one-hop E and composite selected for presentation on 
bearing and polarization displays; successive frames illustrate slip in phase (¢) between rays over one fading cycle—Cont. 


d, ¢=270°; e, ¢=360°. 


At the times the data were acquired, the observed 
MUF for this path for one-hop F-layer transmis- 
sion was above the operating frequency during the 
entire daily transmission period. The MUF for 
normal one-hop E-layer transmission dropped below 
the operating frequency shortly after the beginning 
of the transmission period. Mode identification 
was accomplished by continuously observing the 
signal structure pattern throughout the entire late 
afternoon and evening, from the time when both 
the normal one-hop E- and F-layer modes were 
present, through the period of breakup of the normal 
E-layer mode. The polarization states characteristic 
of the two magneto-ionic components of the one-hop 
F-layer mode were determined by observing the 
polarization display generated by the short pulse 
through that period where the two components inter- 
changed order of arrival as the normal E-layer 
breakup occurred. Only on rare occasions was 
appreciable splitting of the E-layer mode observed. 

The significant observed transmission modes are 
shown in terms of relative time of arrival in figure 





10. The times of arrival are all referenced to the 
one-hop F-layer mode (designated C), since that 
was the only mode regularly present, and always 
identifiable. No means was available for reference 
to the time of transmission. In addition to the 
normal one-hop F-layer mode, transmission by a 
one-hop nighttime E-layer mode (designated A) was 
observed during most evenings, for varying periods 
of time. In fact all the data recorded on E-region 
signals was obtained on such sporadic transmissions. 
The mode designated B was tentatively identified as 
a one-hop reflection from the F; region. Its relative 
time of arrival suits that situation, and its direction 
of arrival was essentially on-path. The other modes 
were not positively identified because their occurrence 
was not sufficiently frequent or prolonged to permit 
many observations to be made. However, the mode 
designated F was probably a 2-hop F-layer reflection. 

Bearing information was read from the film record- 
ings, frame by frame, and punched onto paper tape 
for analysis by the Illiac, the University of Illinois’ 
electronic digital computer. Well over 12,000 in- 
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Figure 7. Typical composite signal with one-hop E and one-hop F modes, one-hop F and composite selected for presentation on 
bearing and polarization displays; successive frames illustrate slip in phase ($) between rays over one fading cycle. 


a, ¢=0°; b, ¢=90°; c, ¢=180°. 
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FIGURE 7. 


bearing and polarization displays; successive frames illustrate slip in phase (@) between rays over one fading cycle 








Typical composite signal with one-hop E and one-hop F modes, one-hop F and composite selected for presentation on 


Cont. 


d, ¢=270°; e, ¢=360°. 


dividual frames were read. Each sample bearing 
determination consisted of readings or “‘snaps’’ from 
a series of frames, ranging from 10 to 286, and 
averaging about 40. Figure 11 gives examples of 
bearing samples of various types. The range of 
variability, classed by propagation mode, of the 
snaps within a sample is shown in figure 12. This is 
a measure of the rapid (second to second) fluctua- 
tion range. Information is given for one-hop prop- 
agation via the nighttime E-layer, the F-layver with 
ordinary and extraordinary components both present, 
the F-layer with ordinary only, and the F-layer 
with extraordinary only. On the graph bars are 
marked the minimum, mean, and maximum values 
of the RMS deviations from the sample mean of the 
snaps within a sample. 

Of interest is the relatively high short-time direc- 
tional stability of the nighttime E-layer mode. 
This was true in spite of the fact that both magneto- 
ionic rays were present. They were invariably too 
nearly coincident to resolve. This stability seems 





indicative of a nighttime E-region having structural 
regularity over a considerable area. The structure 
must be such that it produces a very nearly specular 
reflection, rather than a structure which is patchy, 
thus producing scattered, or nonspecular, reflections. 
The rather sudden appearance and disappearance 
of the E-layer signal indicate relatively rapid changes 
in the E-layer ion density or height, causing rapid 
changes in the “skip” distance. Since vertical 
soundings show little variation in the height of 
reflection from the E-region, ion density variations, 
having structural uniformity over a fairly large 
area, seem the more probable explanation. 

The large variability in the undifferentiated 
F-layer mode, compared to the single magneto- 
ionic components of that mode, is almost surely 
due to wave-interference effects (resulting in cor- 
rugated phase fronts) between the two magneto- 
ionic rays. The paths traveled by the two magneto- 
ionic rays differ considerably more for F-layver trans- 
mission than for E-layer transmission. This illus- 
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Figure 8. One-hop F mode, short pulses only, expanded sweep. 


a, Ordinary ray only; b, Extraordinary ray only; c, Both rays, essentially resolved. 
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a, Vertical components in phase, horizontal out of phase; b, Vertical components out of phase, horizontal in phase; c, Vertical components in phase, horizontal out of 
phase; two rays nearly coincident. 





mode, short pulses only, expanded sweep, ordinary and extraordinary rays both present, partially resolved. 
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Ficure 10. Observed mode structure: note that on a given day 





not all modes observed were necessarily present simultaneously. 





trates the importance of completely decomposing the 
signal to its basic rays before drawing conclusions 
about ray directional characteristics from the 
variability of the bearing indication. 

Distributions of the sample means about the great 
circle path, taken over the entire duration of the 
experiment, are shown in figure 13. Notice how 
tightly the #-layer signals cluster about the great 
circle path. The long-time mean for this mode is 
only 0.6° from the great circle bearing, and the 
RMS deviation of the sample means about this 
long-time mean is only 1.6°. 

The F-layer bearing distribution chart shows an 
interesting split near the great circle bearing value. 
The reason for it can be seen in the distribution 
graphs for the ordinary and extraordinary com- 
ponents separately. The ordinary component shows 
a definite bias to the north. In fact, no samples in 
this mode showed means south of the great circle 
path. The extraordinary component, on the other 
hand, shows a definite bias to the south of the great 
circle path. There are a few entries to the north 
of this value, but examination of the data in detail 
reveals that these instances accompany unusually 
large deviations of the ordinary component to the 
north. Thus, at these times, both components 
showed an unusual northerly shift. The Basic 
Radio Propagation Prediction charts issued monthly 
by the National Bureau of Standards indicated that 
at these times a horizontal gradient of ion density, 
increasing to the north, might have been expected. 

Examination of the phase refractive index of the 
ionospheric region shows that the sense of the bias 
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FicurE 11. Examples of bearing samples; each sample contains a number of snaps, each snap being read from one frame of film. 
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of the two components is as should be expected. 
No attempt to check theoretically the magnitude of 
the bias has been made. Notice that the biases of 
the two components are not equal, the ordinary 
component bias averaging about 9°, and the extraor- 
dinary component bias averaging about 4°. This 
would be expected as a consequence of the tilt of 
the earth’s field at this latitude. Note also that the 
average of the means of the ordinary and extraor- 
dinary components taken separately is only 0.5° 
from the mean for the undifferentiated mode (the 
data for the two magneto-ionic components sepa- 
rately is not included in that for the composite mode). 
Thus if directional measurements on signals in this 
mode were taken over a period of time, they would 
show a northward bias. If a position fixing task 
permitted long-time measurements, such a bias might 
be taken into account, thereby improving the fix 
accuracy. 


3. Ray Resolution by Direction-of-Arrival 


When the component rays of a multimode signal 
arrive from directions which differ sufficiently, they 
can}be resolved by a directive antenna system and 
displayed separately. The Wullenweber direction- 





finding system at the University of Illinois is such a 
system. It consists of an array of 120 vertical antenna 
elements (folded monopoles) disposed on a circle 
994 feet in diameter on the outside of a reflecting 
wire screen 955 feet in diameter. The elements are 
used 40 at a time, those in use being coupled to the 
receiver through a phasing system which forms a 
fan-shaped radiation pattern with a small azimuthal 
beamwidth. A rotary capacitive switch is provided 
to scan the directive pattern in azimuth by succes- 
sively connecting a new element at one end of the 
set of 40, while dropping one at the other end. 
The capacitive switch is motor driven at a rate of 900 
rpm. Thus a complete azimuth scan is obtained in 
ts sec. The switch can also be manually set. 

The output signal from the receiver is visually 
displayed, causing a radial deflection of the spot on a 
‘athode ray tube whose deflection yoke rotates in 
sychronism with the scanning switch. 

The system is designed to operate in the frequency 
range from 4 to 16 Me/s. The half-power beam- 
width ranges from 14° at 4 Me/s to 3° at 16 Me/s. 

The signals from radio station WWYV are particu- 
larly convenient for study because of their simplicity 
and their continuous presence on several frequencies, 
and because of the favorable transmitter location 
with respect to the receiving site. The direct path 
length isabout 1000km. The variety of propagation 
phenomena whose directional characteristics are 
observable is gratifying, and not a little surprising. 
The system has been particularly successful in con- 
veying a description of the angular spectrum of 
signals which are too complex to be handled on 
smaller direction finding systems. Some examples 
of such signals are described in the following para- 
graphs. In interpreting the figures, note that the 
great circle bearing of the transmitter is 93°. 

One of the interesting phenomena observed is the 
nature of the signal bres akup as the MUF drops down 
through the operating frequency. This is presently 
observable on the 10 Me/s transmissions of WWV 
almost every evening. Figure 14 shows a sequence 
of photographs of the azimuth display taken during 
an afternoon and evening. The first in this series 
shows the nature of the display on a single-ray signal. 
It is essentially a polar plot of the array radiation 
pattern for this frequency (10 Me/s). The irregu- 
larity in the nose of the pattern is a commutation 
defect. Multiple ray arrivals at angles appreciably 
off the great circle bearing are evident in the inter- 
mediate photographs. Multir: ay character can be 
observed even where the rays are not resolved by the 
broadening of the response pattern. In the last 
photograph only scatter signal is present. It is, 
however, appreciably above the noise level. While 
the details of the signal behavior vary from evening 
to evening, the general pattern is usually similar to 
that shown here. 

In figure 15 areshown two examples of fairly strong, 
but highly non-specular incident signals. They are 
from WWV on 10 Me/s. Such signals are charac- 
terized by very high fading rates (up to 20 or 30 fades 
per second). 
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Fiaure 14. Sequence of signal 
displays showing signal break- 
up as MUF decreases, pass- 
ing through operating fre- 
quency; signal is from WWV 
on 10 Me/s, great circle 
bearing is 93°. 
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Figure 15. Examples of anomalous signals propagated via spread-F and scatter structures. 


a, Spread-F signal; b, Scatter signal. 


Signals of the type shown in figure 15a seem to be 


associated with the existence of spread-F over eastern 
Canada. Such displays are frequently observable, 
and may last from a few minutes to several hours. 

Energy may be incident overan angularsector ranging 

from 5° to 60°. 

Using some reasonable assumptions, it has been 
possible to estimate the point of reflection for some 

signals of this type. It is often possible to receive a 

normal on-path signal from WWV on 5 Mc/s at times 

when anomalous displays are observed on 10 Me/s. 

By setting the steerable beam manually, it has been 

possible to measure, as a function of direction of ray 

arrival, the difference in time of flight between the 
time ticks propagated in a regular mode and that of 
the ticks propagated in the anomalous mode. Assum- 
ing that the anomalous reflection occurs at F-region 
heights, knowledge of the direction of arrival and 
excess time of flight permits location of the reflection 
point. A number of such determinations are shown 
in figure 16. Each series of points connected by a 
ine represents a number of measurements on one 
display of the phenomenon. A period of the order of 
20 minutes was required to obtain the measurements 
in each set. The different sets were measured on 
different occurrences of the phenomenon. A number 
of parallels of geomagnetic latitude have been added 
to the map for reference to show the tendency of the 
reflection point loci to lie, at any one time, along 
parallels of geomagnetic latitude. 

Figure 15b shows another type of non-specular 
signal having even less regularity and higher fading 
rate than that of figure 15a. Such signals have 
been called “scatter” signals. There are times when 
the individual raylet behavior in these signals is 


Attempts have been made to plot reflection point 
loci for signals of this sort. However, because of 
the very rapidly changing configuration of the dis- 
play, these attempts have met with much less 
success than was achieved on signals such as shown 
in figure 15a. 


All of these phenomena are being made subjects 
for study in greater detail. 


4. Conclusion 


The two experimental studies described demon- 
strate the power of direction-finding techniques in 
studying propagation phenomena. A particularly 
promising extension of this work is the tracking of 
a frequency-swept ionosonde with high-resolution 


DF equipment. An experiment of this type is being 
planned. 


The experiments described above were performed 
by members of the staff of the Radio Direction- 
Finding Research Laboratory, Department of 
Electrical Engineering, University of Illinois, as 
parts of research programs sponsored by the Office 
of Naval Research and the Bureau of Ships. 
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Diversity Effects in Long Distance High Frequency 
Radio Pulse Propagation” * 


S. A. Bowhill * 


(August 15, 1960; revised January 4, 1961) 


Spaced antenna measurements are described, made on a 8,600-km path between Co- 


lombo, Ceylon, and Great Baddow, England, using pulsed radio signals. 


Simple interpre- 


tations in terms of # and F region multiple reflections give good agreement with observed 


delays for the various echoes. 


The correlation between echo amplitudes at spaced antennas was found to be much 
greater for a pulsed signal than for a CW signal, indicating that most of the diversity in 
long-distance CW transmission arises through phase incoherence between the various orders 


of reflection. 


1. Introduction 


A study of the ionospheric irregularities responsible 
for the fading of long distance radio signals can best 
be accomplished by oblique incidence pulse measure- 
ments. Previous workers [Wilkins and Kift, 1957: 
Silberstein, 1958] have shown that very complex 
modes of propagation may be present on long-dis- 
tance paths. Detailed theoretical study of these 
modes has been lacking until recently, when Kift 
(1960) gave a very detailed method for studying the 
entire propagation path. 

This paper describes the basic phenomena in- 
volved, and gives some experimental results obtained 
on 16.16 Me/s frequency, using transmissions from 
Colombo, Ceylon, received at Great Baddow, 
England. 

The complicated pattern of the pulse signal 
received led to a consideration of possible reflection 
mechanisms; third and other order F region reflec- 
tions were found, with F region multiples present 
on occasion. The complexity of a typical pattern 
implies that it is a difficult matter to distinguish, 
in the CW case, between ‘xding due to phase inter- 
ference between different orders of reflection, and 
fading of the reflections themselves. 

The path differences for the various modes are 
calculated in section 2, and the experiments are 
described in section 3. Section 4 compares the cal- 
culated path differences with those experimentally 
observed; the fading analysis is given in section 5. 

! The work described in this paper was performed by the author while at 
Marconi’s Wireless Telegraph Co., Ltd., Baddow Research Laboratories, Great 
Baddow, Essex, England. 

2 Paper presented at the Conference on Transmission Problems Related to 
High-Frequency Direction Finding, at UCLA, June 21-24, 1960. 


- 3 —- Research Laboratory, Pennsylvania State University, University 
ark, Pa. 





2. Ray Theory of Long-Distance Propagation 


A radio wave traveling a long distance between 
points on the earth’s surface, such as A and M 
(fig. 1), may be reflected sometimes by the F layer 
and sometimes by the F layer. In section 2.1, the 
elementary theory (ignoring ionospheric tilts) for a 
thin layer is given; in section 2.2, this is extended 
to include a sporadically occurring layer at a lower 
height than the main reflection. Section 2.3 de- 
scribes how a layer of finite thickness may be quickly 
allowed for. 


2.1. Single Thin Layer 


For any particular path, the distance d of the 
receiver from the transmitter is fixed, and the only 
variables needed to determine all the properties of 
a particular mode of propagation are the height of the 
reflecting layer, h, and the order n of the reflection, 
i.e., the number of “hops” the radio wave makes 
from the earth to the layer and back over its path. 
If the angle of arrival, the angle the wave makes 
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Fiaure 1. Typical reflection geometry. 
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with the horizontal at the earth’s surface, is 6, one 
may write the relations 


d= Rsin(d/2nR)tan(6+d/2nR)—R-+ Reos(d/2nR), 
(1) 


p=2nRsin(d/2nR)sec(6+d/2nR)—d, (2) 


where p=path difference between the reflected wave 
and a wave traveling directly over the earth’s 
surface. 

The eqs (1) and (2) may be regarded as a para- 
metric equation for p in terms of h. 

Numerical application of the formulas gives the 
set of curves shown in figure 2. Three sets of 
curves are shown on this graph. The first set, for 
integral values of n, the order of reflection, shows 
the variation of the time delay p with the height of 
the reflecting layer, and it is evident that the greater 
the height of reflection, the greater will be the spacing 
between successive orders of reflection in the pulse 
pattern. 

Not all the orders of reflection shown are possible 
in any practical case. Firstly, if we neglect diffrac- 
tion effects, it is not possible for radio waves to 
arrive at the receiver with effective angles of arrival 
less than zero. A set of curves has therefore been 
drawn showing the angles of arrival of the different 
orders of reflection at various layer heights. These 
curves of constant angle of arrival happen to be 
nearly straight lines, and also are approximately 
parallel. This implies a relation between the rela- 
tive time delay of an order of reflection and its angle 
of arrival, practically independent of the height of 
reflection. The upper limit is determined by the 
electron concentration of the ionized layer. The 








greater this concentration, the steeper the angle at 
which a radio wave may be incident upon it and still 
be reflected. This is expressed by the relationship, 


t=h sec i, 


where f=radio wave frequency, 
i=angle of incidence of wave at layer= 
(n/2)—o—d/2nR, 
fvu=vertical incidence critical frequency of the 
layer. 


As the graph is constructed for a particular fre- 
quency, namely 16.16 Me/s, curves of constant f can 
be drawn in on figure 2. If the vertical incidence 
critical frequency of the layer is known, the highest 
permissible order or reflection may be directly read 
off. 

It is now necessary to compare these theoretical 
suggestions with actual ionospheric data observed 
at vertical incidence pulse sounders for the period 
in question. The path from Colombo to Great 
Baddow passes, in the main, over a region where 
there are no pulse sounders whose results are regu- 
larly available. However, use may be made of the 
control-point method, which is normally applied to 
predicting the behavior of long-distance HF propa- 
gation paths. 

The western control point is fairly near the 
ionosphere station at Graz, Austria, while the east- 
ern point is at the same latitude as Bombay, and 
not too far removed in longitude. Mean ionospheric 
data for the month of July, 1953 are available from 
these two ionospheric sounding stations, and are 
given in table 1, with allowance for the difference be- 
tween G.m.t. and the local solar time at the control 
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Figure 2. Colombo-Baddow path delays. 
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TABLE 1 
| E. control point W. control point | Mean 
Time | virtual 
G.m.t. | | | height 
Pr. | WF re | ae | 
She eee 
| 
Me/s | km Me/s | km km 
08 90° 4 as 4 | 300 | 366 
09 93 | 432 ee fee 374 
10 9.6 419 5. 2 300 | 360 
We | 25 | oe 52 | 340 | 374 
12 9.3 | 383 5.1 310 346 
13 8.8 361 5.0 | 330 346 
“ | 8&1 327 50 06|l6C Cid (at 
15 7.0 | 324 4.9 300 312 
16 5.9 | 316 5.0 310 | 38 
17 ry ee ee 5.0 a 296 
point. It is clear that the ionization density at the 


western control point is the limiting factor, except at 
1700 G.m.t. 

Corresponding values of the effective reflection 
height and limiting frequency for the western control 
point have been plotted in figure 3 for each hour in 
the afternoon. If a line is drawn vertically down- 
wards from the point corresponding to any hour, it 
should intersect the curves of constant n at the 
experimentally observed values of delay. This 
will be compared with experimental data in section 4. 


2.2. Multiple Thin Layers 


In this section a simple graphical method will be 
described for deriving the theoretically expected pulse 
pattern if multiple reflections are present. It will 
be assumed that the signal is propagated by reflec- 
tion at two sharply bounded layers, situated at 
known heights. The upper layer EIJK on figure 1, 
will be assumed uniform in properties, while the lower 





the horizontal direction. At B and J it is present, and 
reflects the ray; at D, F, H, and L it is absent, and 
the ray passes through the layer without deviation. 
The reflection at B is an “E”’ reflection; that at E 
is an “F”’ reflection; while that at J, J, and K will 
be termed an “‘M”’ reflection. 

Two points must be noted about this type of 
propagation. The first is that, provided the numbers 
of each type of reflection remain the same, they may 
be permuted in any order without affecting the char- 
acteristics of the received signal. The second is that 
the angle of departure of signal at A is equal to the 
angles of arrival at C, G, and M, and also determines 
directly the angles of incidence of the ray at the two 
layers, since they are at predetermined heights. 

The characteristics of the \/ mode of propagation 
are closely related to the EF and F reflections. In 
fact, as can be seen from figure 1, if the extension 
lines JP and JQ are drawn, its transmission distance 
GPQM is equal to that of two F reflections, less 
that of one E reflection. The total transmission 
distance GHIJKLM is also equal to that of a 2F 
reflection less that of a 1F reflection. So M reflec- 
tions will not be treated separately, but simply 
included as 2F+ IE reflections, where I denotes a 
negative order of reflection. For instance, the 
path illustrated in figure 1 is just equivalent to 1#+ 
1F+2F—-1E=3F. 

If propagated solely by reflection at the E layer or 
the F layer, let a ray make Nz or Ny hops, respec- 
tively, in traveling a distance d (the fact that Nez 
and Ny may be fractional presents no analytical 
difficulty; they may be regarded as obtained by in- 
terpolation between integral values of Ng and Nr 
for the particular value of 6). The ray therefore 
travels a distance d/Nz in one E hop, and d/N;, in one 





layer DFHJL will be taken to be inhomogeneous in | F hop. If the ray makes mE reflections and nF 
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Fiaure 3. Thick layer delays for July 3, 1953. 
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reflections, clearly 


m,n ; 
gel. (3) 
+YE Nr 

It is convenient to represent the various orders of 
reflection on a graph like figure 4. Each intersection 


| 


of integral ordinates and abscissas represents a pos- | 


sible mode of propagation. Points along the hori- 
zontal axis represent rays propagating by 1, 2, 3, 
etc., reflections at the # layer alone; similarly, on 


the vertical axis they represent reflections occurring | 


at the F layer alone. Points not lying on either axis 
represent mixed modes of propagation, partly / and 
partly F. 

Now eg (3) is just the equation of a straight line on 


this diagram, since Ng and N;, are constant for con- | 


stant 6. The full lines on figure 4 are for @=0°, 5°, 


10°, 15°, and 20°, for a height of reflection of 120 km 
for the # layer, and 250 km for the F layer. 

The points on the negative side of the FE axis rep- 
resent the M type reflections mentioned previously ; 
for instance, the point P on figure 4 is 6F'+-3E, or 
3 M reflections. A line has been drawn from the 
origin with a slope of —1; no intersections on or 
below this have any physical meaning. 

One may easily find the total time delay of a 
mixed order of reflection. If the ray 
delayed by pg or pr when propagated entirely by 
reflection at the / or F regions respectively, then its 
delay for unit distance is pg/d and p,/d; the delay 
experienced by an mE nF echo is then 


a q 
P=Pr°- rt Pr —* (4) 


The loci of constant p on figure 4 are therefore again 
straight lines. They are drawn as dashed lines, for 
delays of 0.5, 1, 1.5, 2, 2.5, and 3 msec. 


would be | 


| 
| 
| 


The highest order of F region reflection which can 
be propagated is limited by the F region critical 
frequency; the relevant frequency for each angle for 
an F region height of 250 km is inserted opposite the 
F axis on figure 4. 

Since the nature of Fs ionization is to occur in 
relatively intense patches, it will be assumed that the 
E mode is determined by its spatial extent, rather 
than its critical frequency as such. From eq (3), 
E propagation occurs for m/Nz of the path; if it is 
assumed that there must be on an average this 
amount of Es to support the mE nF mode of propa- 
gation, lines of constant fractional coverage of Es 
can be inserted on figure 4; they are shown for frac- 
tions of 0.2, 0.4, 0.6, and 0.8. 

As an example, one may take an F critical fre- 
quency of 6 Me/s and a fractional coverage of 0.2. 
All the points in the quadrilateral ABCD will be 
possible orders of reflection; they are 


3F, 4F, 5F 

LE 3F, 1E 4F, 1E 5F 

1E 4F, 1E 5F, 1E 6F 
2E 6F. 


The delays of these various orders can be interpo- 
lated between the straight lines of constant delay. 

It is evident that when Es propagation becomes 
important each order of F reflection will have a series 
of F multiples associated with it, one set arriving 
sarlier, the other set later, at the receiver. The in- 
terval between these successive orders of F reflec- 
tion varies from 0.15 msec for the 3F reflection to 0.3 
msec for the 7F reflection. If the Hs ionization is 
present over a large proportion of the path, the 


multiples associated with successive orders of F 


reflection overlap, giving a complex pattern. 
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Fiaure 4. Statistical effect of Es coverage. 
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2.3. Extension to Thick Layer 


Assuming the Breit and Tuve theorem to hold, the 
finite thickness of the reflecting layer may be allowed 
for by plotting the h’f curve on figure 3, using the 
limiting frequency scale derived previously. The 
curve ABC on figure 3 represents the h’f curve of a 
parabolic region, with the height of maximum 
ionization, h,,—=300 km, semithickness 120 km, and 
critical frequency 7 Mc/s. The delays and angles of 
arrivals of the various orders of reflection which 
would be observed, if this layer were present through- 
out the path, are given by the intersections of this 
curve with the n=3, 4 etc., lines. In the limit of a 
very thin region, the curve ABC becomes a vertical 
line, terminating sharply at the critical frequency. 

It should be noted that the curve ABC, if extended 


€ 


beyond the right edge of figure 3, will intersect the 





n=4 line again, at a higher angle of arrival. This, 
of course, represents the so-called “high ray’? some- 


times observed on oblique incidence paths. 
3. Experimental Measurements 


3.1. General Description 


The experiments to be described used transmissions 
from an Air Ministry sender at Colombo, Ceylon, on 
a frequency of 16.16 Me/s. The duration of each 
transmitted pulse was 100 usec and a rhombic trans- 
mitting aerial was used, the polar diagram of which 
(see fig. 5) had a maximum at abovt 10° to the 
horizontal. The experiments commenced on the 
pulse transmissions at Baddow in April 1953, and 
continued until the end of January 1954. Apart 
from unavoidable breaks in’ transmission, the 
signals were sent out during daylight hours on one 
day every week. At the end of January 1954, some 
special transmissions were arranged by the Radio 
Research Board on four consecutive days. The 
cooperation of the Radio Research Board and the Air 
Ministry throughout the experiments proved most 
valuable. 
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Fiaure 5. Transmitting antenna polar diagram. 


3.2. Receiving Arrangements 


The signal pulses were transmitted at a rate of 
25/sec; by the time they reached the receiver, the 
presence of different propagation modes had length- 
ened them to abovt 1 msec. Consequently, some 
form of triggered time-base was necessary to display 
the incoming waveforms on a cathode-ray tube 
screen. Various methods using the signal itself to 
trigger the time-base proved unsuccessful since the 
signal frequently faded below local and atmospheric 
noise. As the source of 25 c/s at the transmitter was 
a highly stable crystal oscillator, a similar oscillator 
was later used at the receiving station. 

The block diagram of the receiving equipment is 
shown in figure 6. The signals were received on a 
half-wave dipole antenna by a CR150 communication 
receiver, working on its widest bandwidth of 13 ke/s. 
Under some circumstances it was necessary to use a 
bandwidth of only 8 ke/s, due to the presence of 
powerful CW stations on nearby frequencies; how- 
ever the difference in the observed pulse waveforms 
when these two bandwidths were used alternately, 
suggested that little information was lost by using 
the 8 ke/s bandwidth. The receiver output was 
displayed on an oscilloscope; the 25 ¢/s triggering 
impulses for the oscilloscope time base were derived 
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Fiaure 6. Equipment block schematic. 
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from a 100 ke/s crystal oscillator, the frequency of 
which could be altered by about 1 part in 4,000 by 
means of a frequency pulling circuit. 

Two typical waveforms (instantaneous still photo- 
graphs of the cathode ray tube display) are shown in 
figures 7(a) and 7(b). These photographs are two 
of the series described in section 5.2. It is seen that 
the length of the pulse train is over 1 msec, and that 
it may be quite complex. 


3.3. Types of Echo Observed 


To study the structure of the echo in detail, a 
method of recording was used in which the incoming 
signal voltage modulated the brightness of the oscil- 
loscope trace. A photographic film, moved at right 
angles to the trace, then gave a record similar to the 
group-height-time records familiar in vertical inci- 
dence work. 
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Ficure 7. Instantaneous pulse waveforms. 








Using this method of recording, the only attention 
the apparatus required was occasional adjustment of 
the crystal frequency, to ensure that the pulse pattern 
was kept about in the center of the oscilloscope time 
base. 

Generally speaking, the signal was usually de- 
tectable from 1000 to about 1530 G.m.t.; on some 
occasions, the transmissions did not extend over the 
full period. The signal was usually at its strongest 
at about 1300 G.m.t., but no detailed measurements 
of signal strength were made. The type of pattern 
observed varied considerably with propagation 
conditions. 

Figure 8 shows a simple type of pattern, in which 
the time base was applied vertically, with delay 
increasing upward. The sinusoidal variation in 
delay, with a period of about 10 sec, is due to 
imperfect synchronization of the pulse. The vertical 
white dashes are due to pulses of static, of very short 
duration ; the striations in the trace near the left-hand 
edge are the result of CW interference from a trans- 
mitter on an adjacent channel. 

In the summer, much stronger echoes are received 
on days when Es ionization is present. As shown in 
section 2.2, its effect is to generate a group of closely 
spaced echoes, centered on each order of F' reflection. 
Depending on the geometry of reflection, the spacing 
may be from 0.1 to 0.2 msec. Figure 9(a) is an 
example of predominantly Es propagation, with only 
one order of F reflection present. The four compo- 
nents of the pulse have a mean separation of about 
0.15 msec. Figure 9(b) shows an intermediate case, 
where four F echoes are present, with a small amount 
of Hs causing each to aplit into three distinct echoes. 
In figure 9(c) the Es echoes have completely obscured 
the # echo pattern. 


4. Comparison With Simple Ray Model 


Designating the three pulses in figure 7(a) as P, 
Q, and R, in order of arrival, the time separation 
between pulses P and Q is about 0.24 msec, and be- 
tween Q and RP is about 0.28 msec. To explain these 
separations by a constant reflection height would lead 
to a height of only 200 km, from figure 2: this is 
clearly incompatible with the ionospheric data of 
section 2.1. 

One possibility for explaining the various modes is 
in terms of low- and high-angle rays. The curve 
GHTI on figure 3 represents a parabolic layer with 
height of maximum ionization 300 km, semithickness 
120 km, and critical frequency 6 Me/s. The points 
G, H, and I correspond to low and high angle third 
order F reflections, and a high angle second order 
reflection, respectively. This configuration produces 
a set of three modes with very nearly the observed 
overall separation: however, the two modes H and J 
occur with nearly identical delay, rather than the 
observed .28 msec separation. A departure of the 
layer from a parabolic form, giving a flattening of the 
“nose” of the layer, might displace point J to J, 
thereby removing the degeneracy. The ray would 
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Fiaure 8. Typical simple pattern. 


then arrive at a very low angle to the horizontal. 
On the whole, a flattening involving an equivalent 
vertical incidence group retardation of over 70 km for 
a frequency increment of less than 0.2 Mc/s seems 
unlikely, and Kift (1960) has suggested that high- 
angle modes propagated over great distances are 
unstable. However, this interpretation cannot be 
ruled out, particularly in view of the approximations 
involved in assuming uniform propagation charac- 
teristics over the entire path. 


An alternative interpretation is that two of the | 


three successive pulses of figure 7(a) have the same 
order of reflection but are different magneto-ionic 
components. On a thin-layer hypothesis (fig. 2), if 
pulses P and FR are assumed to be the same compo- 
nent for the third and fourth order reflections, the 
delay between the two, 0.52 msec, corresponds to a 
reflection height of about 300 km, and their angles of 
arrival are about 5° and 11° respectively. If pulse Q 
is assumed to be a third order reflection, its effective 
reflection height must be about 335 km in order to 
account for the additional delay, and its angle of 
arrival must be about 7°. 
counted for by pulses P and Q being the ordinary and 
extraordinary third order reflections respectively. 
Assuming the same parabolic layer as previously, the 
ordinary and extraordinary ray reflection conditions 
have been plotted on figure 3 as ABC and DEF re- 
spectively. It is suggested that the pulses P, Q, and 
R may correspond to the points A, D, and B respec- 
tively. 

Since there appears to be some doubt as to which 
of these interpretations is correct, the pulses will be 
referred to by the letters P, Q, and FR in subsequent 
sections. 


This might well be ac- | 
t=) 


5. Fading Analysis 


Using the analysis of the previous sections as @ 
basis, the principal characteristics of the fading 
observed over this path will now be described. 


5.1. Simple Echo Fading 


On a limited number of occasions, only a single 
order of reflection was present in the pattern. At 
these times, an easily interpretable fading curve 
could be obtained by photographic recording; an 
example is shown in figure 10. The upper and lower 
traces are the fading curves received at two antennas 
of a space diversity pair, 600 ft apart in a line per- 
pendicular to the direction of propagation. There is 
clearly a high correlation between these curves, indi- 
ating a distance scale of much greater than 600 ft 
in the random diffraction pattern. 


5.2. Complex Echo Fading 


When only a single order of reflection is present, 
the fading of a pulsed signal is identical with that of 
a CW signal, and only fluctuations in amplitude are 
important. If a CW signal is composed of several 
orders of reflection, however, the possibility of phase 
interference between the various echoes may give 
both a faster rate of fading and a smaller structure 
size on the ground (greater diversity in a space- 
diversity system). 

This comparison may be readily effected by switch- 
ing the receiving system bandwidth. If a wide band- 
width is used in the receiving system, so that the 
orders of reflection are resolved, each order will de- 
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scribe its own fading curve. However, if the band- 
width is narrowed to such an extent that the pulses 
completely run into one another, thus approximating 
to the CW case, then the resultant fading will be due 
not only to the fading of the orders of reflection, but 
to the effect of phase interference between them. 
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Fiacure 9. Pattern showing Es echoes. 








To illustrate this effect, figure 11 shows four suc- 
cessive stretches of fading record, in which the re- 
ceiver bandwidth was switched from 8 ke/s to 1 ke/s 
instantaneously. There is complete diversity be- 
tween the two 1 ke/s signals in 11B and 11D, but the 
8 ke/s signals of 11A and 11C display marked corre- 
lation. In order to study the rapid fading of the 
various elements of the pulse waveform in detail, and 
to relate this to the behaviour of a diversity receiving 
system, photographs were taken of the pulse wave- 
form at 2-sec intervals. Two examples of these 
photographs were shown in figure 7. In each case 
the upper waveform represents the pulse received by 
a four-fold diversity system, i.e., the mean of four- 
pulse patterns received on four separate antennas 
and receivers. The other pulse shown was that 
received simultaneously on a single antenna. 

To remove the possibility of slight random jitter of 
the pattern, the average of a number of waveforms 
was first calculated, and is shown in figure 12. With 
this as a guide, the individual waveforms were syn- 
chronized, and peak values read off for the P, Q, and 
R reflections. The mean values of these quantities 
are shown with circles on figure 12. They are greater 
in amplitude than the previously measured peaks of 
the three pulses, due to slight lack of synchronism 
in calculating figure 12. 

Typical fading curves for the various pulses for a 
single antenna and for a diversity system are shown 
in figures 13, 14, and 15. 


5.3. Diversity Comparison 


Several interesting points are immediately evident 
from figures 11, 12, and 13. Firstly the depth of 
fading for the three pulses is quite different, being 
much less for the Q reflection than for either the P 
or the R. The best measure of this depth of fading 
is the ratio of the standard deviation of the fading 
curve to its mean value. If the depth of fading is 
defined as D, then 


where x=amplitude of fading curve at any instant. 
Values of D are given in table 2 for the three pulses. 
For a single aerial, the values of D are 0.35 and 0.34 
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Figure 10. 


Spaced antenna fading curves—single echo. 
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TABLE 2 

Pulse Aerial Mean D pi U 
signal 

sec 
P ro 12. 7 0. 353 0. 439 3.0 
Diversity - : 19.1 . 367 . 466 2.9 
Q Single E 29.1 .174 . 450 2% 
Diversity - _- 43.5 173 . 465 2.9 
R Single — 15.5 . 343 . 343 rae 
Diversity - 24.5 . 306 . 065 1.5 


for the P and FP pulses, but only 0.17 for the Q pulse. 
The greatest value of D which would be expected 
experimentally would be the value for a signal with 
a Rayleigh probability distribution, which would give 
a value of 0.52. 

Though the diversity system has virtually no effect 
on the depth of fading for two of the echoes received, 
this does not mean that it would also be ineffective 
for a CW signal; since in CW reception the different 
orders of reflection may interfere destructively. The 
only condition for a diversity system to give an im- 
provement in the communication system performance 
is that the phase relation between the echoes shall 
be different at the four diversity antennas. 
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Spaced antenna fading curves—multiple echoes. 


The difference in fading speed of the three pulses 
The fading periods of the P and Q 
pulses are about 1 cycle in 5 to 10 sec. The fading 
of the RP pulses, however, is very much more rapid. 
The best quantitative measure of the fading speed is 
the correlation function of the amplitude. The first 
ordinate of this is the correlation p, between the 
amplitudes of the same pulse in successive snapshots 
(taken at 2-see intervals). These correlation coeffi- 
cients are given in table 2. If it is assumed that the 
time correlogram of the fading curves has a Gaussian 
form, 1.¢., 
p(r)=e-8.58r?/U? 
p(r)=correlation between amplitudes sepa- 
rated by a time interval 7, 
U'=average time between peaks of fading 
curve, 
= “fading period”’ of curve, 


where 


then the fading period of the pulse amplitudes can 
be calculated. These values of U are also given in 
table 2. 

Another feature of the fading curves is the effect 
of the diversity system on the fading of the various 
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pulses. If the amplitudes of the pulses at the four 


antennas were completely unrelated, then the value 
of D for the diversity combination of these four 
should only be half that for the pulses received on a 
single antenna. In fact, however, as can be seen 
from table 2, the value of D for the P and @Q reflec- 
tions is nearly the same whether a single aerial or a 
diversity system is used, indicating that for these 
pulses the amplitudes received at the diversity an- 
tenna were practically identical. Only in the case of 
the rapidly fading pulse, P, is there an appreciable 
difference in the values of D. 
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Fiaure 15. R pulse amplitude. 


A simple relation is derived in the appendix be- 
tween the depth of fading from a diversity system 
and the original depth of fading of the single aerial 
signal, in terms of the average correlation between 
the signal amplitudes received by the various diver- 
sity aerials. The simple theory also correctly pre- 
dicts the correlation between the single antenna 
amplitude and the diversity amplitude. 


6. Conclusion 


The work which has been described was originally 
intended to clarify the meaning of measurements of 
diversity at high frequencies. The fading of the 
different modes of propagation were studied separ- 
ately, and it was shown that the diversity arose both 
through diversity of the modes of propagation them- 
selves, and also through phase interference between 
them; which of these is the major source of diversity 
depends on the degree of complexity of the pulse 
pattern. 

The need to identify the pulses received led to a 
consideration of the way in which a steady signal 
was propagated over long distances; it has been 
shown that the “sky billiards” approach of geomet- 
rical optics gives an adequate explanation of the 
experimental observations. 
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It is suggested that when the propagation is simple, 
i.e., when there is only one predominant order of 
reflection present, the deversity of a CW signal is 
due to the irregularities in the reflecting layer; but 
when the propagation is complex, with many orders 
of reflection present, the diversity of the CW signal 
is due principally to phase interference between the 
orders of reflection. 


7. Appendix. Theory of n-Fold Diversity 
Combination 


A diversity signal is formed by averaging n signals 
RR, ... Rk, to give a signal Rp. It is required to 
find. 

1. The ratio A by which the standard deviation 
o of the original signals is reduced by the process of 
averaging. 

2. The correlation of the average with any one of 
the signals. It will be assumed that the correlation 
of each signal with every other signal is equal to p. 





‘ n?.R4=(RitRot+ ... +R)? 


=n?R?+no? [1+ (n—1)p] 


sles Oe 


nA? me 5 
n—1 


or p= 











0.306 
0.343 





For the FP pulse of section 5.3, A= =(0.892 


or p=0.73. 








+o? [p+(1—p)/n] 


*. Correlation between R and Rp 


__ lp t+ (1—p)/n] 


Vo?.o*[p+( (I—p)/n] 


So, taking the previous example, 
correlation between R and Rp=0.89 


This compares with the value of 0.887 calculated 
directly from the observations of section 5.3, indi- 
cating the accuracy of this theory. 
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Influence of Ionospheric Conditions on the Accuracy of 
High Frequency Direction Finding *’ 
P. J. D. Gething 


(August 22, 1960) 


The accuracy of fixes obtained by HF direction finding stations has been examined by 
means of a dispersion factor computed for each fix; this factor is a measure of the consistency 


of bearings taken from different stations on the same transmission. 


It is shown that the 


accuracy is significantly lower during times of ionospheric storms than at times when no 
storm occurred, and that the effect of the storm is mainly on fixes involving F2 layer 


propagation. 


1. Introduction 


It is natural to suppose that the accuracy of DF 
bearings taken on transmissions in the HF band 
will depend to some extent on prevailing ionospheric 
conditions. In particular, ionospheric storms might 
be expected to lead to some deterioration of accuracy. 
The present paper describes a statistical analysis 
undertaken to test this hypothesis. 

The bearings analyzed were obtained from stations 
equipped with standard Adcock direction finders. 
They were taken on both fixed and mobile trans- 
mitters, the majority of the transmitters being in 
the distance range 1,000 to 3,000 km from the DF 
stations. The true positions of the transmitters 
were not known in the majority of cases. 

The quality of each bearing was estimated by 
means of the Brooke system. In this system, a 
variance is assigned to each bearing depending on 
several factors known to be related to accuracy. 
The variance components were originally estimated 
from an analysis of the errors in some thousands of 
bearings obtained on targets of known position; 
they are regularly revised on the basis of analysis of 
current data on such targets. The system is de- 
scribed in detail in an accompanying paper by 
E. M. L. Beale. It should be noted that the Brooke 
variances do not take account of variations in iono- 
spheric conditions explicitly, although the com- 
ponents depending on the swing and range of snap 
bearings are probably correlated with such varia- 
tions to some extent. 

The bearings obtained on a given task from dif- 
ferent DF stations do not normally intersect in a 
point. An electronic computer has been programed 
to solve the statistical problem of finding the most 
probable position of the transmitter (Best Point 
Estimate or BPE) and a surrounding rectangle 
representing a 90 percent probability region. The 
Brooke variances are used in the calculation, and 
affect the position of the BPE and the size and posi- 
tion of the rectangle. The computer also calculates 


1 Contribution from Government Communications Headquarters, Chelten- 


ham, England. 
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a so-called ‘dispersion factor,’ the definition of which 
is given in section 3. This factor is a measure of the 
consistency of three or more bearings taken from 
different stations on the same task. The average 
value of the dispersion factor should be (n—2), 
where n is the number of bearings. Under adverse 
conditions the average value would be greater than 
(n—2). 

The average value of dispersion factor has been 
used throughout the analysis as a measure of ac- 
curacy. The advantage of using this factor in a 
quantitative investigation is that it can be calcu- 
lated on all tasks without any knowledge of the true 
positions of the transmitters. 


2. Definition of an Ionospheric Storm 


Various phenomena are associated with ionospheric 
storms, the most characteristic being the change in 
the critical frequency of the F2 layer (fo F2). The 
critical frequency often shows a brief initial rise, 
termed the positive phase, followed by a steep fall 
to a much longer negative phase, during which fo F2 
is below the value expected for calm conditions. 

It is conventional to classify storms as moderate 
or severe. A moderate storm is defined in the 
B.B.C. monthly propagation reports as one in which 


fo F2 lies between 60 percent and 90 percent of the 


fortnightly running hourly mean fo F2, and a 
severe storm as one in which fo F2 is less than 60 
percent of this mean. The B.B.C. reports list days 
on which severe and moderate storms occurred. 

The B.B.C. classification was used in the initial 
stages of the investigation. Later, in order to deter- 
mine the magnitude and duration of each storm more 
precisely, published data from the U.K. vertical 
sounding stations at Slough and Inverness were 
examined. The start of the storm was then taken 
to be the time at which fo F2, measured at Slough, 
fell below 90 percent of the monthly median for the 
corresponding time. The storm became severe when 


fo F2 fell below 60 percent of the same median value. 


The progress of a typical storm is shown in figure 1, 
where the Slough value of fo F2 is plotted against 
time. 
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Fieure 1. 


Indicates discontinuity in readings. 


It is, of course, rather arbitrary to define the start 
and finish of a fairly widespread storm by means of 
soundings taken at a single station. Comparison of 
Slough and Inverness data showed that the broad 
features of the behaviour of fo F2 were the same at 
these stations, but that the smaller features of the 
plotted records were not identical. Measurements 
from other European sounding stations had not been 
published when the analysis was undertaken. Since 
the local times at points of reflection for most of the 
tasks analyzed did not differ by more than about 
two hours from the time at Slough, the behaviour of 
fo F2 at Slough was probably reasonably representa- 
tive of conditions over the area concerned. 


3. Dispersion Factor 


As already explained, the computer calculates a 
BPE of the position of the target transmitter for 
each set of bearings fed to it. The observed bearin 
are then compared with the bearings which wound 
pass through the BPE, and an angular residual 
formed for each contributing DF station. Each 
residual is squared and divided by the Brooke vari- 
ance of the bearing. The dispersion factor is then the 
sum of the resulting numbers for all the bearings 
contributing to the fix. The factor is calculated by 
the computer to one decimal place. 

Simple statistical theory suggests that the disper- 
sion factors for n stations fms. have a x’ distribu- 
tion with (n—2) degrees of freedom. Details of this 
distribution may be found in most textbooks on 
statistics, e.g., Yule and Kendall, Introduction to 
the Theory of Statistics, chapter 20. The mean value 
of the dispersion factors should be (n—2) and the 
variance about the mean should be 2(n—2). 

These results depend on the assumptions usually 
made in DF analysis, that bearing errors are nor- 
mally distributed about zero mean and that the 
errors in bearings taken from different stations on the 
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Ionospheric storm on 7th, 8th, 9th, and 10th, of July 1958, readings at slough. 


same task are not correlated. It is known, however, 
that none of these conditions is strictly satisfied in 
practice. In particular, it seems that the leptokurtic 
distributions of bearing errors which occur in practice 
will lead to a variance rather greater than 2(n—2). 
It is difficult to estimate the effect on the dispersion 
factors of small systematic errors in bearings and 
of correlation between bearing errors. These points 
have not been considered in detail here, since we are 
concerned more with a comparison of factors under 
different conditions than with their absolute values. 

The dispersion factor will be denoted by 2, follow- 
ing the notation used by Beale and others. No 
convenient name or symbol exists for the quantity 
Y»—(n—2)}; the symbol z, will be used in this 
paper. We except from simple theory that >, will 
be zero on average, positive under adverse conditions 
and negative under favourable conditions. It will 
be shown later that, in practice, the average value 
of =, for a given day tends to be positive. 


4. Comparison of Calm and Stormy Days 


A period of six months, from July to December 
1958 inclusive, was selected for study. According 
to the B.B.C. reports, severe storms occurred on 15 
days in this period and there were 49 days with 
moderate storms. 

It was decided to compare values of =, on days of 
severe storms with values on samples of calm days. 
Saturdays and Sundays were not used in either 
sample because computer results for tasks on these 
days had not been calculated. There then remained 
12 days of severe storms for study. Four calm days 
were chosen, one each in July, September, October, 
and November, thus giving a reasonable distribution 
throughout the period. ays on which magnetic 
storms or Dellinger fades occurred were not saute 
for this sample, but otherwise the days were chosen 
at random. 
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The results are set out in table 1 below.: The TABLE 3. =, for moderate and severe phases of storms 
average value of 2, for all stormy days is +1.27 | ==> SS a 
and is significantly greater than the value of +0.57 ‘ Average value of En 

* > O58 
on calm days. ; = sears acest ate 

It is clear from table 1 that propagation conditions — — 
were worse from the point of view of DF accuracy | —— ~~ EE : Pe et - 
on days of severe storms than on calm days, as | isth August... +1 414 
expected. The question next arises as to whether this | {3th September__.-----------.-.-- eons tt +1 
effect applies to all frequencies, or is more marked in | 28th October__--_------— = ha. Eee 

° bs Z oe mm 6 ‘ 4/5th December- ----_- Sees ae ee +1.4 Sei orks 2) 
a particular part of the HF spectrum. Table 2 shows | 17/isth December_.__ ae pa tean 121 Eanes 
the value of =, for 1 Me/s intervals for the same sii 
sample of calm and stormy days. The difference | 9h January. Seihatestcietveeiea +1.0 +2.9 
. * d > 28th February/Ist Mareh_____.__.___- . +0.9 | ae eee 
between the results is more marked on the first and | eo es Sst Sera eee Sea = . +1.8 
second frequency bands (i.e., up to 10 Me/s) than | 3340 MMOs p----------- 2 ey - 
on the third band, but all are affected to some extent. ; a 
° : 5 4 Comilined resmits :. 2... .2.-2....- es +1.17 +1. 66 

In order to investigate further the effect of storms sai 
on accuracy, it was decided to use Slough measure- 
ments of fo F2 to define periods of moderate and The storm on 8th and 9th July 1958 had to be 
severe storms as described earlier. The period con- | omitted because the relevant logs had been destroyed. 
sidered was extended to June 1959. Table 3 shows | The average value of =, shows the expected in- 
values of 2, during the moderate and severe phases | crease under increasingly stormy conditions. 
of the storms considered. It should be noted that several dispersion factors 

were affected by the rejection of bearings showing 
ar Mar Tal = ah ‘4; aT hd achi >» 
TaBLE 1. Comparison of calm days with days of severe storms large de partures from : the BPE. Eithe r mac hine 
coe Ai TE Ae rejection or operator rejection may be applied to the 
‘iat . ' results calculated by the computer; the criteria for 
alm days Stormy days . . : ° ° . 
eal ee eS ta Jo Omyes si rejection will not be described in detail here, but 
cea i | depend in both cases on an examination of the dis- 
Average Average | ; r 2 i - : 
1958 value of | 1958 valueof | persion factor. Rejection of a bearing which is 
aoe Te er Sn: obviously ‘“‘wild” normally leads to a considerable 
eee ree, per oe reduction in this factor. 
22d September.-------.--| +0.6 | 9th July-———-- 2-2. +12 An examination of computer results showed that 
PO 0.7 BF 2 ° ‘ * 
ser seat: le apgamat is (ao a higher percentage of bearings were rejected on 
—_—_—_——_——— 5th September. -..__.-..-- +1.1 : . r dave 5 » ave ave of sever 
pe ee ee pe Ae | tn stormy days than on calm days. On days of severe 
| 24th October-._.....----_- +1.6 storms about 3.9 percent of all bearings were re- 
| 28th October.............- +1.4 . ther mis ee “ 
4th December _...........| +410 jected, whereas on the calm days examined only 
5th December. -___- : +1.3 9 re rare Treiee Mise rac ‘< ‘ —_ 
(Re: osama i 1.2 percent were rejected. This result is to be ex 
| 18th December -------.-- +1.3 pected in view of the increased size of the dispersion 
oS es | 41.27. | factors on stormy days. The greater percentage of 
ee =i aes |__| rejections on stormy days tends to reduce the 
difference between the average value of 2, on calm 
TaBLE 2. Dependence of En on and stormy days; hence the difference would be 
frequency ! even more marked than in table 1 if the same per- 
I dn eee centage of bearings had been rejected on both types 
Calm days | Stormy days of dav. 
Frequency : eS ‘ . . : 
A sort was carried out to decide whether bearings 
| a taken by any particular DF operators gave rise on 
ei silent hethcnninaeneriin average to abnormally large or small contributions 
Mels to the dispersion factors. Results were inconclusive. 
he See 2.7 
2-3 $0.3 +1.9 
3-4 1.5 2.0 . 
+5 +12 rs 5. Dependence on Mode of Propagation 
5-6 +0. 2 | +1.6 
i-6 4.9 | +172 Individual tasks on two calm and two stormy days 
6-7 +3 | +65 were examined with a view to determining the prob- 
— 7 7; able mode of propagation. 
8-9 4:2 +1.3 — 
10 | 0 +14 The range from the BPE for the transmitter to 
610 826]. 87 +1.41 London was taken to be a reasonable approximation 
ae a a oe to the true distance of the transmitter from the DF 
. < ° ° 1h) 
11-12 ae ae station. The maximum usable frequency (MUF) 
a - 3 . . , 16 
13-14 "2 Day for the same path was derived for the /, F'1, and F2 
= Pry rer layers from the U.S. Central Radio Propagation 
ree —— er Laboratory (CRPL) predictions. These frequencies 
ie were then compared with the known signal fre- 
580999613 227 








quency, and possible modes of propagation noted. 
By suitable division of the path it was possible to 
decide whether multiple-hop modes were likely. 

The tasks were then sorted into two groups de- 
pending on whether or not propagation was via the 
F2 layer. The sort was necessarily approximate, 
since neither the predicted frequencies nor the dis- 
tances used were free from error. However, the 
assignment of mode of propagation was expected to 
be accurate enough to reveal any large difference in 
DF accuracy dependent on the layer involved. 

Results of the sort are shown in table 4. =, was 
again used as a measure of the scatter of the bear- 
ings. From the table it appears that direction 
finding on signals propagated via the # or Fl layers 
is less affected by ionospheric storms than direction 
finding on signals propagated via the F2 layer. 


TaBLeE 4. Dependence of =, on probable mode of propagation 








| Average values of 2, 





| Two calm days | Twostormy days 








+0.5 
| + .6 +1.5 





6. Effect of Storms on Systematic Errors 


The errors of bearings taken on targets of known 
position were calculated and averaged for each DF 
station for three calm days and three stormy days. 
On the calm days, the average systematic error was 
+0.6°; on the stormy days it was +0.1°. 








The’nature of systematic errors is not completely 
understood, and it is therefore not possible to draw 
any conclusions from the reduction when conditions 
are stormy. Each systematic error probably rep- 
resents a combination of effects due to instrumental 
errors, site errors not removed by calibration and 
ionospheric effects. 


7. Conclusions 


The investigation has established that DF accuracy 
deteriorates during ionospheric storms and is worst 
during severe storms. The effect of a storm may 
extend over the whole HF band, but is more marked 
below 10 Me/s than above. The main effect seems 
to be on signals propagated via the F2 layer. 
Systematic errors during storms generally become 
less positive. 

The results imply that the Brooke variances should 
be increased during storm periods; such an increase 
would lead to rather larger probability rectangles on 
fixes taken in these periods. 

The present investigation demonstrates the 
advantages of adopting an objective classification 
system for bearings and of using an electronic com- 
puter to reduce DF data on a uniform basis and in a 
form suitable for analysis. 


The author thanks Miss A. Partridge and Mr. R. 
M. Rampling who carried out most of the analysis 
and who made many valuable suggestions. 

This paper is published by permission of the 
Director, Government Communications Head- 
quarters, Cheltenham, England. 
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Vol. 65D, No. 3, May-June 1961 


Phase Difference Observations at Spaced Aerials and 
Their Application to Direction Finding "* 


W. C. Bain 


(August 22, 1960; revised December 15, 1960) 


The measurement of phase differences at two aerial pairs, each spaced by three to four 
wavelengths, was used to give mean bearings over periods of three to five minutes on certain 


transmitters with frequencies near 6 Mc/s. 


The results showed that site errors had prob- 


ably been reduced by this system to less than 0.5 deg, and possibly to 0.2 deg, which is a 


value derived on theoretical grounds. 


deviation was found to be slightly greater than anticipated. 


The variance of the mean bearings due to lateral 


However, the results were 


such as to indicate that the performance of a wide-aperture direction finder should not fall 
seriously below that which would be expected of it theoretically on the basis of existing 


knowledge. 


1. Introduction 


In the winter of 1954-55 a number of observations 
were carried out at the Radio Research Station in 
which phase differences were measured between the 
signals picked up at aerials spaced by several wave- 
lengths. Recently it was noted that these results 
were relevant to the problem of determining the 
accuracy obtainable from a wide-aperture direction 
finder, an instrument capable of improved perform- 
ance in conditions of multimode propagation |Bain, 
1956]. It therefore seemed desirable to present this 
account of the work and its implications here. 

The four aerials used were short unipoles; they 
were placed on the circumference of a circle of 
diameter 180 m so that one pair of aerials was alined 
in the north-south direction and another in the east- 
west direction. Inside a hut at the center of the circle 
each of these pairs of aerials was connected to phase- 
measuring equipment of the type described by Ross, 
Bramley, and Ashwell [1951]. When the receivers 
were tuned to the desired station, the overall phase 
balance of the system was checked, and the phase 
differences between the signals picked up on the 
north-south and east-west pairs were recorded 
photographically at three-second intervals for a 
period not exceeding five minutes, thus giving about 
70 values for each. Ionospheric conditions were 
characteristic of the time in question and were never 
greatly disturbed. Observations were confined to 
the hours of daylight. 

It is important to note that such a system is not 
directly suitable for use as an operational direction 
finder. Not only does it have bearing ambiguities, 
but the bearing calculated from the phase differences 
is not in general correct when two rays are incident 
on it from the same azimuth. 


1 Contribution from Department of Scientific and Industrial Research, Radio 
Research Station, Ditton Park, Slough, Bucks., England. 

2 Paper presented at the Conference on Transmission Problems Related to 
High-Frequency Direction Finding, at UCLA, June 21-24, 1960. 





2. Results 
2.1. Method of Analysis 


In reducing the observations a mean yw of each 
distribution was calculated in the following way. 
If there are n phase differences y% in each distri- 
bution, then 


y= me{ © exp (jy) }. 


This avoids the ambiguities which can arise if an 
ordinary arithmetic mean is used. 

The standard deviation (c) of the distribution was 
taken to be the root-mean-square departure from 
¥. For a completely uniform distribution of phase 
differences the standard deviation would have the 
value 104° approximately. <A finite number of ran- 
dom phase differences will, of course, give a standard 
deviation somewhat less than this, probably about 
94° for the numbers of observations used here. 


2.2. Observed Values 


The observations were concentrated upon a few 
stations which had been found from a previous 
experiment to give widely spread phase difference 
distributions. Of these, the most frequently used 
were the following. The distances and bearings 
shown are from the site of the equipment at Wink- 
field, near Windsor (Berkshire). 

It will be seen that with the aerial spacing used 
here, the diameter of the system will be in the range 
of three to four wavelengths for reception of these 
stations. 

It was indeed found that very large standard 
deviations were occasionally obtained; a few were 
so great that they could have arisen from a random 
distribution of phase differences. A camplete set of 
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the standard deviations obtained on these three sta- 
tions is shown in figure 1. These can be seen to be 
well spread out over the range of possible values. 
However, the north-scuth and east-west standard 
deviations are correlated to some extent, for if one 
is large the other also terds to be large. In the case 
of Arganda, the north-south phase differences were 
much more widely spread than the east-west; the 
mean standard deviations were in the ratio of 2.0 
to 1. This probably happened because the Arganda 
bearing is only 10° away from the north-south direc- 
tion, combined with the fact that the incident 
radiation was spread much more in elevation than 
in bearing. 


3. Bearings Deduced From Phase Measure- 
ments 


If the phase difference between north and south 
aerials is ¢,, and that between east and west aerials 
is ¢:, then the bearing y can be taken to be derived 
from the formula 

tan y= ¢»/¢). 


With widely spaced aerials the values of ¢, and 
¢. obtained from the measurements may differ from 
the true values by some integral multiple of 27. 
However, as the position of the transmitter was 
known in these experiments the ambiguities could 
be eliminated by choosing the multiples of 27 so as 
to bring the bearing as nearly as possible to the true 
bearing of the transmitter. The multiple of 27 
which had to be used never exceeded 3 in any of the 
present instances. 

Mean bearings were calculated by using the mean 
phase differences y from each of the groups of obser- 
vations; the results obtained are summarized in table 
2, which gives results for the mean of all mean 
bearings for each station. 

It can be seen that the results are consistent with 
the hypothesis that there is no systematic error in 
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Figure 1. Histogram of the number of coses of standard 
deviation of the phase difference distributions with values 
as shown. 





the observations. However, it is useful to consider 
in some detail the possible errors which might arise. 
They can be divided into three main categories for 
the present purpose: 

1. Rapidly fluctuating errors, with periods of the 
order of seconds; they may be due to wave inter- 
ference or polarization errors. 

2. Lateral deviation errors, with periods of the 
order of tens of minutes, due to the effect of iono- 
spheric tilts. 

3. Site errors, which will be considered here to be 
virtually constant with time. 

The first two of these errors would be expected to 
give rise to a spread of values in the mean bearings 
obtained for each transmitting station, but the 
third would lead to a systematic error. On the 
assumption that the site error due to a narrow- 
aperture Adcock direction finder is 2.5°, it can be 
estimated from the curves for site-error reduction 
given by Hopkins and Horner [1949] that the site 
error for each of the stations in this wide-aperture 
direction finder would be of the order of 0.2°. Now 
the range in the fiducial limits for the mean as 
given in table 2 is sufficiently great to conceal alto- 
gether an effect of this magnitude. It is still, of 
course, possible that the site error of the system is 
indeed as small as this, but the most that can be said 
from the experimental evidence here is that it is 
probably less than 0.5°. This is certainly a worth- 
while improvement on a narrow-aperture direction 
finder. It should be borne in mind that this result 
was achieved despite the large fluctuations in phase 
difference commonly found with the transmitters 
listed in table 1. 

The magnitude of the spread of the mean bearings 
can also be examined to see if it is in accord with 
what would be expected. To do this some estimate 
must be obtained of the size of the errors due to 
rapid fluctuations and to lateral deviation. The 
error variance (¢3,) due to rapid fluctuations can be 
assessed from a knowledge of the variance of the 
north-south and east-west phase differences (of and 











TABLE 1. Details of stations observed 

Location | Frequency | Distance | Bearing 
os ee eee } tat cell 

| | 
Me's | km | deg 

Allouis, France -- cies sabes Se 6.20 | 520 155.0 
PIED, BOR i. sacs anna nade cksncvcace} 5.99 1270 191.0 
Schwarzenburg, Switzerland___.__.._---- 6.16 | 730 128.1 














TABLE 2. Errors in the consolidated mean bearing 
No. of groups | Bearing error | 90% fiducial 
Station of observations consolidated | limits for con- 
mean solidated mean 
ienhomnashesitangin — 
deg | deg 
EA CORT: 13 | -0.1| 0.8 to +0.6 
LS SERRE TS 5 | —.4) — .8 to .0 
Schwarzenburg..-.-.--..---- 13 | +.2 — .3to+.7 
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o3, say) in each group. If there are n, north-south 
phase differences and n,. east-west phase differences, 
and @ is the bearing of the transmitter, then a 
reasonable estimate of o%, is given, as is shown in 
the appendix, by 


3 —(228 *y 4( 
T 28 nN \ 2B/ m 


sin @\? o? 








where 


p= cos 6, 


d=the diameter of the system, 
\=the wavelength, and 
5=the angle of elevation of the incident wave. 


Data about the magnitude of ionospheric tilts in 
the F region can be used to give the variance due to 
lateral deviation. The figure of 1.2° determined by 
Bramley [1953] has been adopted here as the standard 
deviation of the magnitude of ionospheric tilts in one 
plane. The expected lateral deviation variance for 
each station was found from this, again with the 
assumption that the reflection was from the F2 layer. 

In table 3 the variances of the individual mean 
bearings as obtained in two different ways are com- 
pared; if the estimates made are correct the figures 
in the last two columns should not be significantly 
different. 

The number of observations on Arganda is too 
small to give much hope of establishing a significant 
difference between variances, and in fact the values 
shown in the last two columns for this station are not 
significantly different. However, for both Allouis 


and Schwarzenburg it looks as if the estimated vari- | 


ance may be too small, for in each case the difference 
between observed and estimated variances is nearly 
significant at the 95 percent level. As the rapid- 
fluctuation variance has been estimated from internal 
evidence, it looks as if this difference is associated 
with the lateral-deviation component. However, 
for present purposes it is sufficient to note that the 
observed values are not very far from the estimated 
ones, and there is nothing about these results to 
suggest that our ideas are seriously wrong about how 
to apply known information about the ionosphere 
and direction finding to predict the performance of a 
wide-aperture direction finder. 

















TaBLE 3. Estimated and observed variances 
= ; — - 
| | Estimated 
| Estimated Estimated variance of Observed 
: | variance due | variance due | individual variance of 
Transmitter | to rapid flue-| to lateral! mean bear- | individual 
| tuations,o%, | deviation | ings (sum of | mean bear- 
| yh previous 2 | ings 
| | columns) | 
ee Ae eT ee |e (RE A 8 (ee Aa] eI 
deg? | deq deg? deg? 
Oe 0. 08 1.17 1.25 2.05 
es 10 | 0. 32 0. 42 0.14 
Schwarzenburg. __._- .07 | . 54 61 | 1.08 





| 
| 
| 





4. Conclusions 


The mean bearing indicated by the phase-difference 
measuring system considered here, which is of three 
to four wavelengths in aperture, is likely to be 


' accurate to +0.5°, and may be as good as 0.2°. 


To achieve this accuracy time-averaging over a few 
minutes had to be applied, but in an operational wide- 
aperture direction finder the fast bearing fluctuations 
should be smaller than was the case here and should 
require less time averaging to reduce them to an 
acceptable level. In the practical direction finder the 
time required to reduce the variance due to lateral 
deviation of the rays will probably be the same as 
here, and several hours would be needed to average 
out its effects. In the present experiments the vari- 
ance ascribed to lateral deviation turned out to be 
somewhat larger than expected, but the difference 
was not so great as to suggest that large unantici- 
pated effects were present. There are, therefore, 
grounds for hope that wide-aperture direction finders 
will be capable of the performance expected of them 
on the basis of existing knowledge. 


This work was carried out as part of the program 
of the Radio Research Board of the Department of 


i Scientific and Industrial Research, and this note is 


published by permission of the Director of Radio 
Research. 
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6. Appendix. Error Variance in the Mean 
Bearing of a Group of Observations Due to 
Rapid Fluctuations in the Phase Readings 


The true bearing here will be taken to be given by 
the expression 
y=are tan (¢,/¢)), 


where ¢,=26 cos y, ¢:=28 sin y. 
The error in bearing 6y is then given by 


_cos yé¢2.—sin voor 


oy 2B 
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The north-south phase difference, which is here 
(¢,+6¢,), has a variance of given by 


oj ai (691), 


using the bar to denote an average. If there are 
n, observations of this phase difference, the variance 
of the mean is given by (6¢,)?/n;. Similarly, the 
variance of the mean of the east-west phase distri- 
bution is (6¢2)?/nz. The mean value of the bearing 


error is given by 











— cos yéd;—sin yoo, 
y= ° 


28 





Hence the variance of the mean bearing is given by 


cos y\? o3 , (sin vy\" o? 
oie=( + F 








28 Ny 28 Ny 

This derivation applies only when 6¢,; and 6¢, are 
uncorrelated. However the formula has been checked 
by calculating bearings for each individual pair 
of phase differences for certain groups; the standard 
deviation obtained has been in each case quite close 
to the value given by the above formula. 
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Research at the National Bureau of Standards Appli- 
cable to Long-Distance Location and Direction-Finding 


Problems” 


R. Silberstein 


(September 8, 1960; revised December 30, 1960) 


Research pertaining to radio-location and direction-finding problems conducted at the 
National Bureau of Standards since 1941 includes evolution of a technique for determin- 
ing polarization error, study of non-great-circle bearings, the development of a rapid- 
scanning directional antenna, and the development of Loran-C with its precision timing 


capabilities. 


1. Introduction 


Work in long-distance radio-location and direction 
finding at the National Bureau of Standards began 
in the Radio Section in Washington, D.C., which 
became the Interservice Radio Propagation Labora- 
tory (IRPL) during World War IT and later became 
the Central Radio Propagation Laboratory (CRPL) 
now located at Boulder, Colorado. 


2. Polarization Error Studies During World 
War II 


A direction finder may be responsive to vertical 
or horizontal polarization or to both polarizations 
phased in various ways, but one polarization may 
produce an undesired response which causes an error 
in the reading. This undesired response is a func- 
tion of the geometrical properties of the instrument, 
its mountings and various accessories, such as RF 
cables. A combination of modes in the downcoming 
ionospheric wave results in more or less random 
polarization. From instant to instant the relative 
amplitudes and phases of the wanted and unwanted 
components of the arriving wave keep changing 
causing the polarization error to change. 

It would be a very difficult thing to compare direc- 
tion finders and evaluate their polarization errors by 
observing sky waves, for which bearings depart from 
their true great circle values for a variety of causes. 
Some simple test technique was desired which would 
offer an estimate of the polarization error propensi- 
ties of any particular direction-finder design. The 
first test of this type had been evolved by Barfield 
in England. A balloon with a target transmitter 
and antenna producing a circularly-polarized wave 
at a downcoming angle of 45° was placed near the 
equipment under test. The polarization error ob- 
served was called the “standard wave error.” 





! Contribution from Central Radio Propagation Laboratory, National Bureau 
of Standards, Boulder, Colo. 

* Paper presented at the Conference on Transmission Problems Related to 
High-Frequency Direction Finding, at UCLA, June 21-24, 1960. 





At NBS, a group under the late Harry Diamond 
devised techniques whereby pickup factors for the 
desired and undesired polarizations were determined 
from fixed target transmitters. These pickup fac- 
tors were then treated analytically in accordance 
with a procedure designed for each type of direction 
finder and polarization errors determined for any set 
of conditions. 

In the case of the simple Adcock antenna, the 
desired polarization is parallel to the plane of prop- 
agation (vertical when the downcoming aagle is zero) 
and the undesired polarization is perpendicular to 
the plane of propagation (horizontal). 

The pickup factor for a wave polarized parallel to 
the plane of propagation was called h. That for 
a wave polarized perpendicularly to the plane of 
propagation was called k. If h and k were known 
then the polarization error could be written (in the 
notation used by the Diamond group) as: 

Fi 


tan a. 


h E\. cos y 


where 


e is the polarization error angle, 

FE, is the resultant total field vector perpendicular 
to the plane of propagation (i.e., horizontal 
polarization) after combination of the incident and 
ground reflected wave, 

EF}, is the resultant total field vector normal to the 
ground (vertically polarized) after the combination 
of the ground wave and ground-reflected wave 
polarized parallel to the plane of propagation 

yY downcoming angle 


A figure of merit was proposed for the case of 
equal resultant fields with the direction finder on 
the ground. Here tan ¢,=(k/h), either k/h or e, 
being the figure of merit. 

The advantage of the method was that A and k 
could be measured on the ground. One type of 
error called ‘radiator parallax” occurred in measur- 
ing k in direction finders like the balanced H Adcock. 
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This was due to a net unwanted vertical pickup in 
the system when a horizontal dipole was used as 
a source of perpendicular polarization. In the case 
of spaced vertical loops another type of error called 
“collector parallax” was studied. The forward tilt 
component of the target transmitter’s magnetic 
vector produced an unbalanced voltage in the 
system. Each type of error could be minimized by 
increasing the distance to the target transmitter. 
The tests for determining fh and k were actually 
conducted with both the direction finder and the 
target transmitter on wooden platforms above 
ground. The elevated platforms were especially 
useful in obtaining a strong horizontally-polarized 
ground wave for measuring k. The importance of 
this arrangement is attested to by the author, who 
at one time attempted to measure the horizontal 
field strength over a stretch of swamp where vertical 
polarization was attenuated about 70 db less than 
horizontal. The transmitting antenna was a wire 
dipole. Minor unbalances of this antenna as it 
swayed in the breeze produced great surges of output 
in the field strength meter as incidental pickup of 
the vertical polarization responded to the vertical 
component of the unbalanced antenna system. 


3. Radio Propagation Effects on High Fre- 
quency Direction Finders (World War II 
Period) 


In this project several receiving stations took 
routine bearings on distant stations. Results were 
compared with predicted path MUF’s and field 
strengths. Much attention was paid to an old 
phenomenon which was noted by the British in the 
1930’s [Keen, 1947] and which has recently been 
attracting attention again, namely, propagation of 
signals by a quasi-scattered non-great-circle mode 
after failure of the regular path MUF. 

The British results were duplicated again in various 
ways during World War II. For instance, when 
a BBC broadcast beamed at North Africa was 
monitored in Virginia, the bearings taken during 
the hours of predicted MUF above the operating 
frequency indicated that the bulk of the received 
energy was arriving via the great circle path. After 
predicted MUF failure of the eastern control point, 
however, the signals exhibited reduced fluctuating 
field strength with widely swinging bearings suggest- 
ing a scattering region in Africa. It was assumed 
that the MUF was still high enough to the south 
of England for propagation to the scattering region 
from which energy was being received in Virginia 
by regular propagation over a more southerly path. 
In those days it was assumed, following T. L. 
Eckersley, that high- frequency backscattering and 
sidescatteriag occurred in the FE region rather than 
at the ground where most of it is now known to 
take place, but the general assumption of propagation 
from a scattering region was valid. 

Frequently European stations in the lower part 
of the high-frequency band, from which great-circle 
propagation via the regular layers was not possible, 








were observed to have a northerly bearing deviation 
when field strength was weak, suggesting a mode 
involving an Auroral Zone gradient. These studies 
are called to mind by the results of recent Canadian 
experiments on 40 Mc/s BBC television in which 
both types of bearing deviations were seen. 

Similar observations of signals from United States 
stations were readily explained, but a group of 
widely deviated bearings of a Singapore station, 
observed at Sterling, Va., were never fully understood. 


4. Recent Developments in Techniques 


A precise system of direction finding, developed 
for use on spherics at VLF, has potentialities for 
other uses [Hefley, Linfield, and Davis, 1961]. The 
name “Ephi’” was chosen for the system because the 
bearing of the transient signal is determined from 
the relative phase (phi) of the vertical electric field 
(E) received at spaced antennas. In the system 
three vertical antennas are erected at the corners of 
an equilateral triangle of sides 1/5 to 1/10 wave- 
length. A transient signal from a given direction 
arrives at each antenna at a different time. These 
arrival times are not compared directly; instead, 
sets of fixed time delays are inserted in each channel, 
each set defining an azimuthal sector. Whenever a 
spheric occurs in such a sector the pulses arrive 
nearly simultaneously at the outputs of their re- 
spective channels and trigger a coincidence circuit. 
Thus the number of spherics in any sector may be 
counted. The system is superior to the conven- 
tional crossed-loop system in that it works on only 
vertical polarization. Also it lends itself to auto- 
matic readout much more readily than a system 
employing oscilloscopes. 

Testing of a new type of wide-aperture direction 
finder using electronic scanning has just been com- 
pleted at CRPL [Cottony and Wilson, 1961]. The 
model tested had a beamwidth of 614° and scanned 
over an are of 45°. Although the present operating 
frequency is near 40 Me/s special application of the 
basic features to high frequency work appears to be 
feasible. 

A theoretical paper was written on a proposed 
experimental technique for measuring the angle of 
arrival, azimuth and polarization of a downcoming 
wave by the use of a combined 4-element Adcock 
and crossed-loop direction finder [Wait, 1959]. 


5. Loran C as a Source of Precise 
Timing 

A system of Loran stations on 100 ke/s, known as 
the East Coast chain has recently been installed with 
a master station at Cape Fear, N.C., and two slave 
stations at Martha’s Vineyard, Mass., and Jupiter 
Inlet, Fla. [Doherty, Hefley, and Linfield, 1961]. 
These stations employ a recently declassified system 
of navigation formerly known as Cytac. 

In Loran C the difference between the differentiated 
signal envelope and the signal itself yields a wave 
crossing the zero axis at a point so well-defined that 
it can be used to eliminate ambiguity of a single cycle 
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in timing a point on the wave. Synchronization is 
accomplished entirely by computer circuits. 

A slave station can be held to 0.1 usec or better. 
The timing accuracy using the Loran C ground-wave 
emissions 1s 1 ysec or better over land for ranges to 
about 1,500 miles and 0.1 usee over sea for ranges to 
about 2,000 miles. Much greater distances can be 
covered by sky-wave synchronization with errors of 
only a few microseconds. Because of the sporadic 
impulsive nature of atmospheric noise and the fact 
that the Loran C pulses arrive at a 20-cycle per second 
rate, it is possible to synchronize oscillators with a 40 
db ratio of average noise power to signal power. 

The advantage of a positioning system with an 
accuracy 1,000 times greater than that for independ- 
ent clocks compared with WWV sky-wave signals 
(assuming a 1-min integration time) is obvious. 
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Design for Spinning Goniometer Automatic Direction 


Finding 


1, 2,3 


W. J. Lindsay and D. S. Heim 


(September 29 


This paper discusses some aspects of ins 


goniometer radio direction finder essentially automatic in operation. 


, 1960) 


trumentation design in making a spinning 
Also discussed is the 


application of narrow band synchronous post detection filtering for improving the bearing 


sensitivity. 
1. Introduction 


In this paper various aspects of the problem of 
automating a spinning goniometer radio direction 
finder (hereafter abbreviated SGDF) are discussed. 
These considerations involve the handling of the 
directional information and do not concern other 
control operations such as tuning, switching, etc. 
Such a system is desirable for use as a research tool 
and to study the feasibility of applying various con- 
cepts operationally. Although the work was done on 
a AN/TRD-4 unit, the concepts may be applied to 
any SGDF. The availability of automatic operation 
does not preclude the operation of the equipment 
either partially or entirely in the usual manual 
fashion. 

In a narrow aperture SGDF the simulation of a 
spun antenna is achieved by balanced amplitude 
modulation of the incoming wave at a rate set by the 
spin frequency of the goniometer. If sense voltage 
(from an omnidirectional antenna) is added to this, 
the effect is to restore the carrier component and we 
have a conventional amplitude modulated signal. 
These two modes of operation will be referred to as 
the DSB and the AM cases respectively. The AM 
wave may be expressed as: 


e(t)=E (1+k cos wt ei (1) 
Taking the real part: 


, rE 
Rele(t)|=E cos ee [cos (wou) t+ cos (wy—a) t] 
(2) 
The detected envelope of the AM wave is; 


ég=(— (1+k cos af)] (3) 





While the detected envelope of the DSB wave 


' Contribution from Cooley Electronics Laboratory, Electrical Engineering 
Dept., University of Michigan, Ann Arbor, Mich. This work was supported 
by the U. 8S. Army Signal Research and Development Laboratory under Con- 
tract No. DA 36-039 se 78283. 

> The presentation of this paper was sponsored by the Office of Naval Research. 
Reproduction in part or in whole is permitted for any use of the United States 
Government. 

3 Paper presented at the Conference on Transmission Problems Related to 
High-Frequency Direction Finding, at UCLA, June 21-24, 1960. 
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(AM wave minus the carrier) is: 


Kh 


Tv 


yi, 


- 


7 
ae 


In the above expressions w is the modulation fre- 
quency of the spinning goniometer while w is the 
center frequency of the incoming radio frequency 
wave. It will be noted that in (2) the sideband 
frequencies are displaced from center frequency only 
by the spin rate of the goniometer, which in the 
present case is 30 c/s; also the fact that in the detected 
envelope of the DSB wave only even harmonies of 
the modulation frequency appear, there is no funda- 
mental component. The detected envelopes of both 
the DSB and the AM waves are shown in figure 1. 


cos 2 t— 3 5 


€a 
3.0 


cos 4a,f—. . | (4) 


— ONE CYCLE OF ROTATION 





OSB WAVE 


(a) 





AM WAVE 
(b) 


Waveforms for detected DSB and AM mode signals. 


FIGURE 1. 








In this type of system the bearings are determined 
by the position of the null of the antenna pattern 
with respect to a point of reference, usually North. 
It is seen that in figure 1a there are two such points, 
thus the familiar 180° ambiguity, while in figure 1b 
there is only one point at which the signal goes 
through the null in a specified direction of slope of 
the wave. It should be noted that the null point 
has been shifted by decoupling to the dotted line. 
This brings about an effect, which at present is the 
most significant detriment to the bearing accuracy 
obtained, of a shift in the apparent bearing when 
the character of the receiver is altered even slightly. 
For example when the gain control is changed over 
its range the bearing will change by two or three 
degrees. Since, as explained in a subsequent section, 
the reconstruction of the visual display and the zero 
crossing comparison for measuring relative phase 
both depend on having a constant average output 
from the receiver for a given input signal form but 
at different levels, it is evident that a receiver is 
needed that has low envelope distortion in order to 
eliminate this effect. There is indeed a great need 
for receivers designed for instrumentation purposes 
as well as communications, particularly with regard 
to the phase character of the receiver. If the signal 
is pulsed then the problem of maintaining the average 
is more difficult because d—c components are intro- 
duced by the pulsing action of the signal. These 
will average out over a period of time, but for short 
periods cause havoc in the bearing determination. 
In viewing the DSB and AM situations certain 
things come to mind with regard to performing (1) 
Instrumental measurement of the bearing angle, and 
(2) Improvement of the signal to noise ratio by 
applying narrow band synchronous post detection 
filtering. 

In the DSB case the null point region is that of 
the lowest signal to noise ratio in the entire cycle of 
rotation; thus in a noisy situation, which is the usual 
case rather than the exception, the null position is 
imperfectly defined, and indeed the whole structure 
of the wave, that is the collection of harmonies which 
define the structure, is needed in order to determine 
just where the average null point is located. A 
similar effect may be had by filtering out the second 
harmonic term in (4) and noting the relative phase 
of this wave to a reference wave. 

In fact an increase in signal to noise ratio compara- 
ble to the AM case may be achieved but since an 
additional ambiguity is introduced this is not as 
attractive as using automatic sense if one wants to 
use narrow band filtering. If it is desired to keep 
the DSB information as is, except to filter the noise 
from it, then a comb filter may be used. This may 
be either the discrete frequency type or more likely 
the commutated type. In either case it is consider- 
ably more involved than a simple narrow-band single 
frequency filter. In the AM case the phase measure- 
ment is made on a sine wave (for a continuous signal) 
which has the same fundamental period as the 
goniometer rotation therefore the ambiguity is 
removed, and one may now apply filtering at the 


| amplitude. It 





basic scan frequency rather than at twice the scan 
frequency. A description of the means of obtaining 
automatic sense, the design of the bearing readout 
system and a discussion of narrow band post detec- 
tion filtering with some operational illustrations 


follow. An important aspect of the design is mini- 
mum cost with reasonable performance and 
simplicity. 


2. Automatic Sense 


The addition of a carrier or sense signal to the 


| goniometer output changes the receiver input to a 


standard amplitude-modulated signal, modulated 
at the spin rate. The phase of this modulation 
relative to the phase of the rotation gives an unam- 
biguous bearing. It can be shown that if the phase 
of the injected carrier is within plus or minus 90° 
of the correct value it distorts the envelope, but only 
to the extent of adding harmonics (which may be 
filtered out) and causing a decrease in fundamental 
does not affect the relative time 
phase of the fundamental. Thus the generation of 
the sense signal is not so critical as is often supposed. 
If the sense signal is larger than the directional 
signal and is held to within + 45° of the proper value, 
adequate performance may be obtained. 

The simplest means for obtaining the sense signal 
is to insert an antenna in the center of the DF array. 
Since the goniometer output has suffered some loss, 
as well as being the difference between two signals, 
the injected sense voltage can easily be made suffi- 
ciently large. If the distance across the array is small 
compared to a wavelength, a fairly simple network 
will maintain the phase margin within the required 
limits [Keen, 1947]. In fact, with rather heavy 
filtering, satisfactory operation has been obtained 
in the case illustrated with no additional phasing 
over the band of operation. The R-390 receiver is 
conveniently designed with both a balanced and 
unbalanced input so that the sense signal may be fed 
directly into the unbalanced antenna input jack. 

The instrumentation for producing the proper 
pattern on the screen of the CRT is depicted in 













































figure 2. This unit serves two functions: (1) it 
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Figure 2. Automatic sense circuit. 
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full-wave-rectifies the sine wave output of the re- 
ceiver to give the usual propeller-type display, (2) it 
causes a blanking signal to be applied to the grid of 
the CRT in such a manner that only one half of the 
propeller appears on the screen. In the AN/TRD-4 
equipment the display is filled in with a 72 ke/s 
carrier. Figure 3 shows the usual input to the 
indicator goniometer during one complete revolution 
during DSB type operation. It can be seen that the 
propeller is described on the face of the tube twice 
during each revolution. During period 7, the tube 
must be blanked during the negative excursions of 
the 72 ke/s carrier. During period 7}, since the 
goniometer has rotated 180°, the blanking must occur 
on the positive excursions if the same half of the pro- 
peller is to be displayed. It can be seen from the 
circuit that with no signal from the receiver both 
triodes will conduct equally, and, since their inputs 
are 180° out of phase, cancellation will occur in the 
plate circuit. In the presence of a signal from the 
receiver (ideally a sine wave), the tubes will alter- 
nately conduct fully and be cut off. It is this 
alternating operation which causes the proper balf of 
the display to always appear as discussed above. 
The resultant 72 ke/s signal, when applied to the 
CRT grid, blanks during the negative excursions and 
brighteas or enhances the desired half of the pro- 
peller during the positive excursions. 


3. Bearing Readout Computer 


The basis for the computer design is that of zero- 
crossing comparison, as has been used by others 
[Florman and Tait, 1949; Hatch and Byatt, 1958]. 
Once the zero-crossing of the reference and DF signals 
have been determined they may be used to activate 
either an analog-type coincidence detector or a 
digital-type bearing computer. It was concluded 
that for visual monitoring purposes the analog was 
more suitable, and for detailed processing of the 
information the direct digital readout was more 
suitable. For the purpose at hand there was no need 
to eliminate the 180° ambiguity from the analog 
readout. 

The analog coincidence detector is fed into one 

channel of a two-channel Brush Recorder. The 
second channel is used to record signal strength. 
The particular Brush Recorder used has a frequency 
response up to about 50 c/s. It is desirable to have 
an adjustable readout rate in the digital system so 
that the bearing readout rate can be adjusted to fit 
the study being made. 
_ For studying the fast-fade variations, Bain [1955] 
indicates that in order to study the phenomena one 
should be able to take at least several bearings per 
second. Of course the rate for slow fluctuation 
studies should be much less than this, thus there 
exists the need for variable readout rate. 

Reference is made to figure 4, a block diagram of 
the readout system. The goniometer bearing scale 
was removed and machined out under the thumb rim, 
and 1° slits were cut on the outer surface of the rim 
with an indexing head. Thus, with a light in the 
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FicurE 3. Input to display goniometer during conventional 
DSB operation. 
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Ficure 4. Block diagram of automatic bearing computer and 


readout. 


wheel and a phototransistor on the outside, a tone 
wheel is available. A properly adjusted slit in front 
of the transistor is a necessity. The output of the 
phototransistor is fed into a two-stage transistor pre- 
amplifier and a pulse-forming circuit, then to the 
counter. 

The reference generator is a small permanent-mag- 
net-type alternator which was built on the end of the 
goniometer and is coupled directly to the shaft. 
The output from the reference generator is squared 
and formed into the start pulse in the zero-crossing 
detector. The output from the DF receiver is 
similarly treated in the zero-crossing detector unit. 
The squaring circuitry outputs from these units are 
fed to the analog coincidence detector and thence to 
the pen recorder for the visual record of bearing. 
The pulse output from the reference channel starts 
the counter, and the pulse output from the receiver 
stops the counter. The censor unit at present is to 
be used as a yes-no stop pulse control depending on 
signal level. The start pulse is delayed slightly in 
the gate generator in order that the same signal may 
be used as a reset control at North position. When 
the censor unit says “no” then the counter will 
cycle on a 360° count until the censor unit says ‘“‘yes”’. 
The reference point, North in this case, may easily 
be moved to some other point around the compass 
if desired. The problem here is that of averaging 
through the cut point. For example, if the bearings 
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were symmetrically distributed about North, than 
the average of these would be South, obviously the 
wrong answer. One simple solution to this is to 
have two cut points with two counts and pick the 
output that is farthest removed from a cut point. 
The counter output is then fed to the tape punch 
driver, which is controlled by the censored stop pulse, 
and thence punched onto paper tape. 
unit is to be transistorized the Tally 8 channel 
punch was chosen. Eight channels allow a sector 
coverage of 255°. In order to cover the full 360° 
sequential channel coding must be used which still 
allows a maximum of 30 readings per second. 


4. Narrow Band Synchronous Post Detection 
Filtering 


The application of NBSPD filtering has been 


Since this | 


discussed in the literature [Busignies and Dishal, | 


1949], however little seems to have been done to 
exploit its potential to improve bearing sensitivity. 
In fact, as Busignies and Dishal state, it appears 
that the capability of utilizing this is one of the 
fundamental advantages of this type of DF system 
over the co-operative navigational type system. 
This stems from the simple fact that in the great 


majority of situations the determination of the | 


average direction of arrival over period even as 
short as just a few seconds is a very narrow band 
process. From an informational point of view the 
channel capacity required is in terms of bandwidth 
the approximate frequency of variation of the bearing 
one wishes to observe. This is of course determined 
exactly by the reading time of the specific receiving 
system in use. This means that for most purposes 
a bandwidth of a few cycles (say less then 10) will 
be adequate, unless there is some wideband informa- 
tion in addition to the bearing information, or if 
one is interested in studying the instantaneous rapid 
fluttering effects caused by wave interference, or 
for other specialized types of signals. There is 
more to be gained from a given factor of bandwidth 
reduction in the predetection circuitry than at post 
detection. By the same token this is much more 
difficult to achieve. In applying post detection 
narrowing if BW (IF) is the bandwidth of the IF in 
the receiver and BW(A) is the modified audio 
bandwidth then the following relation gives the 
input carrier ratio for equal output in terms of 
bearing determination. This is done for purposes 
of comparison to the operational illustration shown 
later. 

BW(IF)]'“ 


Carrier Ratio=K 


and substituting for BW(IF)=2.0 ke 


BW(A) =5. 0¢/s 


9 1/4 
CR=K | —K (400)"~K (4.475) 








This means a difference in terms of db of: 


CR Diff (db) =20 log 10 (4.475) 
= 13.0 db 


Results from operation in this manner with auto- 
matic sense are given for purposes of illustration in 
figures 5 through 11. <A single tuned circuit with 
center frequency at the rotational frequency of the 
goniometer is placed after the detector. This unit 
has adjustable bandwidth with the narrowest po- 
sition at about 5.0 c/s. This is the condition under 
which the pictures were made. This gives a reading 
time of about 140 msec. which allows observation 
of bearing changes up to about 7.0 e/s. 

Figure 5 shows a strong stable signal in the broad- 
cast band, station WJR in Detroit. The double tip 
is due to imperfect foldover in the automatic sense 
circuit, which is being corrected. 

In figure 6, (a) and (b) are, respectively, the con- 
ventional TRD-4 display and the modified display, 
in which the target signal, at 6.93 Me/s, was adjusted 
in level until one could just estimate a bearing, at 
around 105°, on the conventional TRD-4 display. 
These shots are all 4-sec time exposures. 

In figure 7, the signal level is reduced by 3.6 db, 
at which point the TRD-4 display is just barely 
discernible as a bearing indication. 

In figure 8, the signal level is reduced from where 
it was in figure 6, by 9.6 db, which leaves nothing 
in figure 8(a), although it can be seen that there is 
still quite a good bearing indication in figure 8(b). 
The limit is reached by reducing the signal level by 
another 6.0 db, producing figure 9. One can still 
get a fair estimate of the bearing, and the pattern 
which resembles a one-sided halo, is developing into 





Ficure 5. Bearing indication on modified AN/TRD-4 for 
strong stable signal. 
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b 


Figure 6. Comparison of operational results from (a) conventional AN/TRD-4 and (b) modified AN/TRD-4. 


b 


FicurE 7. Comparison of operational results from (a) conventional AN/TRD-4 and (b) modified AN/TRD-4. 
Signal Level Reduced by 3.6 db from that in figure 6. 











FIGURE 8. 


b 


Comparison of operational resutts from (a) conventional AN/TRD-4 and (b) modified AN/TRD-4. 


Signal Level Reduced by 9.6 db from that in figure 6. 





Figure 9. Threshold condition for the modified AN/TRD-4. 


the characteristic background for this mode of opera- 
tion. With no signal there is a faint green back- 
ground level on the CRT face. When a signal is 
tuned in the propeller half suddenly rises out of the 
hazy background as though by intensity modulation. 
The difference between the input signal level for 
figure 6, and that of figure 9, is 15.6 db. 





Figure 10, depicts the appearance of a typical, 
fairly-weak FSK signal at 6.901 Me/s. 

Figure 11, shows the target signal at 6.93 Me/s. 
being pulsed at 20 c/s. The integration effect of the 
film favors 10a more so than the eye in practice. 


5. Summary 


It has been demonstrated that one may partially 
automate and otherwise improve the operation of a 
SGDF with respect to certain prescribed criteria with 
rather simple and inexpensive methods. One must 
of course keep in mind what the particular objectives 
for the system are and the resulting advantages to 
be gained by making certain changes and the price 
that sometimes must be paid to make those changes. 
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b 
FicureE 10. Typical, weak FSK signal, as obtained from (a) conventional AN/TRD-4 and (b) modified AN/TRD-4. 


b 
FicurE 11. Pulsed signal, as cbtained from (a) conventional AN/TRD-4 and (b) modified AN/TRD-4. 
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Antenna arrays which are designed to utilize correlation techniques can result in direc- 
tivity patterns with very narrow beamwidth. However, analysis of resolution capabilities 
of these arrays indicates a marked change in expected performance in the presence of two 
or more signal sources. These effects are analyzed for single frequency signal sources and 
for randomly varying signal sources. It is shown that optimum results occur when the 
nonlinear processing of the antenna voltages is limited to a single multiplication. Under 
these conditions the correlation array has a directivity equivalent to that of a linear array 
of twice the length. 


1. Introduction 


Increasing attention has been given in recent years to the merits of correlation techniques 
as a means of improving the resolution capabilities of passive antenna systems. The possi- 
bility of making a decided saving in the overall size of directional antennas, even though this 
saving must be purchased at the cost of increased complexity of the antenna circuitry, has a 
distinct appeal in certain antenna applications. It is the purpose of this paper to outline the 
mathematical analysis of the resolution characteristics of such an antenna array and to com- 
pare these results with those which could be expected from a familiar linear additive array. 

A correlation array will be defined as one in which voltages induced on the elements are mul- 
tiplied together and the resulting voltage is averaged over some prescribed time interval to give 
a desired output voltage. There are, of course, two fundamental variables in this definition. 
In the first place, when there are a number of elements in the array there are a great many 
possible combinations of the element voltages. And, secondly, when there are undesired 
fluctuations in the output voltage, the time interval which is available for averaging this 
voltage can have a primary influence on the resolution capability of the array. 

For this discussion of the mathematical analysis of resolution characteristics the first varia- 
ble will be avoided by limiting the calculations to specific examples of correlation arrays. 
The second variable will enter the analysis and will appear in the final results. 

It will be assumed that the problem at hand is the resolution of two signal sources sepa- 
rated by some angular displacement 6. Both single-frequency and band-limited signal sources 
will be considered. The signal sources are assumed to have identical power spectral density 
over the frequency band of interest. They cannot be resolved by frequency selectivity in the 
receiving system, but must be resolved by appropriate operations on the voltages induced in 
the elements of the array. 


2. Single-Frequency Sources 


In order to analyze the correlation array without becoming too deeply involved in the re- 
quired calculations, a relatively simple array of four isotropic elements will be first considered. 
The pattern of the response of this array as a function of @ can be shown to be equivalent to 
that of a much larger linear additive array. This particular aspect of the capabilities of 
correlation arrays has been reported in the literature [Berman and Clay, 1957; Jacobson, 1958; 
Drane, 1959]. These articles have demonstrated that, in response to a single signal source, 
a correlation array with properly selected element spacing has a pattern which is exactly 
equivalent to that of a much larger linear array. When two signal sources are present, how- 
ever, an interaction of voltages produced by the sources occurs in the multiplication processes 
of the correlator, and this equivalence between linear and correlation arrays has to be reanalyzed. 


1 Contribution from Department of Electrical Engineering, University of California, Berkeley 4, Calif. 
2 This research was supported in part by the Office of Naval Research under Cont ract Nonr-222(74). ase. . 
3 Paper presented at the Conference on Transmission Problems Related to High-Frequency Direction Finding, at UCLA, June 21-24, 1960. 
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The signal source (see fig. 1) located in the plane of the array and at an angular displacement 
6 from the main lobe axis, induces a voltage V; on the 7-th antenna element, where 


V,=A cos (wt+7;) 


(d; is the distance from a reference point in the 


d; : 
7,=— sin 0 . ; . 
gies” line of the array to the ith element.) 








output 


FicurE 1. A four-element correlation array. 


If the bandwidth of the integrating circuit is established to reject the second- and fourth-order 
frequencies resulting from the multiplications, the output voltage is 


s 1 : 
Vour=—~ A‘ [cos w(71—72—73+74) + C08 w(71—T2+73—74) +008 w(71-+72—73—T;) J. 


And if the relative spacing between elements is established with di,=D, d.3;=D, and d34,=2D 


; ae | Sa 
(and letting seals’ ® sin 6): 
- 1 4 a « Ps ry 
V out=5 A! [cos Y+ cos 3X+ cos 5X]. 

But this, except for the constant term, is the voltage pattern for a six-element linear additive 
array with a uniform element spacing of 2D. 

So in this sense the pattern of the four-element correlation array of overall length 4D and the 
voltage pattern of the six-element linear array of overall length 10D are equivalent. 

Berman and Clay have discussed this equivalence between correlation and linear arrays. 
In general, directivity patterns for a single monochromatic source have the following mathe- 


matical equivalence: 


Correlation array Linear array 
4 elements, length 4D 6 elements, length 10D 
6 elements, length 8D 15 elements, length 28D 
8 elements, length 16D 32 elements, length 62D 


However, in resolving two similar sources these patterns cannot be used directly, but must be 
treated with some care. For example, assume source A(V,=E cos wat) is located on the main 
lobe axis and source B with the same amplitude (Vz=E cos wet), is displaced through 6. The 
output voltage of the six-element linear array would be: 


Vour=6E cos wt+2E (cos X+cos 3.X¥+ cos 5.X) cos wat. 
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Whether or not the sources are coherent, the voltage output for arbitrary movement of the 
sources is simply the sum of the voltage patterns of each source considered individually. 

The correlation array, on the other hand, will contain cross product terms resulting from the 
multiplication processes whose form will depend on the coherence of the two sources and, in 
general, on the number of multiplication processes. 

For the example of the four-element correlation array considered here, the output voltages 
for two conditions of coherence are: 


1. For two coherent sources: 


“Ww 


« 
> e 


J al 


1 3 a - 
Ets E*(cos X+cos 3X+ cos 5X) +: E*(3+11 cos X 


Col 


+10 cos 2X+7 cos 3X+6 cos 4X +3 cos 5X+2cos 6X). 


2. For two sources of slightly different frequencies: 
3 4 1 4 r rs r ~ 
Vout=5 E t+ E*(cos X+cos 3X + cos 5X) 
1 ; ; . _ eae ‘i 
+; E‘*(4 cos X+4 cos 2X+2 cos 3X+2 cos 4X) +5 E*(f(cos nX-cos 6t) 


where 6t is the frequency difference between source ''A and source ‘B, and f(cos nX-cos ét) 
represents a number of terms all containing this beat frequency component. 

In both of these equations the first two terms correspond to the linear array pattern, while 
the remaining terms arise from the nonlinearity of the correlation array. In the case of two 
sources with identical frequencies, the cross-product terms are constant with time and time 
averaging of the output voltage will not alter the result. When the sources have different 
frequencies, however, time averaging can be employed to reduce the beat frequency (cos ét) 
part of the cross-product terms; but even in this case the resulting pattern will not be precisely 
equivalent to that of the six-element linear array since some of the cross-product terms will 
remain constant with time. 

Calculation of the resolution characteristics of the correlation array is then more involved 
than the similar calculation for the linear additive array. 


: : ‘ r ; : : 
Assuming a basic spacing D=>; the four-element correlation array of length 4D) will re- 


solve two sources of slightly different frequencies at approximately 19.5°. This is equiva- 

Jent to the resolution of a uniform linear array with an aperture of 9D, about twice as long. 
A slightly more complex example of a correlation array is one discussed by Drane: 

The element on the left in figure 2 represents a uniform linear array with aperture a,, while 

that on the right is a simple interferometer with aperture a). 





Outpiit 


Figure 2. A correlation array employing a linear array 
and an interferometer. 
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sin X 


— where 





The directivity pattern of the uniform linear array is proportional to Y 


py MO Gh . . ’ , 
x=“($) sin @. If a,=a, and there is no spacing between the right end of the linear array 
e\2 
and the left element of the interferometer, the directivity pattern of the correlation array is 


sin 4.X 


4X 


J avi x 





which is the same as that of a uniform linear array of twice the length of the actual array. 
However, if two coherent sources are present the output is 
sin X, cos Xz , cos X, sin ad 


43 sin 4X, , sin 4X, ‘ . | 
y = ) a 
J out & 4X, + 4X, + cos (X A + X B 4X, + 4X, 








The third term again occurs because of the nonlinearity of the correlator. 

Calculation of the resolution capability of this array shows that it is equivalent to a uni- 
form linear array one and a half times as long. 

These have been only two examples of the effect of cross-product terms in a correlation 
array. In each case considered, the effect of these terms would vary, depending on the types 
of signal voltages emitted by the sources and by the number of successive multiplicative proc- 
esses in the correlator. If the sources were emitting complex signals which could be decomposed 
into a number of fixed frequency components (a pulse system with a fixed repetition rate, for 
example) then the cross-product terms resulting from the multiplication would become quite 
involved. 

It should be noted that if the correlator contains only one stage of multiplication, then the 
cross-product terms occurring from two sources of different frequencies appear only as beat 
These terms occur as a low-frequency a-c signal appearing with the 
desired d-c measuring voltage. In this case the cross-product terms can be minimized by time 
averaging the output voltage. However, if more than one stage of multiplication occurs 
between the antenna element and the output of the array (as was the case in the example of 
the four element correlation array) the cross-product terms will occur as low-frequency a-c 
terms and also as d-c terms. In this case the equivalence between the correlation arrays and 
the linear arrays is not immediately apparent but must be determined by calculation of reso- 
lution characteristics for the particular type of signal encountered. 


frequency components. 


3. Randomly Varying Sources 


In those situations in which there is at best only incomplete information about the manner 
in which the voltage produced by a signal source will vary with time, it is necessary to consider a 
suitable statistical model which will impose bounds on the expected results and will provide an 
average description of the voltage variation. The most general such model suitable to this an- 
tenna problem is that of a band-limited voltage with a normal distribution. The source can 
then be described as one which produces a randomly varying voltage whose Fourier series rep- 


resentation becomes: 


N 
V(t)=>Dd 5 ¢, cos (w,t+,). 


n=1 
V(® is distributed normally with mean zero. 


c, has a Rayleigh distribution with e=2W(f,)af, where W(f,) is the power density over a 
frequency interval 6f centered at /,. 


®, has a uniform distribution over (0, 27). 
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wo is the lower edge of the band of width 24f ¢/s and w= wot Fr 


The directivity pattern of the correlation array can be stated in terms of the correlation 
coefficient of the voltage produced by the signal source. For a four-element correlation array 
the expected value of the output voltage is: 


Vi(t) Vo(t) Va(t) ValE) 
and, since these are normally distributed voltages, this becomes: 
MaOVOIWOVOlFAMOVOMW OVO] +MOVAOIV:OV30)] 
=R (71) R (r34) +R (113) R (724) + R (14) R (r23) 








where R(r;;) is the correlation coefficient of the voltages induced on elements 7 and /. 
If the signal source is band limited, f,+ Af, and the signal has a uniform power spectral 
density over this pass band, then 


sin 2xAfr 


R()=RO) Sr, 








cos 2nf,r. (autocovariance) 


With proper spacing of the correlation array elements this again can be put in a form resembling 
the directivity pattern of a uniform linear array. 





Vi (t) Vo(t) V(t) Va(t) =3R (O)? [A cos X+B cos 3X+C cos 5X] 


ee : ip sin 2rA fr 
where X=2rzf, > sin 6 and the coefficients A, B, and C are combinations of —5 — terms 
oT. T 





which approach the value 1 as the receiver pass band is decreased. 

If there are two such sources separated by an angular displacement 6, the output voltage 
of a linear array could be suitably employed to resolve the sources. An instantaneous voltage 
observation at various positions as the array is rotated past the sources would generally not 
be an adequate procedure due to the assumed randomness of the sources. Some sort of time 
averaging procedure (as in an average power measurement) or a procedure requiring a definite 
observation interval (the peak voltage occurring in the interval, for example) would be 
necessary. 

As an example of this use of a linear array, assume that the output of the array will be 
squared and resolution obtained by power measurement. As before, assume there are two 
sources, source A on the main lobe axis and source B displaced through 6°. Then before 
squaring: 


Vour=6V4+2V2(cos X-+ cos 3X+cos 5X) 


N N 
Vout=6 >> Cn COS (w,t +¢,) +2 >) d, cos (nt+¢,) (cos X+cos 3X+cos 5X) 
n=1 n=1 
After squaring: 
V2..=36V4+24V4Ve (cos X+cos 3X+cos 5X)+4V}2(cos X+cos 3X+cos 5X)? 


N N 
V24n=36 | >) D>) erm COS (w,t+¢,) cos (ont +n) | 


n=1 m=1 


+24 = DD) Cndm COS (Wnt +,) COS (Wmt+bm) (cos X+cos 3X-+cos 5X) | 


+4] 53 DD didm cos (wpt+¢,) COS (Wpt+m) (cos X+cos 3X+cos 5Xx)'| 


n m 
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The first and last terms of this expression have xX? distributions with mean values equal to the 
power from the signals; the middle term has a bivariate normal distribution with mean zero. 
So the average value of the squared voltage term will be equal to the power received from the 
two sources. The instantaneous value will fluctuate about the average, and time averaging 
can be used to minimize the effect of the fluctuation. 

An expression for the effect of an averaging device can be obtained from general filter 
considerations [Middleton, 1960]. If h(¢) is the effective weighing function of a linear measuring 
device and z(t) is the function to be measured, a measurement M,(7) made at time ‘=T after 
z(t) has been introduced at t=0 can be expressed by the convolution 


r 
M,(1)={ h(u)x(T—u)du h(u)=0, u<c0 
0 
where (0,7’) is the observation interval. 14,(7) will vary from observation to observation, 


fluctuating about the expected value M,(7) with a variance o?=M,(T)?—[M,(T)]’. 
In the general situation, z(¢) is at least wide sense stationary, and 





—_—_———_—— T a 
MT=| h(u)a(T—u)du 
0 





4 T 
MiTy= | i} h(u)a(T—u)2x(T—v)h(v)du dv 
0 0 


T of 
=| q h(u)¥,(u—v)h(v)du dv 
0 0 


where y, is the autocovariance coefficient of the function to be measured. For an ideal inte- 


grator h(u)=7 O0<u<T. And 


M,(1)=z(t) 


MAT I=> | * (1-5) eds 


The correlation coefficient of the power function can be obtained in terms of the correlation 
coefficient of the voltage. The voltage has a normal distribution; the correlation coefficient 
of the power. 

This is, of course, a direct indication of the fluctuation of the output power about the mean 
value. For this case of two statistically similar sources, the variance determines, in terms of 
confidence intervals, the difference in mean power which the array receives when it is directed 
at one of the sources and when it is directed between the sources in order for the sources to be 
consistently resolved. 

Now, with this description of the resolution process by a linear array it is possible to 
investigate the resolution process in a correlation array and to compare the resultant effects. 
As in the previous case, it is necessary to determine the mean value of the output voltage and 
to determine the correlation coefficient of this voltage in order to establish bounds on the 
expected fluctuation of the measurement. 

The expected value of the voltage produced at the output of a four-element correlation 
array by two independent sources (source A and source B) can be calculated either directly 
or by the characteristic function method and can be expressed in terms of the correlation 
coefficients of the individual voltages. 

With two sources present, the voltage on the ith element becomes: 


Vi(t)=Va(t+7;)+ Valt+71). 
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M.\ T) = if ¢ (t) \ >(t) | 3(t) if “(t) — 1 [Ra (ri) R4 (734) “« i R.4(713) ee (T24) —" ies (714) R, (723)] 
. [Pg (t12) Ra (734) Fa R3(713) Rp (724) a Ral 714) Rp (723)] ei [Rs (734) Rp (T12) “- Rs ( 712) Rp ( 734) 
+ 5 Rs (rT24) Rp (713) = R, (T23) Rp (714) aS Rs (714) Re (723) i Rs (713) Ral 2 )] pe 





a (iv)? . ; Stas ; 
V(t)?V (t+s)?, will be the coefficient of 5 In the expansion of the characteristic function 


of V(t),V(t+s). 





C(u,v,s)=exp [iwV (t)+2eV (t+s)]=exp [—4R(0)(u2?+ 0?) — R(s)ur] 
where /?(s) is the correlation coefficient of the voltage. Then 
¥(s)= R(0)?+2R(s)?. 


If it is assumed that the receiver circuits have a rectangular pass band, f.+A/, and that the 
power spectral density of the signal sources is uniform over this pass band, then 


R(0)=average power from the signals 


sin 2rAfs, oo -. 
R(s)=R (0) - 5) nafs cos 2rf,s. 


Substituting these values into the equations describing the ideal integrator: 
M,(T)=R(0) 


9 


r se) ; 
Magy=< (1-5) (R (0)?-+2R (s)2)ds 





am QeAfey gn 4 vay 
fate (- ee mae ) (cos nfs) as | 


which, for large values of 27A fT, becomes 





MTP =D (1+ 5-479) 
The variance of this averaged power then is 


ota [1,(T) = 


fT R(0)?. 


2 aT 
The first and second terms of this expression give the expected voltage due to each source 
individually. The third term contains the cross product components which occur as a result 
of the two stages of multiplication that the element voltages undergo. 

The variance of this output voltage could be calculated directly; however, the large 
number of terms in the final expression would make this quite laborious. A good approxima- 
tion to this variance can be made by calculating the variance of the voltage produced by a 
single source located on the main lobe axis. The correlation function of this voltage can be 
determined readily by the characteristic function method: 


v(s)=9 R(O)'+72 R(O)?R(s)?-+24 R(s)*. 
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If we again assume an ideal integrator and a band limited process, the variance becomes: 


» (14 op 
73 xAsT) [M(T)}*.  (2"afT)>>1. 


This variance is seen to be somewhat more than four times that found for the single multipli- 


cation required in the power calculation. And, in general, if more elements are added to the 


correlation array, this variance will increase approximately by a factor of four for each pair 
of elements added to the array. 

This increase in variance indicates an increase in the fluctuation component of the out- 
put voltage as the number of multiplicative processes are increased. This results in an in- 
crease in the integration time required to reduce this fluctuation to some prescribed level, 
and, therefore, makes this time averaging interval an important part of any discussion of the 
resolution capabilities of these arrays. 

Band-limited random signals 


Resolution at the 95% confidence level 


Four-element correlation array 
length 4D 
en Aperture of equivalent 

uniform linear array 


.6D 





Resolution Integration time 


19° 2rAfT= 56 9 

18° 2rAfT= 303 9.4D 

17. 5° 2rAfT=2790 9. 2D 
ny 


4. Summary 


An analysis of the expected performance of antenna arrays which utilize correlation 
techniques indicates the possibility of a marked saving in antenna size. In the general situa- 
tion, the multiplicative processes in the correlator will introduce cross-product terms which 
will appear as low-frequency fluctuations in the antenna output voltage. Additionally, these 
cross-product terms may contribute to the d-c output voltage. The appearance of the cross- 
product terms complicates the calculation of resolution capabilities, and generally necessitates 
the definition of a time-averaging interval in any discussion of this resolution. In the ex- 
amples considered here, it was evident that the directivity pattern alone did not describe the 
performance of a correlation array in resolving two signals. In these examples the relatively 
simple correlation arrays had resolution capabilities equivalent to those of uniform linear 


arrays of about twice the length. 
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Instrumentation for Propagation and 
Direction-Finding Measurements ° ” 


Edgar C. 


Hayden 


(August 26, 1960; revised January 3, 1961) 


Limitations imposed on radio direction-finding systems are discussed in terms of a 


generalized representation of such systems in the form of a block diagram. 


Factors affecting 


these limitations include: (1) considerations of signal-to-noise ratio in the early part of the 
system, (2) receiver bandwidth requirements for adequate selectivity, (3) width of spectrum 
generated by prereceiver encoding and computing processes, (4) restriction to linear processes 
in multisignal portions of the system, and (5) availability of operational devices suitable for 


use in low-signal-level portions of the system. 


A generalized representation of a radio direction- 
finding system can be given in the form of a block 
diagram or flow chart (fig. 1) showing the operations 
or functions which are necessary in translating the 
basic measurements on the incident signal field into 
the values of the desired incident field parameters. 
Limitations are imposed on the system at each stage 
in the process. Factors leading to these limitations 
include: (1) Considerations of signal-to-noise ratio 
in the early part of the system, (2) receiver band- 
width requirements for adequate selectivity, (3) 
width of spectrum generated by prereceiver encoding 
and computing processes, (4) restriction to linear 
processes in multisignal portions of the system, 
and (5) availability of operational devices suitable 
for use in low-signal-level portions of the system. 

The development of the detailed configuration 
of any direction-finding system is based on some 
concept of the character of the incident signal field. 
This character is described in terms of a mathematical 
model, the parameters of which are the quantities 
to be evaluated by the system. To this point, the 
incident signal field models which have been used 
as bases for direction-finding system development 
have been rather simple, far too simple in fact to 
satisfactorily describe the actual signal fields. One 
of the frontiers in direction-finding system develop- 
ment lies, at the moment, in the area of improving 
the incident field model and of developing more 
sophisticated systems based on the improved signal 
field models. 

Another frontier lies in the area of applying modern 
computing techniques to direction-finding system 
design. One of the major functions of a direction- 


POSSIBLE MULTIPLE CHANNELS 





re — = ———— \ ———E _ 
COMPONENT | ANTENNA | _| processor |‘ | recewer | | processor L { output 
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Figure 1. Generalized direction-finding system. 





1 Contribution from Radio Direction Finding Research Group, Department 
of Electrical Engineering, University of Illinois, Urbana, Ill. 





? Summary of paper presented at the Conference on Transmission Problems 
Related to High-Frequency Direction Finding, at UCLA, June 21-24, 1960. 


finding system is, after all, the performance of a 
computational operation. To this point all systems 
in use are based on some simple-minded analog 
computation scheme. This is an _ unnecessarily 
restrictive limitation, especially considering the 
tremendous strides which have been made in recent 
years in the field of electronic computation. 

Two specific direction-finding systems are pro- 
posed, both intended primarily for use with antenna 
arrays of wide aperture. In one system, intended 
principally for circular arrays, a commutating device 
is used to encode the output of the antenna array 
in a form suitable for transmission through a single 
receiver channel. The output of the receiver would 
be analyzed, using sampled data techniques, to 
obtain the required bearing information. The 
second system makes use of a twin-channel receiving 
system to avoid the necessity for relatively slow 
scanning, thus permitting essentially instantaneous 
acquisition of bearing data (by “instantaneous’’ is 
meant a time small enough to permit obtaining a 
bearing on the shortest pulse transmittable through 
the receiver). The unusual feature of both systems 
is the use of an electronic digital computer to perform 
the bearing interpolation between planes of symmetry 
of the array, and to provide the final numerical 
output of the system. The conventional “indicator” 
is retained in modified form as an operator aid, and 
to permit the operator to “censor’’ information fed 
to the computer, or to “backstop” the computation 
system in situations where it would fail. 


Nore: A more detailed discussion of the material 
presented in this paper will be available in printed 
form during the first quarter of 1961. It will 
appear in Interim Engineering Report No. 9, Wullen- 
weber Direction-Finding System to be published by 
the Radio Direction-Finding Research Laboratory, 
Department of Electrical Engineering, University 
of Illinois, Urbana, Illinois, for the Bureau of Ships, 
under Contract Nobsr 64723. 

(Paper 65D3-127) 
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Brooke Variance Classification System for DF Bearings ** 
E. M. L. Beale 


(August 22, 1960; revised December 27, 1960) 


This paper describes the advantages of having an objective classification system for 


DF bearings. 


The Brooke system is described in some detail, and the problems involved 


in setting up a system on these lines are considered. 


1. Introduction 


For many years we used a simple ABCD classi- 
fication system for HF DF bearings taken with 
U-Adcock direction finders. This was not entirely 
satisfactory for various reasons. The fundamental 
difficulty was that, although in theory the classifi- 
cation depended in a well-defined way on the con- 
sistency of the observed bearings, in practice it was 
largely subjective. The Brooke Variance Classifi- 
cation System is based on the Ross-Barfield system 
developed during the Second World War and 
described by Ross [1947]. The present paper is an 
attempt to describe the system from the point of 
view of a group that might be considering the intro- 
duction of such a system. 

Section 2 of this paper reviews the reasons for 
having a bearing classification system. The sources 
of error in DF bearings are described briefly in 
section 3. The history of the Brooke system is 
reviewed in section 4. Section 5 is devoted to a 
discussion of the problems involved in introducing a 
similar system into another network. The solutions 
found for the Brooke system are indicated. Section 
6 describes the statistical analysis required to esti- 
mate variance components from check bearing data. 


2. Why Does One Need a Classification 
System? 


The first question is, why does one need a bearing 
classification system? It is impossible to give a 
short answer to this question, since it depends on 
the answer to an even more fundamental question: 
What use is one trying to make of one’s DF bearings? 
One will presumably plot the bearings reported on 
any task on a gnomonic chart. One can then either 

(1) gaze in wonder at the set of plotted lines, and 
finally emerge with a point estimate for the position 
of the transmitter, and a circle of arbitrary radius 
centered on this point such that the true position 
“probably” lies within the circle, | 
or (2) use one of the many unsound plotting methods 
that have been proposed from time to time, 


1 Contribution from Admiralty Research Laboratory, Teddington, Middlesex, 


England. be 
Paper presented at the Conference on Transmission Problems Related to 


High-Frequency Direction Finding, at UCLA, June 21-24, 1960. 


| or (3) use an electronic computer to find a more 
precise least squares solution, i.e., to find 
(a) a Best Point Estimate (or Most Probable 

Point) defined as the point minimizing the weighted 

sum of squares of the angular errors, the weights 

being inversely proportional to the assumed 
variances of bearing errors at the station concerned, 
and 

(b) a region, normally a circle or a rectangle, 
approximating to the region where this weighted 
sum of squares exceeds its minimum value by less 
than some given constant, such as 4. The region 
may be called a 90 percent (say) Probability 

Region, or Confidence Region. Its purpose is to 

give an appreciation ef the probable accuracy of 

the point estimate,’ 
or (4) use a method based on some valid approxima- 
tion to the least squares solution. 

In cases (1) and (2) it may be helpful to have some 
rough measure of the accuracy of any particular 
bearing, but the traditional ABCD classification 
based on the consistency of repeated observations, 
or some other subjective method, may be good 
enough. 

In cases (3) and (4) it is important to have a 
realistic measure of the probable accuracy of any 
particular bearing. Ross [1947] wrote 


“It has long been customary for bearings to be classified 
somewhat arbitrarily by the direction-finding operator. This 
classification, intended to give an indication of the probable 
reliability of the bearings, is usually based on the quality 
(sharpness) and steadiness of the minimum. This practice 
is a legacy from the days when radio communication was 
conducted chiefly on the long and medium wavelengths, and 
it was not unreasonable to expect the operator to estimate the 
reliability of bearings in this way. In the short-wave band, 
however (20 to 100 m approximately), where propagation is 
chiefly controlled by ionospheric reflections, conditions are 
far from simple. Bearings are generally in a state of con- 
tinuous variation and it becomes increasingly difficult for an 
operator to classify them by any more or less intuitive 
process.” 


3 The conventional justification for this as a Confidence Region is based on 
the assumption that all possible bearing lines from a given station are approxi- 
mately parallel. Beale [1960] provides a method of justifying it even when this 
approximation is unacceptable. Beale’s paper shows that the approximate 
confidence region is justified if a quantity Ng, called the “intrinsic nonlinearity 
of the model” is smaller than about 0.1. It seems likely that the DF problem 
nearly always satisfies this condition; even though one can easily imagine situ- 
ations with a “‘narrow base line’”’ when the problem is decidedly nonlinear when 
expressed in terms of natural parameters, such as the latitude and longitude o 
the transmitter, 
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Ross proceeds to describe an objective classifica- 
tion system which is the basis of the Brooke system. 
The details of such systems will be described in 
section 4. Meanwhile we should consider another 
advantage in having a rational and objective 
classification system: it is almost indispensable to a 
thorough quantitative analysis of the magnitudes of 
the various errors arising in practical DF. It is 
virtually impossible to get a set of observations that 
is statistically balanced with respect to all relevant 
factors. So, for example, when studying the effect 
of distance one must allow for the difference in fre- 
quencies of transmitters at different distances; and 
such a process requires the estimation of variance 
components due to the different errors making up 
the observed bearing error. 


3. Sources of Error in HF DF Bearings 


The sources of error in HF DF bearings are dis- 
cussed by Ross [1947] and Bowen [1955]. They can 
be summarized as follows: 

(1) Instrumental errors due to small errors of 
antenna balance, and similar causes. These errors 
can be greatly reduced by a local calibration, but 
there may be an appreciable residual error varying 
with radio frequency and direction of incidence (in 
azimuth and elevation) and from site to site. 

(2) Distant-site errors, due to irregularities of, and 
obstacles on, the terrain around a high-frequency 
direction finder but outside the near site or calibra- 
tion area. The values of these errors vary with 
radio frequency and direction of incidence (in 
azimuth and elevation). The variances vary from 
site to site, and tend to decrease as the radio fre- 
quency increases. 

From our present point of view, errors of types 
(1) and (2) are similar. They are virtually inde- 
pendent of time. For a given transmitter working 
on a given frequency the errors may vary a little 
in that the angle of elevation of the incoming ray 
will depend to some extent on ionospheric conditions, 
but there will be an appreciable systematic compo- 
nent of such an error. 

(3) Ionospheric lateral deviation errors, due to 
large scale irregularities in the ionosphere. These 
errors vary only slowly with time (with periods of 
the order of tens of minutes). Their variance does 
not depend on the DF site, but depends strongly on 
the distance of the transmitter. Not much is known 
about the effects of time of day and radiofrequency 
on lateral deviation. It seems possible that errors 
are somewhat greater at night than by day and some- 
what greater on low frequencies than on high fre- 
quencies. If so, this is probably due to the increased 
number of multiple reflections possible at night and 
on lower frequencies. Large systematic and random 
effects can be observed around sunrise. Later 
studies have not essentially changed the views ex- 
pressed by Ross [1949] on the subject of lateral 
deviation. 

(4) Wave interference errors and polarization 
errors. These errors may be quite large, but they 








generally vary rapidly (with periods of the order of 
seconds). Their effect can therefore be greatly 
reduced by taking the average bearing (preferably 
weighted according to signal strength) over a period 
of 10 sec or longer. 


4. History of the Brooke System 


The Brooke Variance Classification System arose 
out of the system developed during the Second 
World War by Mr. W. Ross and Dr. R. H. Barfield, 
and described by Ross [1947]. The basic assumption 
underlying the system is that the error in any ob- 
served DF bearing is a sum of three uncorrelated 
errors: 

(1) An instrumental and distant site error, whose 
mean square value can be determined as a function 
of the radiofrequency and the DF site. 

(2) A lateral deviation error, whose mean square 
value can be determined as a function of the distance 
from the transmitter to the DF site. This component 
must of course be estimated at the plotting center. 

(3) A wave interference and polarization error. 
An estimate of the mean square value of this error 
can be determined from the mean ‘‘swing”’ of the 
observed bearing over a period of 1 or 2 sec, or from 
the scatter of repeated observations made at intervals 
of about 5 or more sec, or both. The distance from 
the transmitter to the DF site can also be used to 
throw some light on the probable magnitude of this 
error, since wave interference due to multipath 
arrival is apt to be more serious on the more distant 
targets. 

When after the war attention was again directed 
to the classification of HF DF bearings, the Ross 
system was examined but considered to be unsatis- 
factory in its existing form. The users also con- 
sidered the subjective method of classifying bearings, 
then in use, to be unsatisfactory. Mr. Norman 
Brooke, while a member of the Admiralty Depart- 
ment of Operational Research, carried out an investi- 
gation into the problem of classifying HF DF 
bearings, the outcome of which was a proposal to 
modify the Ross system, and to extend it to apply 
to visual, as well as aural null, equipment. He 
found that it was necessary to give the variance 
components values that were rather larger than 
those proposed by Ross. 

The modified version was found to be satisfactory. 
It was decided, with Mr. Ross’s concurrence, that 
the system should be called the Brooke system. 

Before the system was introduced a statistical 
analysis was made of a large number of check bear- 
ings on known transmitters, in order to derive 
realistic estimates of all the variance components. 
Further, it is known that station performance 
changes to some extent with time, so the analysis 
of check bearing data continues on a routine basis 
in order to keep track of these changes. It is found 
that the variance component depending on the DF 
station has to be changed by 1 unit (i.e., 1° squared) 
from time to time—perhaps once every 2 years on 
the average. The other variance components are 
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reviewed from time to time, but changes have not 
been found necessary. 

So we may conclude that the Brooke system works 
satisfactorily. The essential reason for this seems 
to be that the system combines physical theory and 
empiricism in about the right proportions for our 
present state of knowledge about HF DF. We 
certainly do not maintain that no further improve- 
ments in the system are possible. 


5. How Could a Similar Classification Sys- 
tem Be Introduced Into Another Organ- 
ization? 


Suppose that some authority has decided that it 
would be worth introducing an objective bearing 
classification system on these lines into his HF DF 
network, what steps must be taken to implement 
this decision? 

The first task is to define a set of factors on which 
the variance is assumed to depend. It would seem 
natural to include 

(a) The DF station, 

(b) the radiofrequency, 

(c) the distance of the transmitter from the 

station, and 

(d) some measure of the quality of the observed 

bearing. 

Some discussion on this last item is perhaps 
required. A cursory reading of Ross [1947] might 
give the impression that an objective classification 
system is only possible if one takes repeated snap 
observations to determine the bearing. This is not 
so. The true situation is as follows: 

Admittedly if one takes repeated snap observations 
(which may be mental averages over a period of 1 
or 2 sec), then one is likely to obtain a more accurate 
bearing than if one simply gazes at the display for 
about 30 sec, because experience has shown that it 
is almost impossible to average mentally over a 
period of more than about 5 sec. Furthermore, the 
scatter of the individual snap bearings provides 
additional information concerning the reliability of 
the mean bearing. But if the transmitter is not 
active for long enough to take repeated observations, 
or if the bearing is required so urgently that there 
is not time to take and report the results of several 
snaps, then one can still use an objective classification 
system based on other factors. Indeed both the 
original Ross-Barfield system and the Brooke system 
include means of assessing a variance when it is 
impossible to obtain more than one snap. 

The next task is to postulate a formula for the 
variance as a function of these factors with a finite 
number of unknown parameters. 

It seems natural to write 


B=V,+ Vit V,, 
where B (for Brooke) denotes the total variance (i.e., 


the total mean square error if one proceeds as if all 
errors had zero means). 





V, denotes a component depending on radiofre- 
quency and the DF station, representing instrumen- 
tal and distant site errors, 

V2 denotes a component depending on the distance 
from the transmitter, representing lateral deviation 
error, and V, denotes a component depending on 
the quality of the observed bearing, and the type of 
equipment and display used, representing wave in- 
terference and polarization error, and also observa- 
tional error due to imperfections in the DF operator. 

We must still specify the forms of V;, Va, and V, 
to reduce them to a finite number of parameters. 
It may be helpful to consider each in turn with 
reference to the existing Brooke system. 

The Brooke system assumes that 


V,=a,+byf, 


where a; and b; are parameters estimated separately 
for each DF station, and f denotes the frequency 
band. The definition of these frequency bands is 
being reconsidered. There is some evidence that 
the variance is not very sensitive to radiofrequency 
above about 9 Me/s. The following bands may 
therefore prove appropriate. 


f=1 for frequencies between 2.000 and 3.999 Me/s, 
f=2 for frequencies between 4.000 and 8.999 Me/s, 
f=3 for frequencies above 9.000 Me/s. 


For the distance component, the Brooke system 
assumes that 


7,=0.6/d? for d<0.8, 
7.= for 0.8<d<3.5, 
=! for 3.5<d, 


where d denotes distance from the transmitter to 
the DF station, measured in thousands of kilometers. 

For the quality component, the Brooke system 
assumes that 


Vi.=V,+V,., 
where 


V,=0 for swings between 0 and 8°, 
=1 for swings between 9 and 13°, 
=2 for swings between 14 and 18°, 
=4 for swings between 19 and 23°, 
=6 for swings between 24 and 37°, 
=9 for swings over 38° or unmeasured, 
and V,=72/nd,?, 


where n is the number of snap bearings taken, r is 
the range of these snaps, and d, is the average range 
of n observations from a Gaussian population with 
unit variance (tabulated by Tippett [1925)). 

The theory underlying the use of these formulas is 
that V, might be expected to give a valid estimate 
of the variance of those components of error (due to 
polarization errors and wave interference) that vary 
rapidly. But in practice the operator cannot help 


257 








being biased by the knowledge that all his snaps 
refer to the same transmission. He will therefore 
tend to make his observations more consistent than 
they should be, and an additional variance compo- 
nent must be introduced to allow for this. The 
bias is likely to depend on the difficulty of the task 
of taking a thes which is measured by the 
average “swing’’. For an aural null equipment the 
swing is defined as the width of the arc about the 
minimum over which there is a just detectable change 
of signal. For any equipment with a visual display 
it is defined as the angle through which the instan- 
taneous value of the bearing moves during the time 
taken to read its mean value for an individual snap. 
There is of course no logical reason why the formula 
for V, should be the same for different types of 
display, even though all the equipments have the 
same antenna system. But the statistical analysis of 
check bearings suggests that there is no important 
difference in this respect between aural null, spinning 
goniometer, and twin-channel c.r.d.f. displays.* 

If only one snap is taken, increased values of Vz 
and V, are used to compensate for the absence of 
the V,, term. 

Different organizations may find it convenient to 
estimate V, differently. This rust depend on what 
information the DF operator finds it convenient to 
report. If he simply reports a letter classification 
(A, B, C, or D), then one could put V,=0 for an A 
class bearing, and estimate 3 constants, being the 
values of V, for B, C, and D class bearings 
respectively. 

One might also use the signal strength as an 
indicator of bearing quality. This does not in fact 
enter into the existing Brooke system. 

Having defined the form of the bearing classifica- 
tion system, one must estimate the parameters. 
One can get preliminary estimates of these from 
physical theory and from the values used in existing 
systems, but these estimates must be corrected (or 
confirmed) using data obtained by the DF organiza- 
tion in question. These data may be either genuine 
operational data or check bearing data where the 
true positions of the transmitters are known to the 
analyst. 

As far as I know, no serious attempt has yet been 
made to estimate variances from operational data 
when the true positions of the transmitters are 
unknown. The use of such data poses the funda- 





4 There is a small point of some statistical interest about the estimation of Vq. 
Elementary theory suggests that ¥ V, is an unbiased estimator of the standard 
deviation of errors varving randomly from snap to snap, for a given value of this 
standard deviation. But this does not prove that ¥ V, provides the best estimate 
of this standard deviation, given the range r. Indeed one might expect that a 
better estimate would be obtained by attaching some weight to the estimate 
obtained from the long-term average value of this standard deviation—thereby 
making the used estimate less dependent onr. After holding forth on these lines 
to my colleagues on more than one occasion, I eventually persuaded them to doa 
detailed analysis of the check bearing data to investigate this point. They 
found that, if anything, the variance component V, should be made more de- 
pendent on r. It has been suggested that this is because the dividing factor 
d,2 is too large because the error distribution is not normal. (And indeed even 
from a normal population a smaller dividing factor is required to give an unbiased 
estimate of the variance as opposed to the standard deviation.) But I prefer to 
regard this as suggesting that the range has additional value as an indicator of the 
quality of the bearing. More light could be thrown on this point by an analysis 
of the individual snap observations, but it is of little operational importance. 














mental difficulty that one has to lean rather heavily 
on the assumption that the bearing errors are uncor- 
related. If one is prepared to accept this assump- 
tion, then it is fairly easy to estimate the absolute 
values of the variances of a set of DF stations if we 
know their relative variances. For in this case we 
can determine the true Best Point Estimate; i.e., the 
hypothetical position of the transmitter minimizing 
the sum of squares of the bearing errors divided by 
their variances. We can then use the fact that (if 
the error distributions are normal with zero mean) 
the minimum value of this sum is distributed as x? 
with n-2 degrees of freedom, where n is the number 
of DF stations contributing to a task. In particular, 
the mean value of this sum is n-2.’ 

In practice the relative variances are unlikely to 
be known. The problem then becomes considerably 
more complicated. It is discussed in a companion 
paper, Beale [1961]. 

It is much easier to estimate variances from check 
bearings, provided satisfactory check bearings can be 
obtained. 

Check bearings may be used both to verify that 
the equipment is functioning properly and to indicate 
the probable accuracy of the DF station on opera- 
tional tasks. A check bearing program to check the 
equipment does not need to be controlled nearly as 
carefully as a check bearing program to estimate 
accuracy. Bowen [1955] writes 


“The use of check bearings on known transmitters in 
establishing the performance of a high-frequency direction- 
finder is a widely accepted technique. It has been found, 
however, that two problems present themselves: 

(a) There is difficulty in arranging a programme in which 
the radio frequencies, ranges, azimuths, propagation paths, 
transmitter powers and other parameters are adequately 
sampled and distributed similarly to those which will define 
the normal task of the direction finder. 

(b) Unless great care is taken to hide the identity of the 
transmission from the operator, both as regards its true 
bearing and its identity with previous transmissions, con- 
siderable operator bias will occur. 

“The methods by which these problems are overcome will 
generally be peculiar to the task with which the d.f. organi- 
zation is concerned. It may be that (a) will require the 
deliberate setting-up of a transmitting system, although care- 
ful selection of known trans nitters may obviate this; (b) can 
be overcome either by careful organization of the way in 
which check-bearing tasks are fed to the operator, or by a 
system of random scale displacements at the direction-finder.” 


Even if one fails to overcome these problems 
completely, it would seem unreasonable to use this 
as an excuse for not introducing an objective classi- 
fication system. Indeed the results from an objec- 
tive system should be very much better then those 
derived from overall station performances based 
on such data—since the system will give some 
rational method of extrapolating to conditions met 
only rarely in the check bearing program. 

The statistical problems involved in estimating 
variance components from check bearings are dis- 
cussed in the next section. 


§ Provided, as will usually be the case, that the quantity N¢ introduced by 
Beale [1960] is small. 
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6. Estimating Variance Components From 
Check Bearing Data 


The fundamental statistical problem in this work 
can be expressed as follows: 

Given n independent observations, each from a 
normal population with zero mean, such that the 
ith observation x; has a variance 


B= Deus 
j=1 


where the a; are all known but some at least of the 
v, are unknown, estimate the unknown »,. 

Before showing how to solve this problem, I 
should perhaps indicate its relevance. The quantity 
az, represents the actual error in the ith bearing 
observation, B; represents its assumed variance, and 
the v; represent the parameters in the classification 
system. For example, suppose that we have a 
simple model with one component that can be 
assumed known, and unknown components depend- 
ing on the radiofrequency and quality, assumed to 
depend on 4 parameters, such that 


(6.1) 


V= V1 +foe, 
where f denotes frequency band, either 1, 2 or 3, and 
V,=0; if the swing <8°, 
V,=y if the swing >8°. 


Now suppose that the first few check bearings 
have the following characteristics: 


Known Vari- Frequency 
Bearing No. ance Component Band Swing 
1 6 1 3°, 
2 3 2 13°, 
3 5 7°, 


etc. 
Then, if we write formally 73;=1, (6.1) is satisfied 
with 
a%,=1, ao=1, a3=1, ay4=0, a5—6, 


Ay =1, A22=2, a23=0, ay=1, a=, 


a3,=1, A=, a33—=1, O34=0, 035=5, 


“9 


the parameters to be estimated being 1, %, v3, and v4. 

We assume that our bearing errors are measured 
ni degrees. Variances are therefore measured in 
degrees squared. 

Now B, is intended to represent the mean value 
of z7, so it seems natural to try to minimize the sum 
of squares of the deviations of the x? from the corre- 
sponding B;. But we must remember that the ex- 
pected magnitudes of these deviations depend on 
the variances of the z?. In fact if x; is normally 
distributed the variance of 2? is 2B}. It therefore 
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seems reasonable to try to minimize 


n 
240i (2i— B,)’, (6.2) 


where the weights w; are regarded as constants but 

where their numerical values are given by © 

w,= 1/B?. (6.3) 

Differentiating (6.2) with respect to v,, regarding 
w; as constant, we have, using (6.1) 


Wie: (x37— B;) =0, (6.4) 
i=1 


for each 7 such that v; is to be estimated. 

In particular if for some j, a;; is always either 
0 or 1, the corresponding equation from (6.4) 
becomes 

>jw.(77—B,)=0, 


where 2; denotes summation over all observations 
such that a j=1. 

The problem is therefore to find estimates 0, for 
the unknown 2»; satisfying (6.1), (6.3) and (6.4). 
This will have to be done iteratively. Choose trial 
estimates for the v;, compute the corresponding 
values of B;, from (6.1), and hence compute the 
weights w; from (6.3), substitute these values of 
w, in (6.4) and hence estimate new values of the 2, .6 
The whole cycle can be repeated if necessary, but 
it should be remembered that small (relative) errors 
in the w; will be of little significance. 

It can easily be shown that these estimates are in 
fact maximum likelihood estimates for the v; The 
above derivation is less rigorous but seems more 
intuitive.’ 

One important proviso should be made. It is 
undesirable to have weights for individual bearings 
of the same order of magnitude as the sum of the 
weights for all other bearings put together. In the 
analyses for the Brooke system, this is avoided by 
putting w;='% whenever B;<V¥2. So the used 
version of (6.3) reads 


w,=min (14, 1/B?). 


It is of interest to consider the probable ac- 
curacy of the final variance estimates. It is easy 
to find the variance of an individual component, 
say v;, if we assume that all the other components are 
known exactly, and also that each w; equals the 
reciprocal of the square of the true variance of the 


6 In practice it may be more convenient to improve the estimates of the param- 
eters referring to one variance component, say Vy, assuming that the others 
have their present trial values. This reduces the number of equations that have 
to be solved simultaneously. : 

7 Some statisticians may point out that for this problem the method of maxi- 
mum likelihood gives an ‘‘inadmissible”’ estimator, since it may produce negative 
variance components. Wiiha fair-sized sample this is unlikely to happen; but, 
if it does, one must replace the negative component by some more plausible value 
and re-estimate the other components accordingly. 
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corresponding 2; For (6.4) then becomes 


n 
Dy Wiltz (2-5 aust; = 0, 
i=1 j=1 
so that the estimate 6, of x is given by 
n 
b= Zy Wean(ti—D ayes) /Bwee'r. 
i= j#k 


Now var 23=28?, and w,=1/B}, so 


var =>) a 2B / > ms 
i=1 i=1 37 
=2/ 5 (s 
aI DI wie. (6.5) 
/ 
If a,=0 or 1 for all 7, (6.5) reduces to 
var 7,=2/2,0;, (6.6) 
where ~;,' denotes summation over those check 


bearings for which a,=1, i.e., for which the given 
variance component is relevant. 

Errors in estimating the weights w;, and also 
errors in the other variance components, will not 
normally have a major effect on the variance of 
dx. On the other hand the fact that the observa- 
tions may not be strictly independent, and that the 
error distribution may not be strictly normal, may 
have a considerable effect on the validity of these 
formulas. Our experience suggests that, presumably 
because of these causes, formulas (6.5) and (6.6) 
underestimate the variances of the #, by a factor of 


about 2. The formulas 
var },~4/2w,ai,, or var },~4/2,wi, 


are therefore recommended in practice. 

There is one pecularity about these variance 
component estimates that one should be aware of, 
although it does not affect the operation of the system 
in any way. There is a degree of indeterminacy 
about the system in that a constant can be sub- 
tracted from all values of, say, Vz and added to all 
values of, say, V; without altering any of the Brooke 
variances. The only restriction is that none of the 
variance components can ever be negative. Our 
practice is to put the minimum values of V; and V, 
equal to zero, and to put the ‘‘unattached variance”’ 
into V;. This is purely for operational convenience. 
In fact we believe that the minimum value of the 
effect of lateral deviation, which is represented by 

7a is about 1° squared, as indicated in section 5 
above. 

We have now considered the estimation of variance 
components in the first instance. 








The other vital element in an objective bearing 
classification system is a method of keeping the 
system up-to-date. For routine corrections to be 
timely and not unduly subject to sampling fluctua- 
tions, the corrections should be based on as simple a 
statistical model as possible. We therefore assume 
that all variances of bearings taken at a particular 


station are underestimated by a constant small 
amount. The best estimate of this amount is then 
approximately 


M4 (x?7—B,)/Zwi, 


where summation extends over all bearings taken 
at this station. This quantity is called the ‘‘apparent 
correction”, short for ‘‘apparent Brooke variance 
correction”’, and its variance is approximately 4/Zw,. 
If we wait until 2w,>20, the variance of the 
apparent correction is less than about 0.2, so the 
correction is unlikely to be in error by as much as 
one unit. Nevertheless the apparent correction is 
rounded down in absolute value to the nearest unit, 
since some weight should be given to the fact that— 
if the system has been running for some time— 
previous data suggested that no correction was 
necessary. 

In practice this works as follows. For each check 
bearing the “weight” w,; and the “indicator” 
w,(z7—B,) are recorded, and these are summed for 
all check bearings taken in a month at each station. 
These data are inspected and combined with those 
for enough previous months to make 2w,> 20 for the 
station in question. Then 2w,(z?7—B,)/Zw;, rounded 
down in absolute value to the nearest unit, is added 
to all values of V, for this station; i.e., to the value 
of a; in the formula V,=a;+bf. 

After any adjustment to the value of a,, the 
monthly totals of Dw,(z3—B,) must obviously be 
reduced by Zw; multiplied by the addition to a; 
before they are combined with subsequent monthly 
totals to decide whether further changes should be 
made. But it should never be necessary to re- 
compute the B; and w; for individual bearings. 

In practice our estimated variance components are 
not subject to frequent oscillations. 

In addition to these routine adjustments to the 
overall variances for each station, it is important 
that the other variance components should be re- 
viewed from time to time, by computing apparent 
corrections from the formula 2w,(7?—B,)/Dw,; with 
summation extending over all check bearings asso- 
ciated with a particular value of the component 
under review. The desirability of introducing fur- 
ther terms into the expression for the variance can 
be tested in the same way. 

This may be as good a place as any to mention 
the problems of systematic errors and wild bearings. 
Many DF stations exhibit small systematic errors 
in their bearings. These have been studied for some 
time, but corrections to allow for them are not ap- 
plied because the physical basis of the errors is not 
understood. If the systematic error is ignored, one 
is acting as if the error distribution had mean zero. 
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and the Brooke variance is really an estimate of the 
mean square error. 

DF stations are apt to produce a certain number 
of wild bearings that are a long way off the true 
bearing, and which one hopes to be able to reject 
in plotting by their inconsistency with the other 
bearings. Obviously one does not want to include 
these in one’s check bearing data for the determina- 
tion of variances. But the problem is, how large 
must the error be before one 1s justified in rejecting 
the bearing as wild on the grounds that it would 
probably be recognized as such even when the true 
target position was unknown? This is not an easy 
problem. The proper answer presumably depends 
on the variance of nonwild observations, and on the 
number, location, and accuracy of the other DF 
stations in the organization working on the same 
tasks. The rejection of bearing errors greater than 
10 or 15° may often be reasonable. 


Although it has fallen to me to present this account 
of the Brooke system, I hope I have made it clear 
that most of the credit belongs elsewhere. I have 
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been able to refer to some of the relevant work by 
Messrs. Ross, Brooke, and Bowen; but many people 
have contributed in various ways—notably those 
who organize the check bearing program, and collate 
and analyze the data. 


7. References 


Beale, E. M. L., Confidence regions in nonlinear estimation, 
J. Roy. Stat. Soe. (B) 22, 41-88 (1960). 

Beale, E. M. L., The estimation of variances of position 
lines from fixes with unknown target positions, J. Research 
NBS 65D (Radio Prop.) (May-June 1961). 

Bowen, K. C., Sources of error in U-Adcock high-frequency 
direction-finding, Proc. I.E.E. 102 Pt. B, 529-532 (1955). 

Ross, W., The estimation of the probable accuracy of high 
frequency radio direction-finding bearings, Journal I.E.E. 
94 Pt. 3A, 722-726 (1947). 

Ross, W., Lateral deviation of radio waves reflected at the 
ionosphere, DSIR Radio Research Special Rept. No. 
19 (1949). 

Tippett, L. H. C., On the extreme individuals and the range 
of samples taken from a normal population, Biometika 17 
364-387 (1925). 


(Paper 65D3—128) 








BOAT SS RU aE RIA a Je 





OR ESP 





aes 


PRN TRIMER Ee: 


Pa Sete 


é 
oe 
° 
















JOURNAL OF RESEARCH of the National Bureau of Standards—D. Radio Propagation 
Vol. 65D, No. 3, May-June 1961 


Estimation of Variances of Position Lines From 
Fixes With Unknown Target Positions ”” 


E. M. L. Beale 
(August 22, 1960; revised December 27, 1960) 


Formulas are derived for the estimation of the variances of position lines from fixes 
with unknown target positions. Two approaches are considered, (1) that presented by 
Daniels [1951], and (2) an analysis of the squares of the errors in the position lines assuming 
the target is at the least squares estimate for its position. 


1. Introduction 


An important problem in position fixing is to estimate the accuracy of one’s position 
lines. This can be done by assigning a variance to each line, see for example Daniels [1951] 
and Beale [1961]. These variances are usually estimated from fixes on targets whose true 
positions are known to the analyst. There are two major difficulties about estimating vari- 
ances from fixes where the true position is unknown: (1) The estimates depend critically on 
the assumption that errors in different position lines are independent; and (2) The statistical 
problem is difficult, and any valid method seems to involve a substantial computing effort. 

In spite of these difficulties, the problem deserves attention, because it is sometimes 
impossible to set up a satisfactory program of check fixes on known targets. Two approaches 
are considered in this paper. One follows Daniels [1951], the other is based on an analysis 
of the squares of the errors in the position lines assuming that the target is at the least squares 
estimate for its position. 

Section 2 of this paper presents the basic assumptions common to both these approaches. 
The next 3 sections are concerned with Professor Daniels’ approach to the problem. The 
approach is presented in general terms in section 3. Specific formulas for the important 
special case of 4 lines per fix are developed in section 4; and the application of these formulas 
is considered in section 5. The formulas for the alternative approach are developed in sec- 
tion 6, and their application is considered in section 7. 

Finally, in section 8 some suggestions are made for an artificial sampling experiment to 
try out both these approaches. 


2. Basic Assumptions 


The basic assumptions made in this work are: (1) The earth is flat near the true position; 
(2) the position lines are straight lines; (3) an error of observation displaces the line parallel 
to itself; (4) the errors have zero means and are statistically independent; and (5) the errors 
are normally distributed, and we have rough estimates of their variances. 
Given these assumptions, we can take Cartesian coordinates, and denote the j™ position 
line by the equation 
x sin 0;—Y COS 0;=Dpj, (2.1) 





where 6, is a known constant, and p,; is regarded as a random variable with mean 





— sin 0;—7 cos 4), 


1 Contribution from Admiralty Research Laboratory, Teddington, Middlesex, England. : 
2 Paper presented at the Conference on Transmission Problems Related to High-Frequency Direction Finding, at UCLA, June 21-24, 1960. 
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and variance o3, where (é,7) are the unknown coordinates of the true position. In our theo- 
retical work, we assume that £=n=0, so that the mean value of p,; is zero; but we must re- 
member that p; is then not a directly observable quantity. 

Some discussion of these basic assumptions may be helpful. 

The first 3 assumptions are introduced to make the model linear, i.e., to make the value 
of each observation in the absence of errors a linear function of the unknown parameters 
and yn. It is clear that in some DF situations, notably when the DF stations have a narrow 
base line, the problem is decidedly nonlinear when expressed in terms of obvious parameters, 
such as the distances of the target east and north of some origin. But there is evidence to 
suggest that the “intrinsic nonlinearity” as defined by Beale [1960], is nearly always small 
for DF problems. This implies that any results derived from the standard linear theory will 
be valid, provided that they are expressed in terms independent of any particular system 
of parameters (and provided of course that any other assumptions made about the error 
distributions are valid). 

Using the language of Beale [1960], we must define all the quantities we use in terms of 
sample space and the solution locus, and not solely in terms of parameter space (or possible 
target positions on the Earth, which is a representation of parameter space). Now a position 
line can be defined in sample space as the intersection with the solution locus of the hyperplane 
where one coordinate, i.e., one observation, is constant. All we require in this work is an inter- 
pretation of the quantities p,;, 6;, and oj. This can be obtained by taking some specific point 7 
on the Earth near the true target position, and imagining a parameter system that coincides 
locally with distances east and north of T. Now let 8; denote the actual 7" observation. (In 
the DF problem this is the bearing from the 7“ DF station.) Let 8;7 denote a hypothetical 7“ 
observation giving a position line passing through 7, and let p,(@;) denote the (signed) distance 
on the Earth of the j™ position line from 7 as a function of 8). 

Let \; denote d p,/dB;, evaluated where 8;=8)r. 

Then we can write p;=A ;(8;—8;r), o7= Ajvar (8;), and 6; is the angle that the hypothetical 
7 position line passing through 7 makes with some arbitrary z-axis passing through 7. 

Note that it is not necessarily legitimate to interpret 6; as the angle that the observed 

position line on the Earth makes with some z-axis; or to interpret A; as dp,/d8,; evaluated for the 
observed 6;. It is important that the same point 7’ should be used for all position lines in the 
fix. (This approach can be used to derive the standard (Gauss) iterative procedure for finding 
the least squares estimate for the target position when the variances of the 6; are assumed 
known. The point TJ must then be taken as the trial value of the least squares point at each 
stage.) 
Our 4th basic assumption, that the error distributions have zero means and are inde- 
pendent, is very important. In certain circumstaaces it will not be satisfied in practice, in 
which case the least squares point will usually be a less accurate estimate of the target’s true 
position than standard theory suggests, even if the variances are correctly estimated; and 
furthermore the variances will be underestimated by the methods described in this paper. If 
we had numerical values for all correlations involved, we could allow for them, but in general 
there seems to be no practical alternative to assuming independence, hoping for the best, but 
realizing the possibility of being misled by correlated data. 

Our 5th assumption, that the errors are normally distributed and we have rough estimates 
of their variances, is not so critical. It is only required to give appropriate weights in our least 
squares estimation of the variances, and to estimate the accuracy of our final variance estimates. 


3. Daniels’ Approach in General 


The problem considered in this paper was considered by Daniels [1951]. In section 9, 
Daniels writes: 
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‘‘When the true position corresponding to each cocked hat is unknown, . . . it is found that the method of 
maximum likelihood is useless for estimating the variances of the lines. It gives in general inconsistent answers, 
and sometimes no answer at all . . . . The method [breaks down because] the number of incidental parame- 
ters to be estimated increases with the number of cocked hats used. 

‘Nevertheless a simple method of estimating the variances does exist requiring no assumptions about the 
error distributions, provided parallel displacement is assumed and the errors are independent. Suppose there 
are N three-line cocked hats, not necessarily for the same true position, the lines having variances aj, 03, o3 
to be determined. The expression 

U=Pp, sin (6;—63) + p2 sin (6;—6;) + p3 sin (0;—@:) (3.1) 
has the same value whatever the point from which the perpendiculars p;, p:, and p3 are dropped on to the lines, 


and wu is easily measured by taking the point at a vertex. . . . If the point is chosen at the unknown true po- 
sition it follows that 

E(u?) =o7 sin 2(.—63) +03 sin? (@;—6;) +03 sin? (0,—6@2). (3.2) 
If N>3, the N equations with the observed values of u? replacing E(u?) may be solved by least squares. Since 
for normal errors var (u?) =2{ E(u?) }?, arbitrary weights have to be used in the first solution and correct weights 
approximated to in subsequent iterations. The estimates are, however, unbiased for any error distribution 
whatever the weights used.”’ 


But it should be noted that with fixed DF stations the normal equations derived from (3.2) 
to estimate the variances are liable to be very ill conditioned, and may even be singular. If the 
simplifying assumption about parallel displacements were strictly correct, then the values of 
6,, 62, and 6; would be the same for all tasks, so that the right hand sides of all equations of the 
form (3.2) would be strictly proportional, and one could never hope to estimate more than the 
given linear function of oj, o3, and oj. One might hope in practice to be saved by the error in 
the approximation; since @, may be effectively constant for all reasonably possible observa- 
tions on a given target, but not constant over all targets. But this will only be so if the targets 
are well distributed about the DF stations. 

It is fairly easy to see that this difficulty is fundamental to the problem, and is not simply a 
defect in the present approach. For if the 6; were strictly constant, one could for example 
assume that o;=02.=0, and the signed distance of the point of intersection of the first two po- 
sition lines from the third will have a certain probability distribution (normal if the bearing errors 
are normal) with mean zero and variance given by (3.2). This variance is constant as long as 
the true values of oj, o3, and oj are constant. So data of this type can never disprove the 
hypothesis that o;=02.=0. 

In the analogous one-dimensional problem, if instruments 1 and 2 give independent obser- 
vations x, and 2, of an unknown scalar quantity £, and z; has mean é and variance oj, then 
E(x,—22)?=o7+03. With an arbitrarily large number of pairs of observations one can therefore 
get an arbitrarily accurate estimate of oj-+03, but cannot estimate oj or oj separately. Butifa 
third independent instrument is available, then one can estimate oj+03 by (x1—<22)*, oj 
+ 63 by (a,;—23)?, and 03+ 03 by (x.—23)*. Hence of is estimated by the mean value of 3(2;—22)? 
+4(a,—23)?—4(a2—43), i.e., 43 —11%.—L1%3+ X23, as pointed out by Pearson [1902]. 

If one adds a fourth position line to the 2-dimensional problem, then one can obtain four 
expressions analogous to (3.2) by taking each set of three lines in turn. It turns out that these 
four expressions do not suffice to estimate the four variances o7, 03, 03, and 07; because the deter- 
minant of the left hand side of the four equations vanishes identically. But Professor Daniels 
has pointed out in a private communication that one can get more expressions by considering 
the mean values of the products of u-statistics for different triangles. Thus if 


, Ur23= Pi SIN (0:—83)+ D2 sin (6;—6,)+ ps sin (6,—62), 
anc 
Ur24=P; SiN (6.—44)+ pz Sin (6,—0,) +p, sin (0,—82), 
then 
E(uyo3tt24) = 0? sin (0.—63) sin (6.—0s) + 63 sin (@;—6;) sin (@,—#;). 
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FiGurRE 1. A possible set of 4 position lines. 


Attention must be paid to signs. It is convenient to regard each position line as a directed 
line, pointing in the direction @ from some apex at an arbitrarily large distance from the scene 
of action. One can then define the corresponding displacement p as positive for points whose 
bearing from the apex is greater than @ (by an arbitrarily small amount) and negative on the 
other side of the position line. This will ensure that the signs of the terms on the right hand 
sides of the expressions analogous to (3.1) come out as expected. But it will still be necessary 
to study the configuration of lines to see which u-statistics are positive and which are negative. 
In the situation illustrated in figure 1, #3, and Ug, are positive, and U3; and 4 are negative. 

Given fixes with n>3 position lines per fix, one can obtain U=n(n-1)(n-2)/6 u-statistics 
from each fix, and U(U+1)/2 derived statistics by considering the squares and products of 
the u-statistics. 


If n=4, this gives 4 u-statistics and 10 derived statistics. 
If n=5, we have 10 w-statistics and 55 derived statistics. 


This approach is therefore not very practical in its present form for n>5. But the case 
n=4 is important, as it is the smallest n for which unique variance estimates can be obtained 
from a set of fixes, each with essentially the same values of 6;, though not all on the same tar- 
get. We therefore explore this case in more detail in the next section. 


4. Daniels’ Approach With n=4 


With 4 position lines one can form 4 triangles, and obtain 4 basic quantities 2123, Uy24, Uis4; 
and W234, which have the following expressions 
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U123= P1823 + P2831 + P3812 











Uj24 = P1824 P2841 + P4812 4.1) 
U134—= P1834 + P3841 + P4813 on 
x U34= P2834+ P3842 + P4803, 
. where s,;; denotes sin (6;—4;). 
From these one can derive a vector vy of 10 derived observations, such 
that 
I=23 a) 
b where 
e Yi=Ujos, Yo= UrogUi2s, Ys=UrosUliss, Ys= UrogUoas, 
3 Ys=Ujos, Yo= UrosUisa, Y7 = Ujq4Uo34, (4.3 
Ys=Ujs4, Yor UrsiUoaa, am) 
Y10= U5s43 
C=Pi, 2=PiP2, Ws= PDs, G1=PiPs, 
=P, de=P2P3, I7= P2P4; (4.4) 
ds=P3, Yo= PsPs, 
qio= Pi; 
and A is 
[82s — 2813823 2812823 0 Si3 —2819813 0 st, 0 0 7 
823824 (—t,—ts) 812824 $12803 813814 —$12814 —Sy2S13 0 S19 0 
823834 —S813834  (ti— ts) 813823 0 $13814 — 8i, — 812814 812813 0 
0 823834 —823824 833 —813834 (ti tt,) —S813823 —812824 812823 0 
Sou — 2814804 0 2819804 Sie 0 —28)0814 0 0 Si, 2 
824834 —814834 —S14504 (4, +#,)0 Sis —813814 0 — 812814 812813 (4.5) 
0 $24834 — 84 823824 —— 814834 §14824 (t;—ts) 0 — 812824 812823 
sz, 0 — 2814834 2513834 0 0 0 87, —2813814 Siz 
0 S34 — 824834 82383, 0 — 814834 813834 $1 4824 (—t.—ts) 813823 
9 0 0 0 834 — 2824834 2823834 834 — 2823824 833 J 
where t, = 812834, f2= 813824 and t3= 814823. 
Then 
Ey) =H(a/’, 03, G3, 0%)’, (4.6) 


where H is the matrix formed from the Ist, 5th, 8th, and 10th colunms of A, since the p; are 
independent and have zero means when referred to the true target position as origin. 

But in order to determine rational, as opposed to arbitrary, least squares estimates for the 
unknown parameters oj, ¢3, 3, and oj, it is necessary to consider the covariance matrix of y, 
We rely on the fact that the p, are independent and have zero means, and find that 


Cis= COV (Ys) = Ly iy) — (Ey) (Ey) = Didar, (4.7) 
where 

"=e, l= 0{03, 1'3= 9003, r= 010%, ) 

rs=204, re=o}03, roe, 
| (4.8) 
. 5 Nei Oo ee Sa 
% ’g=203, To= 9303, 
S r 4 

Tigp= 2034. J 









(The coefficients 2 in 7, 7s, ’s, and 74) should be replaced by 2+-«, if the error distributions have 
a fourth cumulant «, different from zero.) 
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The covariance matrix C=(e,;) of the derived observations is therefore obtained from 
(4.7), or equivalently by defining a new matrix 
F=(f,), 


where f,;=a1;rj?, 


so that C=FF’, 

Now the standard least squares results can be expressed in matrix notation as follows: 
If Y=X-+e, where Y is a vector of observations, X is a known matrix, 8 is a vector of unknown 
parameters, and ¢ is a random vector whose components are uncorrelated and have means 
zero and a common variance o’, then the least squares estimates B of 8 are given by choosing 
B to minimize 

Q=(XB—Y)’(XB-Y), (4.9) 

Assuming that the errors are normally distributed, the likelihood is proportional to 

exp (— Q/20’), and in any case 


B=(h's)“*2’T, (4.10) 
and the covariance matrix for the estimates B is given by (X’ X)~'!o?. (4.11) 


If now y=GY, we have y=GX6+Ge, and therefore the covariance of y; and y; is given by 


DIG no. 


i 
Hence we can identify G with the matrix F, putting c’-=1, and we know that, if B=(o7, 
2 2 27 
02, 93, gi) ? 


y=H¢6+ Fe. 


~ ~~ ~ 


and 
Q=(F-\(Hg—y) }’{E-(H—y)} 
=(Hg—y)’F-F"(ig—y) 
= (Hg—y)’(FR”)""(Hs—y) 
=(Hg—y)’C-'(Hg—y). (4.12) 
Further 


p— UE EH) WE Ey 
=(H'CH)H'Cy, 4.13) 
and the covariance matrix of 6 is given by 
(H’C-"H)-". . (4.14) 
These results are due to Aitken [1935]. 
5. Application of the Formulas of Section 4 


To apply these formulas to a situation in which one has N fixes, i.e., sets of four position 
lines, each with a constant variance and a constant value of 6,;, one proceeds as follows. 
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One first sorts out one’s sign convention, as indicated at the end of section 3. This will 
produce signed values of 2123, Ui24, W134, ANA 34 for each of the N fixes. From these one computes 
derived observations y;, . . ., Yio for each fix from (4.3). Averaging over all fixes, one obtains 
a vector 

y=(, me - - + Sar. 


We next compute the matrix A from (4.5). 

By using the guessed values for the unknown variances oj, o3, 03, and oj, one can derive 
values of 77, 73, . .., Tio from (4.8), and hence compute the elements c¢;; of the matrix C 
from (4.7). : 

This matrix must be inverted to produce C™. 

We must also consider the matrix H whose columns are the Ist, 5th, 8th, and 10th columns 
of A, and our estimates of of, 03, 3, and o? are the 4 elements of the matrix 

(H’C'H)“'H’C“'y. (5.1) 

If our guessed values of the «7, used to define r, were correct, the covariance matrix for 

these estimates would be ~ 
1 
N 

One could try iterating this procedure, using the estimated variances to produce a new 
vector r, and hence a new C, and hence revised estimated variances, and so on. 

I am doubtful of the wisdom of this, though it might be interesting to try it on an artificial 
sampling experiment. If one wants a reasonable estimate of the accuracy of the estimates one 
must compute a revised C to use in (5.2). But I can see little merit and some danger in revising 
the weights to produce revised estimates: Note that the iterated estimates are not unbiased. 

This procedure may on occasion produce negative estimated variances. Steps must be 
taken at least to ensure that negative (or even very small positive) variances are not used in 
the definition of r. Theoretically the least squares problem should be set up as a quadratic 
programing problem when such negative estimates are produced by the standard method. 
But this is too large an issue to face at this stage. 


(H’C~"H)~". (5.2) 


6. A Direct Approach 


In this section we derive some equations concerning the apparent error in the position 
lines, assuming that the target is at the least squares estimate for its position. 
Suppose that the j7™ position line is 
x sin 0;—y cos 0;=p,, 


with the origin at the (unknown) true target position. Suppose further that we have some 
rough estimate sj for 05. 
We write 
o3=s7?(71+6,), and 1/s?=w,;. 
Now the sum of squares function, whose minimum value defines the least squares 


estimate, is 
Q=D5w,;(z sin 0;—y cos 6;—p;)?=a #?+2hary+ by’?+2gr+2fy+e, 
j 


where 
. 9 
a=) w; sin? 6;,, b= > jw; cos? 0,, c= > )wyp}, 
j j b 


; : (6.1 
f= 20s cos 0, g=—D wy p; sin 6,, h= — 210s sin 6; cos 6;. ~_ 
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Note that a, 6, and h are known constants, while f, g, and e¢ depend on the random 
variables p;. 
It is convenient to rotate axes, so that in the new coordinates (X, Y), 


Q=AX?+ BY?+2GX+2FY +. 
This involves finding an angle a such that 


$;=6;1 @, (6.2) 
and the jth position line is 
X sin ¢;—Y cos $;=p,, 
i.e., 
X sin (0;+a)—Y cos (0;+a)=py;, 
so that 
z= X cos a+ Y sin a, 


y=—X sin a+Y cos a. 


And since axz?+2hry+ by? +2g¢+2fy+e=AX?+2HXY+BY’?+2GX+2FY+C, where H=0, 


we have 





A=a cos? a—2h cos a sin a+b sin? a, : 

0=a cos asin a+h(cos? a—sin? a)—b cos asin a, 

B=a sin? a+2h cos a sin a+b cos? a, q (6.3) 
° 3) 

G=g cos a—f sin a, 

F=g sin a+f cos a, 

C=e. | 

This implies that 

tan 2a=2h/(b—a). (6.4) 


This defines a set of values of a differing by multiples of 90°. It does not matter which 


is taken. 
The values of A, B, C, F, and @ can be obtained either from (6.1) and (6.3), or alter- 
natively from the formulas 


A=) jw; sin? ¢;, B= > jw; cos’? ¢,, C= Dd wypi, 
j j j 


F= 20s cos oy, G=— twp; sin $y), H=— 220; sin $; cos ¢;=0, 6.5) 
where ¢; is obtained by substituting for a from (6.4) in (6.2). 
The transformed coordinates of the least squares estimate are then 
Xo=—G/A and Y»>=— F/B. (6.6) 


Now let d, denote the signed displacement of the k” position line from the least squares 
estimate. This quantity can be observed in practice. Then 








G. F " ‘ 
h=PtZ sin oj— B cos ome” > Pi (w;/wy) re, (6.7) 
where 
sin? dy 2 
y= 1— wy (SB OSH), (6.8) 
Nn=Au=— (wyw,)”? (SEB: fn Op COS 0008 or, (iXk). (6.9) 
Note that 
n 
2 A=N— >, sin? d,/A—>)w, cos? ¢/B=n—2. (6.10) 
=1 
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Further 











Be sin? ¢, i w; Sin? g; sin? dx W; SIN gd; COS ¢; SIN ¢; COS ¢; 
2,—1—2w, be 2 — : 
aN me ( a pe tw, > ping AB 
; COS* d; Cos? sin? s 
42 s $8 ox \ polcm (= nt tr, 
So 
n . 
24 Nix =x. (6.11) 


Further, if kl, 


= SII ; Sin COS ¢; COS sin ¢; sin OS g; COs 
Do Nadura 2der (wy)? 9 w(° oe de 4 4 OS 9; 08 *\(= 1 di —- COS d sete 





z B 
iy tne” [= gd; SIN ty 008 dx COS *) from (6.5) =x. 
. A B 
So (6.11) can be generalized to read 
b> NiAir=An- (6.12) 
i=1 
From (6.7) we deduce that 
w.di= => > (wyw,) ried jx Pi Pj- (6.13) 
Now 
Epi=o7= (1+6,)/w;, Ep;ps=0 (i¥)), 
So 
Ewdi= Nix (1 +6;) =A+>d_ 6;A3,, from (6.11). (6.14) 
i=1 


It follows that 
E p» wedi =>) hut 6; Ds Men 24+ 39 bi, from (6.10) and (6.11). (6.15) 
If 6,=6 for all 7, we have the well-known formula 
E SSwd?= (n—2)(1+8), (6.16) 
(which can easily be deduced from first principles). 


It is also of interest to consider the covariance of w,d? and w,d?. 
We have 


n 
cov (wedi, wid?) =D) Wididi,2(1+-6,)?/wi 
i=1 


+2 P Dy WW Aidpdiadji(1+6;) (1+6,)/ww; by analogy with (4.7) 


=1 jx¥i 
n n 
=2 ym >» Nid peAitAj,(1 +6;) (1+6,). (6.17) 
i=l j= 
It is of particular interest to consider the covariance matrix when all 6;—0, since this 


corresponds to the originally estimated covariance matrix. 
Using (6.12) we find this reduces to 


cov (wedz, wd?) =2n3,. (6.18) 
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If n=3, it can be shown that j4;/A2}=Ajo/A22—=Aj3/Ag3, ete. For example, to prove that 
AyyA22— AZo =0, one substitutes for the ,, from (6.8) and (6.9) and multiplies out. The terms 
in 1/A? and 1/B? cancel, so the expression can be written as a quadratic function of w,,w, and 
w; (using the definitions of A and B) divided by AB, and H proves to be a factor of the 
numerator. 

This is an algebraic proof of something we already know on statistical grounds: that with 
only 3 bearings the random variables w,d?, w.d3} and w,dj are totally correlated, and it is 
impossible to estimate more than one linear function of the unknown variances. 


7. Application of the Formulas of Section 6 


To apply the formulas of section 6 to a situation in which one has N fixes, i.e., sets of n 
position lines, each with a constant variance and a constant value of 6;, one proceeds as follows. 
One starts with some assumed variances sj for the position lines, and one uses those to 
compute a least squares point estimate for the target for each fix. One then computes w,d?, 
i.e., dj/s;, being the weighted square of the displacement of the 7 position line from each fix. 


One then averages over all N fixes, obtaining n quantities w,d?, . . ., w,d?. 
One must also compute A;; for 7=1,. . ., n, and dj; for 7, 7=1,. . .,n. These quantities 


are the same for all fixes. 
Note that d,; is defined by (6.8) and A,; by (6.9) if 747. 
It may be more convenient to use the formula 








2 __(W, Sin? d;)(w; sin?oy) 2(w, sin ¢; cos ¢,) (w; sin ¢; cos $,) 4 (w; cos" ;) (W; cos’ $;) » (i%j). (7.1) 


“ me AB B 


Then, from (6.14), we solve the equations 
na —_—_— 
> 36 A7,= Wd? — x, (e=1, cosets n). (7.2) 
i=1 


and the resulting values of 6; define the estimated variances, since 
of =8,7(1 + 6;). 
To find the estimated covariance matrix for these estimates 6;, we form the matrix H 
such that 
(hih=Mi, 
and the matrix C such that 
(e4)=2,,;?_ (from (6.18)). (7.20) 


The estimated covariance matrix is then given by (5.2), if the 6; prove to be small, and 


; 4 
this reduces to => C7}. 


N ~ 
If one wants to combine these data with other data, then the contribution of these data to 
the sum of squares to be minimized is then NQ, where Q is given by (4.12) with B=(6,, . . ., 








6,)’ and y=(w,d,", . . . , Wade’)’. 


8. Proposed Computational Program 


Some computational experience with both the schemes proposed here would be very 


valuable. 
In the first instance, it would be of interest to inspect the covariance matrix (4.14) for 


estimates obtained by Daniels’ approach, i.e., with H defined following (4.6) and C defined 
by (4.7), for various values of 6; and o,?. 
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It would also be of interest to inspect the corresponding matrices for the direct approach, 
i.e., with (hj;) = (Aj?) and (¢;;)= (2,77). In this case (4.14) reduces to4 C~!. The (77) element 
of this covariance matrix for the 6 must be course of multiplied by ses? to represent the covari- 
ance matrix for the variances. It will then be directly comparable with the other. 

One may hope to throw considerable further light on these procedures by an artificial 
sampling experiment. One must first choose values of 6;, o;?, and (for the direct approach) 
s;. These will define values of a, 6, and A from (6.1), and hence of A, B and a from (6.3) and 
(6.4). One can then compute the d,, from (6.9). 

One then takes sets of values of p; as pseudo-random normal deviates with means zero 
and variances o,7, and use these to generate synthetic w-statistics from (4.1), and values of d; 
from (6.7). One can then compare the estimated variances obtained by the two approaches 
(a) after one iteration, and (b) after using the results of each iteration as starting values for 
the next (i.e., to compute the 7; for Daniels’ approach and the w; for the direct approach) 
and continuing until the variances that go in come out. 

One can also consider the effects on the two approaches of using poor starting approxi- 
mations to the variances. 

It seems likely on general grounds that Daniels’ approach should give better results, since 
it uses more information. But if the direct approach is at all satisfactory with n=4, one may 
hope that it will be still better with larger values of n. 
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Statistics of a Radio Wave Dittracted 
by a Random Ionosphere" *’ 


S. A. Bowhill 


(August 15, 1960; revised December 27, 1960) 


For some purposes, particularly in connection with the study of the random structure 
of the lower ionosphere, using very low frequencies, it is necessary to find the detailed 
statistical properties of a random signal diffracting in free space. Mathematical tools for 
evaluating these parameters have been developed, and are applied in this paper. Allowance 
is made for the effect of sphericity of the wave incident on the ionosphere, and anisotropy 
of the irregular variations of signal is permitted. The case of oblique incidence of a wave 
on the ionosphere is also considered. 


1. Introduction 


There has of late been an increasing interest in the irregular structure of the lower iono- 
sphere, particularly in connection with scatter propagation. In order to interpret irregular 
fluctuations in radio signals observed at the ground in terms of irregularities in the ionosphere, 
it is first necessary to find how much the nature of the signal changes as it propagates from the 
ionosphere to the earth. This propagation, which takes place in free space, will be termed 
“diffraction,” in the sense used in a previous paper [Bowhill, 1957]. 

Most of the methods so far proposed for dealing with this statistical problem assume that. 
the signal at the layer is like the signal produced when an infinite plane wave passes through 
a thin diffracting screen. The field at the screen is then resolved into a number of plane waves 
propagating in various directions, their amplitudes being given by the Fourier spectrum of 
the field variation. The signal at any distance from the screen is then given by the vector 
sum of the signals from the various plane waves, combined with the appropriate phases. 
Booker, Ratcliffe, and Shinn [1950] stated this principle, but Hewish [1952] was the first to show 
quantitatively how the magnitude of the irregular field variations changed with distance from 
the screen. He assumed a rectangular distribution of spatial wave numbers in the screen. 
Feinstein [1954] has used a rather different approach based on Huygens’ principle, which has 
also been exploited by Chernov [1960]. Ratcliffe [1956] has summarized the physical prin- 
ciples involved in these approaches. 

Pitteway [1958] has recently given an analytic method for calculating the field due to any 
particular assembly of irregularities; this, however, is difficult to apply in a statistically random 
medium. 

It is important to distinguish at this point between the various types of random screen 
postulated. The two types most commonly considered are screens which modulate the phase 
or the amplitude of the signal. A phase screen is said to be “deep” or “shallow” depending 
on whether the phase excursion of the signal is substantially larger or smaller than one radian. 
The deep phase screen problem has been approached by Hewish [1951] by using a sinusoidal 
variation of phase with distance, but the extension of a detailed theory to this case presents 
considerable difficulties. The shallow phase screen, on the other hand, is exactly equivalent 
to a shallow amplitude screen [Bowhill, 1959]. This interchangeability of phase and amplitude 


1 Contribution from Ionosphere Research Laboratory, The Pennsylvania State University, University Park, Pennsylvania. 
2 The research reported in this paper was sponsored by the Geophysics Research Directorate of the Air Force Cambridge Research Center, 


Air Research and Development Command, under Contract AF 19(604)-1304. 
3 Paper presented at the Conference on Transmission Problems Related to High-Frequency Direction Finding, at UCLA, June 21-24, 1960. 
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is exploited in the phase-contrast micrescope, in which variations in phase are changed into 
variations in amplitude by artificially retarding the phase of the specular component by 7/2. 
It is well known that this only gives an improvement in visibility when the phase variations 
are small (i.e., the phase screen is shallow). 

The results of the analysis in this paper will be applied principally to the random diffrac- 
tion of very low frequency radio waves. The fading for these frequencies is usually shallow 
indicating that a shallowly modulated random field at the ionosphere must be responsible. 
Attention will therefore be confined to shallowly modulated screens. 

To illustrate the method of analysis, it will first be applied, in section 2, to a one-dimen- 
sional random screen. The extension to two dimensions is made in section 3. Section 4 deals 
with the correlation between amplitude and phase in the signal, and section 5 with the effect 
of wave sphericity. The work is extended to oblique incidence in section 6. 


2. One-Dimensional Diffraction Problem 


In treating the spatial correlation function and spatial frequency spectrum of the signal 
variations, it is much more convenient to use purely real quantities, rather than the complex 
use of complex Fourier transforms to re- 


variables of amplitude and phase. This enables the 
angular spectrum of plane waves which 


late field variations over the screen to the complex 


they produce. 
For an amplitude screen, which alone will be considered (in view of the equivalence 


indicated above), the electric field # at the screen is given by 
E(x, 0) =A (x) 

where A(z) is a purely real quantity. The field at a distance z from the screen is then defined 
to be 

E(x, 2)=[Ei(2)+7E,(2)] exp (24iz/d) (1) 
where £,(x) and F,(x) are both purely real, and represent the phase components of E(z,z) in 
phase and in quadrature with the undiffracted portion of the incident plane wave. For a 
shallow screen, F(x) gives the amplitude variation, and F2(r) the phase variation. In the 
analysis which follows, various statistical parameters of the quantities ,(x) and E,(2), which 
are actually observed in experiments at ground level, are evaluated in terms of the corre- 


sponding parameters of the original amplitude variation A (a). 
Using the notation of Booker, Ratcliffe, and Shinn [1950], we have for the angular 


spectrum at the screen, in terms of s=sin 6, where @ is the scattering angle, 


>»(s) = ¥ A(x) - exp (—2mrisay/dr) dx (2) 


v~—@ 


where \ is the wavelength. The diffracted field E(z,z) is given by the Fourier transform of 
this, together with the phase factor exp (27icz/A), where c=cos 6. 


Et,2)=5 | exp {22(sx+cez)/A}ds + { A(a) + exp (—2risxo/A)d zp 


JJ. A(%) - exp {2nt [(x—4)s+e2]/A }dayds. (3) 
By the definition of /,(x) in equation (1), 


Ba) =A{ 5 [ - A(%) - exp {2ri [(x—2)s— (1—e) 2Ijdrais \ 


=x i. A(a) [exp {27 [(x—a)s— (1—e) z]/A} + exp { —2xi [(x—a%)s— (1—c) 2]/A}] dads 
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since A(z) is purely real. The spatial spectrum P;(v) of ,(z) is now defined by a similar 
expression to (2): 


i ; l | ad ae : 
P= | E(x) «exp (—2niva/\)d2=3~ | {| A(x) - [exp {222 [(a—2)s—vae— 2z(1—e)]/X) 


+exp {—2zi [(x—a)s+va—2(1—e)]/A} |daodads. 


Since the integrand is a well-behaved function, and all the limits of the integration are 
—o to +, the integration can be carried out in any order. In fact, it is convenient to in- 
tegrate with respect to first s and then z. With the approximation (valid for small s) that e= 
1—s?/2, P;(v) is easily evaluated as 
> l . feé M | 2/9\ /\ 1 I ‘ of 2/6 ) 

PiQv)=5 A(ao) [exp { 2ri(—vao+2v?/2)/X} + exp {212 (—vao— zv?/2)/d} |dxo 


=cos(rzv?/n) | A(x) - exp (—2mivry/dA)dxo 


or 
P,(v) =P (v) - cos (rzv?/d) (4) 
from equation (2). This equation gives the spatial spectrum of /\(x) at any distance z from 
the screen. The power spectrum W,(v) of this variation is defined by 


Wi) =Pi(?) - PT?) (5) 


the asterisk denoting the complex conjugate, and the bar denoting averaging of the result 
‘“‘over systems”—i.e., the mean for a number of statistically similar screens. This quantity 
is essentially real, and is a relatively smoothly varying function of v, because of the process of 
averaging. It can be related by equation (4) to the spatial power spectrum W (vy) at the screen: 





W1(v)=Po(v) - P§(v) - cos? (wzv?/d) 
or 
Wi (v) =Wo(v) - cos? (wzv?/d) (6) 
where 


Wo(v)=Polv) - Po). 


Similar quantities P,(v) and W,(v) can be defined for the quadrature component EF, of the field 
(see equation (5)); a similar analysis to the preceding gives 


P;(v) =Po(v)-sin (w2v?/d) a 
W,(v) =W (rv) -sin?(azv?/d) @) 

Equations (6) and (7) enable the spatial correlograms of /,(x) and £,(z) in the diffrac- 
tion pattern to be determined. Assuming that the correlogram of A(z) is given, the Wiener- 
Khintchine theorem can be used to give W,(v). The power spectrum W,(v) is then calcu- 
lated from equation (20), and the correlogram of F(x) found by the inverse Wiener-Khintchine 
theorem. 

This one-dimensional analysis, in which the random variations in signal are taken to 
occur in the z-direction only, represents an unrealistic physical model. However, as some 
existing work [Hewish, 1951, 1952; Jones, Millman, and Nertney, 1953] has assumed this 
type of variation, some results worked out for this case will be quoted in the next section, 
for comparison with the two-dimensional results. 
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3. Two-Dimensional Diffraction Problem 


In this section the diffraction pattern of a two-dimensional random screen is considered | 
using an analysis closely similar to that employed in the previous section for the one-dimen- 
sional case. The modulation is once again assumed to be in amplitude only, as there is again 
an equivalence between shallow phase and amplitude screens. The effect of anisotropy of 
the screen will be taken into consideration by defining a spatial power spectrum W)(v,x), 
related to a two-dimensional spatial correlogram Po(a,8) of the field A(z,y) at the screen by 
a two-dimensional form of the Wiener-Khintchine thereom. As the field A(a,y) is purely 


real, the correlogram po(a,8) is given by 





A(x,y) A(x+a, y+8)—A(ey)" 


po(a,8) = A*(x,y)—A(x,y)’ 





The field A(z,y), the mean field over the whole screen, has a constant phase at all points since 
A(z,y) is real, and propagates as a single plane wave. It is convenient to imagine it removed; 
A(z,y) then has both positive and negative values, and will be defined to have a mean square 





value of unity; so 


A(x, Yo) =0 
(8) 


A?(%,Yo)=1) 


and the correlogram is given by 





po(a,8) =A(2p, Yo) A (%o+a, Yot B) . 


The single plane wave, which will be called the ‘‘specular component” (by analogy with 
specular reflection from a mirror, compared with irregular reflection from a rough surface), 
can be reintroduced when necessary into the diffraction pattern as a constant added to F,(z,y). 

Proceeding exactly as in section 2, the two phase components of the diffracted field are 


defined by the vector field at (x,y,z): 
Ea, y,2)=(A@,y) + 122(z,y)] exp (2miz/n). (9) 


The phases of the diffracted waves at z are given by purely geometrical considerations. Let 
the plane of one wave intersect the (z,z) plane at an angle sin~'s, and the (y,z) plane at an 
angle sin~'r. Then the angle between the wave front and the (z,y) plane is 


tan“!{ (s-?—1)-!4 (r-?—1) 3}? 
and the phase of the wave at a distance z from the screen is 
Qriz{1+(s-?—1)7!4+ (r-?—1)7}-¥), 
If s and r are small, this can be expanded to give, to a first order, 
Qriz(s?/2+7r?/2)/r 
The corresponding equation to (3) for £(z,z) is therefore . 


E(z,y,z) exp (—2iz/), 


= {ff i A (a, Yo) «exp {2ri[(x— 29) 8+ (y—Yo)?— 8? 2/2—r? 2/2] /A} ds dr dao dyo (10) 


and by similar steps to the previous analysis, it can be shown that 
P,(v,X) = Po(v,x)-cos{ r2(v? +x?)/d} (11) 
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where Po(v,x) is the two-dimensional spectrum of the field at the screen. Similarly 
P,(v,X) = Po(v,x)-sin{ w2(v? +x?) /r} 
and the spatial power spectra of /\(z,y) and (x,y) are given by 
W, (v, x) =W, (», X).cos?{w2(v?-+x?)/A} 


(12) 
Wz (v,x) =W (v,x).sin® {w2(?+2") /r}. 


These spectra are proportional to the area of the diffracting screen; this will be shown below 
in deriving (16). The spatial correlogram p,(a,8) of E,(z,y) is given by the usual relation 





Ey (2,y) FE (a+a,y+ 8) n 
p (a8) = (13) 
EX (x,y) 
since the mean value of /,(z,y) is zero. On comparing this definition with (9), it will be 
seen that (13) has the term /}(z,y) in the denominator, whereas, from (8), A?(%,Yo)=1. In 
fact, the mean square fluctuation in A(a9,yo) at the screen appears, at some distance away, as 
a fluctuation partly in F,(z,y). It is easily shown from (12) that 





43 (2,y) +3 (2,y) =A?(%,Yo) = 1, (14) 


the sum of the mean square values of /\(z,y) and /,(z,y) remaining constant at all distances 
from the screen. These mean square values represent, in the case of an amplitude screen, the 
amount of fluctuation of the amplitude and the phase respectively at a distance z from the 
screen. 

The autocovariance of /, can be calculated from the spectrum W,(y,x) in (12) by the 
Wiener-Khintchine theorem; /}(z,y) is then its value for a=8=0. Before this can be cal- 
culated, the power spectrum W,(»,x) of the signal A(2o,yo) at the screen must be found. Some 
assumption is needed for the spatial characteristics of the signal. They will be assumed to 
be as follows: 


i. A(z,y) has a Gaussian correlogram in all directions. 

ii. The contours of the two-dimensional correlogram are elliptical. 

iii. The ellipses have their principal axes oriented along the z and y axes of the co-ordinate 
‘system. 

The p=0.61 contour on the two-dimensional correlogram is used to define two structure 
sizes D, and D, of the screen in the z and y directions. The correlogram of A(z,y) is then 


po(a,8) =exp { — (a*/2D}) — (6?/2D3) }. (15) 
The angular spectrum at the screen is given by an equation like (2) 


Pus)=|{ A(a,Yo) exp {—2ri(sxy+ryo)/A}daodyo 


and the power spectrum is given by 





W, (s,r) =P,(s,r) ‘P3 (s,r) =—=mean ot { ff A (Xo, Yo) -exp { — Dri (Sa+ryo)/A}dapdyo 
xf [" AGx)- exp (2ni(sxi-+ryi) Ada } 


= | J J J A (ato Yo) +A (29,40) exp { Qmi[ (1% — a0) 8+ (Yo— Yo)? |/} dad yod rod yo 
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where 2 and yp are variables of integration, analogous to x and y, but differentiated from 


them to avoid confusion. So 
Wels =f[ [J — AGoy)-Alo+ ayo+8).exp (2ri(as+Br)/A}dzalydla dB 


where a=2)—2%, B=Yo—Yo. By (14), it is evident that 





A(20,Yo)-A(to+ a, Yot B)=po(a,B) 
and therefore 


woer= fff ; po(a,8)-exp{2ri(as+ Br)/dA}dxodyoda dB. 


Now neither zp nor y) appears in the integrand; so the integration with respect to 2» and yp will 
give a W,(s,r) which is proportional to the total area S of the diffracting screen: 


W,(s,r)=8- | | po(e,8)-exp {2mi(as-+Br)/d}dadB. (16) 


Substituting for po(a,8) from (15) and integrating, 
W(s,r) =2a2SD,D,.exp { —22?(D,?s? + D,?r?)/d} - (17) 
The spatial power spectra of E,(z,y) and E,(z,y) are found by substituting from (17) into (12): 


W,(v,x) =2aSD,D,.exp { —22°(D,?v? + Dz? x?) /d?} .cos?{ rz(v? + x?)/d} } (18) 
W2(v,x) =22SD,D?.exp{ —22?(D,2v? + Dz2x2)/d?} sin? { r2(v?+ x?)/d} " 


It is now possible to find the numerator of (13); the covariance of E,(z,y) and EF, (x+a,y+8). 
The inverse transform corresponding to (16) is 





E\ (a, y) E\(x+e, VT D=—y \{ : W,(v, x) - exp { —2r7(va+xB)/d}dv dx (19) 


with the value of W,(v,x) given by (18). The integral is evaluated by putting 
cos*{ r2z(v?+x?)/A} = 1/24 (1/4)exp{ 2mriz(v?+x*)/dr} + (1/4) exp{ —2aiz(v? +x?)/d}. (20) 


As (19) is the inverse of (16), the term 1/2 gives simply (1/2)p9(a,8). The two remaining 
terms in (40) are identical except that the sign of z is reversed. The integral can be simplified 
further by separating completely the integrals with respect to v and x. This gives 


E\ (a, y) THetayt pa: {I(a, D,, z)I(8, D2, 2) +I(a, Di, —2)I(8, D2, —2)}+ (1/2) po(a+8) 


(21) 





where 
I(a, D,, 2) -| exp {2mi[av+ (2+7iD}/d)v?|/d} dv 


=i ta)" exp {—(1/2) tan! a1—(1—iai)a?/2D4(1+a8)} (22) 


V2nD, 
where the important substitutions 
. AZ . AZ 
=->3) |) — 
DP? DP 
have been made. These parameters will appear as the measures of distance on all the graphs 
depicting the behavior of the various statistical parameters of the diffraction pattern. 
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Finally, on substituting from (22) and (21), 





Ey (2,y)- Ey (2+ a,y+B) = (1/2)exp{ —a?/2D,?— B?/2D,"} 
+ (1/2) (1+ a,?)—!4(1+4,”)~“exp{ —a?/2D,?(1 + a,?) —8B/2D3(1 +43) } 
X cos { aya?/2D,?(1 +4,”) + 426?/2D,?(1 + a,?(1 + a2?) — (tan~'a,+ tan~!a2)/2}. (23) 
By putting a=6=0 | 
EP(a,y) =1/24+ (1/2) (1+ a,2)-"*- (1+ a,2)-*-cos{ (tan-ta, + tan~!a,)/2}- (24) 
# 
The correlogram p;(a,8) is given by the ratio of (23) to (24). The expressions for the auto- 
covariance of (x,y) are identical, except that the positive signs, denoted above with asterisks, 


become negative. Some special cases will now be investigated. 
(i) Near the screen a,=a,=0, and 





FE, (2,y)- Ey (2+ a,y+ 8) = exp(— a?/2D,?— B?/2D,*) = po(a,8) 





E,?(2,y)=1, E.?(x,y)=0. 


as would be expected, since E,(x,y) = A(2,yo) at the screen. 
(ii) Very far from the screen Qy,.>©, tan~'a,=tan~‘a,.=7/2, 
and 





Ey (a,y) E\(a+a,y+ B) = (1/2)exp{ —a?/2D,?—B?/2D,?} 





t1°(x,y) = be*(x,y) = 1/2 
pi(a,8) = p2(a,8) =exp{ —a?/2D,°—6?/2D,’}. (25) 


The second result in (25) means that the space correlation of both amplitude and phase at a 
great distance from an amplitude screen is just the same as the correlation at the screen. This 
is not true at intermediate distances, as we shall shortly see. The first result in (25) shows that 
the fluctuations occur equally in amplitude and phase under these conditions. The phase- 
amplitude locus of the instantaneous signal at large distances therefore forms a circular ‘“‘cloud” 
around the vector denoting the mean carrier signal, as in figure 1. 

(iii) Lsotropie screen D,=D,.=D, and a,=a.=4. 





Ey(x,y) Ey (e+ ay +B) = (1/2)exp{ — (a+ 6") /2D*} 
+ (1/2)(1+@?)-"?-exp{ — (a?+ 6?) /2D?(1 +a") } -cos{a(a’+ B*)/2D*(1+a*)—tan~'a}. (26) 


The quantities a and 8 appear in this expression only as (a?+*), so the contours of correlation 
are circles at all distances, and the diffraction pattern is never anisotropic. Also 


E,*(2,y)=1/2 t (1/2)(1+a?)-! 


LOCUS OF RANDOM SPECULARLY 
COMPONENT AS ONE { REFLECTED Figure 1. A possible phase—amplitude diagram for 
ae eee —— the diffracted signal a long way from a random 
screen. 
A E2(x,y) 
— £, (x,y) 
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the asterisk signifying, as previously, the sign which is to be reversed to give the corresponding 
value for E3(z,y). H?(x,y) and £,?(z,y) are plotted on figure 2. It can be seen that the mean 
square values of E, and EF, are virtually equal at distances from the screen greater than that 
corresponding to about a=4. It will be shown later that, under these circumstances, the 
diffraction pattern is not correlated with the signal at the screen. 

Slight rearrangement of (26) gives 





Ey(z, y) Ex\(w+ea, y+8)=(1/2). exp{—(e?+ 6*)/2D"} 





4+ (1/2)(1-+42)-!, exp{ —(a?-+6?)/2D*(1-+a?)} [cos{a(a?-+ 6*)/2D*(1 +a?) } 
+ asin {a(a?+?)/2D?(1+a?)}]. (27) 


(iv) One-dimensional screen D,=D; D,>@. So a-—0. 
From (24), 
E(x, y)=1/2+ (1/2) (1+4?)-“*-cos{ (tan a)/2} 
* 


=1/24(1/2){11+ (1+a2)"]f1 +a7}-1/2}", 


This also is plotted in figure 2, for comparison with the curve for the two-dimensional case. 
Evidently the mean square values of F,(2,y) and E,(2,y) become nearly equal much closer to a 
two-dimensionally irregular screen than a screen which is irregular in one dimension only. 
This illustrates that a one-dimensional model may give a quite inadequate representation of two- 
dimensional diffraction phenomena. For this case, 





Ey (a, y) E,(a+a, y+8)=(1/2) exp(—a2/2D?) 
+ (1/2)(1+a?)-"* exp { —a?/2D?(1+a?) }-cos{aa?/2D?(1+<a?) 
* 


— (tan! a)/2}. 
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It will be noted from (27) that the correlogram is not Gaussian in form, though /(x,y) has a 
Gaussian correlogram for a=0, and both F,(z,y) and £,(z,y) have Gaussian correlograms as a 
—>o (see equation (25)). This introduces some difficulty as to what is to be called the “struc- 
ture size” d of the diffraction pattern, for comparison with experiment. The following defini- 
tion is adopted 


d=Structure size in x direction = (—0°p/0a’) 71/7. 


=D, for the correlogram po(@,8). 


This definition bases the structure size on the Gaussian curve which is the best fit to the correlo- 
gram, near the origin of a. It can best be found by expanding (23) in powers of a” and 6?, and 
neglecting terms of order higher than the second; using equations (13) and (24), 


I+ +a?) —/4(1+@3)—"/4 cos { (tan7! a,+tan—! a,)/2} 


G/R=—___—_.— Se 

: 1+( a?) ~*/4(1 +43) -4 [cos (tan-! a,+tan—! a,)/2—aq sin { (tan! a;+tan~! a,)/2}| 
1+ (tai) "4(1 +a3)~/4 cos {(tan™ a,+tan™! a,)/2} 

de/Ri= - oa 

° 1+ (1+aj)~ 1 +-a3) 4[cos (tan-! a,+tan~! a@)/2—a, sin { (tan! a,+ tan~! a2)/2}] 


(28) 


To help in understanding the significance of these results, a number of special cases are 
examined below. 
(v) Near the screen, a, and a,—0, and d,—>D,, d,—>D, for E,(x,y), as found”previously under (i). 
For the quadrature component £;(z,y) however, it is easily shown that 


K*+2K?/3-+1 ) 
=D ROK TS 


“41 QI2/ 1/2 
, d,<t { SRO at 


(29) 


Kf 2C+1 


where K=D,/D). 

(vi) Isotropic screen D,\=D,=D, K=1. 

Near the screen, d=D for E,(z,y) as found previously. 

From (29) d=D/y3 for E,(2,y). 

This result implies that the structure size of the shallow phase variations near an amplitude 
screen (or of the amplitude variations near the phase screen) is 1//3 times the structure of 
the screen itself. In addition, the correlogram of F,(z,y) is far from Gaussian; it can easily 
be shown by expanding (23) in powers of @ that, near the screen, 


p2(a,8)=exp{ —(a?-+f?) /2D"}-{1—(a2-+6?)/D?+ (a?+ 6°)2/8D*} (30) 


exactly. This is plotted in figure 3, and has an oscillating form. At any distance from an 
isotropic screen, (28) reduces to for F,(z,y): 


Oana) 
d= D{' a ST 


a 1+<a? 1/2 
i=D{ stay 


These relations are shown plotted in figure 4. The structure size of F,(z,y) actually has a 
maximum at intermediate values of z 





for E,(x,y): 
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Figure 3. Correlogram of the shallow phase variations 
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(vii) One-dimensional screen D,=D; D,->@, a0, K->0. Near the screen, 
d=D for E,(z,y), d=D/v5 for E,(2,y), 

showing the effect described in (vii) to an even greater extent. At any distance, 

for E,(z,y): 








i=D{ 1+(1+<a?)-'/* cos{ (tan~!a)/2} - 
ie 1+ (1+a’)~** [cos{ (tan~‘a)/2}—a sin{ (tan~'a)/2}] 
for E,(z,y): 
a=v{ 1—(1+<a?)~—/* cos { (tan~!a)/2} a 
1—(1+<a?)~** [cos{ (tan~‘a)/2}—a sin { (tan~'a)/2}] 


In the general case of an anisotropic screen, (28) can be used to find the structure sizes d, 
and d, in the z and y directions. To clarify the kind of behavior that may occur, figure 5 has 
been prepared. This shows how the contours of the correlation ellipse change shape with 
increasing distance from a random amplitude screen. The two screens shown are 

(a) isotropic, with D,=D,=1 

(b) anisotropic, with D,=2, D,=1. 

The distances from the screen were chosen as 0, », and rD?/x. It is interesting to note that 
the phase correlation ellipse for the anisotropic screen does not keep the same shape as the 
screen is approached—its minor axis is reduced by a greater factor than its major axis. This 
is because as it approaches the limiting case of an anisotropic screen—namely, a one-dimen- 
sional screen—the major axis of the phase correlation ellipse near the screen approaches the 
true structure size D, of the screen, while the minor axis is reduced by the factor 5 established 
above in (vil). 

The reason for dealing with this topic in such detail is to demonstrate that the structure 
size of the irregularities in the diffraction pattern produced by a shallow phase or amplitude 
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modulated screen is not always the same as that of the screen itself. The preceding lengthy 
mathematical argument is necessary to establish the precise effect of the diffraction process, so 
that existing experimental data may be interpreted with confidence. 


4. Cross-Correlation Relationships for E, and E, 


A good deal of attention has been paid by various workers (for instance, Jones, Millman, 
and Nertney [1953]) to the question of whether the signal in the diffraction pattern is the same 
as the signal at a random screen directly opposite to it. It must be emphasized that the 
question is not whether the signals are statistically similar, i.e., have the same structure size or 
fading time, but whether the detailed changes in signal are similar at the two planes. Expressed 
in more exact language, the cross-correlation is to be evaluated between, for instance, the 
amplitude of the signal at the screen and the amplitude at a distant plane. 

Another cross-correlation coefficient of considerable importance is that between the 
amplitude and the phase of the diffracted signal [Bowhill, 1957]. Jones, Millman, and Nertney 
[1953] assumed that this correlation arises as a result only of correlation between the amplitude 
and phase of the signal at the screen. It is shown here, however, that they may be correlated 
even for a screen which changes the phase only or the amplitude only. 

The theory is developed for a shallow amplitude modulating screen, and is extended to the 
case of a phase screen. The quantity /\(z,y) therefore represents the amplitude of the signal in 
the diffraction pattern (with the specular component removed) and F,(z,y) represents the phase 
of the signal. The quantities 








A(a,y)E\(z,y) and F,(z,y) FE, (2,y) (31) 





are to be calculated; the known values of £3(z,y), £3(2,y) and A?(z,y) then give the cross- 
correlation coefficients p,, between A(z,y) and E,(z,y), and p,. between F,(z,y) and £2(z,y). 
The quantities (31) are computed from the relations 





A(a,y) E, (x,y) = (1/2) E}(a,y)+ (1/2) A2(2,y) — (1/2)[A@,y) -Ai@y) P (39) 








(x,y) Ho(a,y) = (1/2) #7(a,y) + (1/2) A*(a,y) — (1/2) [Er ay) — Fa(a,y)J’. 
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The last terms in these equations represent the areas under the power spectra of 
[A(z,y)—E\(z,y)] and [F,(z,y)—E,(z,y)], since the mean square value of any random variable 
with zero mean is given by the integral of its power spectrum with respect to frequency. 

The power spectrum of [A(z,y)—,(z,y)] can be found from (11); 


P,(v,X) = Po(v,x) - cos { r2z(v? ++x?)/d}. (11) 
The spectrum of [A(z,y) —E,(z,y)] is therefore 
Po(v,X) — Pi (v,x) = Po(v,x)[1—cos { wz(v? +Xx?)/d}] 
and its power spectrum is given, exactly as in (6), by 
W(v,x)[1—cos { r2(v?+x?)/A} ]P?=4 [Wo(v,X)]-2—[We(v,x)]. 
and substituting in (32), 


A(z,y)E,(2,y) = (1/2) [E?(@,y)].+ (1/2) —2 [F3@,y) 0+ L/2EF@,y)].=2 EI@yYIe2—-1 





a relation which holds for any autocorrelation function at the screen. Using /7(z,y) from (24), 
the correlation p,, is given by 





= ee aT _y2- (1+a?/4)-/4(1+a3/4)-"* cos {[tan~! (a,/2)+ tan! (@,/2)]/2} 
pu=A (zy) Fi (a,y)/Ei(ay) {1+ (L+a?)~"/4(1+a3)-* cos {[(tan~! a+ tan! ay)/2]} }1/? 





for a screen with a Gaussian correlogram. 
For an isotropic screen, a;,=d.=a, and this reduces to 


f/( 2 2/9) VU: a 
pu= { (1+a?)/(1+?/2) }”?2(1+@?/4)-1. 
This function is plotted as the full curve in figure 6. The correlation is perfect near the screen, 
since then E,(z,y)=A(z,y). It decreases progressively to zero with increasing distance, 
reaching a value of 0.1 when a=7.2. 


The corresponding function for a one-dimensional screen, a;=d, d2=0, is given by 


oe 2(1+-a*/4) —"’* cos {[tan7! (a/2)]/2} 
Te + a?)-/* cos { (tan! a;t+tan! ay) } }1/? 





This function is plotted as the broken curve on figure 6. The correlation decreases much less 
rapidly with increasing distance from the screen than for the two-dimensional case. 
To find the correlation p,2., (32) is rearranged to give 








ky (x,y) E; (x,y) = 1/2— (1/2) LE, (x,y) — E(2,y) |? 





Figure 6. Variation of the correlation of the diffrac- 
tion pattern amplitude with the amplitude at the 
screen as a function of distance for one- and two- 
dimensional screens. 
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and the power spectrum of [/;(2,y) —2(z,y)] is given from (11) by 





[Pi(v,x) —P2(v,x) | [Pf (v,x) —P3(v,x) ] 


=W,(v,x)[cos {r2(v?-+x?)/X} —sin {r2(v?+Xx?)/} ]?=Wo(v,x) —Wo(v,x) sin { 2x 2(v?+x2)/d}. 
? ? J 


N 


O 





Ea) Ex) = (1/2028) | | W,(v,x) sin {2r2(v?-Ex2)/A} dnd 


= (eDDy2in) | [ [exp {2riz(v?+x?)/A} — exp { —2riz(v?-+x?)/d}] 


-exp {—2x?(Div?+ Dix’) /d?} dvdx. 


=(Dib,2in) { | exp { 2riv?(z+miD?/d)/d jaf exp { 2rix?(z+7tD3/d)/v} dx 


a 
eo 


exp {271x?(— s+ riDyp)/njax \ 


- | exp {Qxie"(— 2-+-eiDn) Aldo { 


and on integrating this expression by the usual methods 





FE, (2,y) E2(a2,y) = (1/2) +a) ~/4(1 +9) ~/*4 sin { (tan7!a,+tan7! az)/2}- 


Also, from (24) 





4. E?(2,y) - E3(2,y)=1—(1+4})-”7(1+43)-? cos? { (tan! a, + tan a) /2} 
and the correlation p,. between F,(z,y) and E,(z,y) is given by 


= (1+a3)-"/4(1+-a3)-"* sin {(tan7! a,+tan=*@,)/2} (38) 
pr [1 — (1-Fa3)— 721 Fa) -? cos? { (tan-!a;+ tan! a,)/2}]! i 





For an isotropic screen, @;=d.—a; this reduces to 
pu=(1+a")-? 

which is plotted in figure 7. The correlation between F,(z,y) and F,(2,y) has a maximum 

value of 1//2 at the screen, and decreases to zero with increasing distance. For an aniso- 

tropic screen, (33) must be expanded in powers of a; and a2. The correlation as a; and a,—0 
is given by 

a+ 1+k? 
1s "Tad P mc Q]1/3 
p12 (3a3-+ 2aa,+ 303]? [3K4+2K2+3]"? 





(34) 


where K=D,/D, expresses the degree of anisotropy of the screen, as in (29). Whatever the 
value of K, (34) shows that this maximum value of the correlation pz lies between the values 


pie=1/V2 for an isotropic screen a,=a, 





0.8 
0.6 
q 0.4 Figure 7. Variation of amplitude-phase cross-cor- 
relation with distance from a_ two-dimensional 
0.2 amplitude screen. 
0 2 3 + 5 6 8 





287 











and 


pi2=1/ 3 for a one-dimensional screen a,;> a». (35) 


The cross-correlation coefficients p,, and p;. between the phase components F,(z,y) and 
E,(z,y) must now be related to the amplitude and phase characteristics of the diffraction 
pattern. For an amplitude screen, it is apparent from figure 1 that for a shallow screen, 


Amplitude in pattern=+(specular signal) + £,(z,y) 
Phase in pattern — F,,(z,y)/(specular signal) 


since the positive phase direction is opposite to the positive , direction, taking an increase 
in phase to mean a retarding of the wave. This means that: 

Correlation between amplitude in pattern and at screen=p,,, and correlation between 
amplitude and phase in pattern=—p,2. Similarly, for a phase modulating screen 


Amplitude in pattern=+(specular signal) + F,(z,y) 
Phase in pattern (z,y)/(specular signal) 
and also 
Correlation between phase in pattern and at screen=p, 
Correlation between phase and amplitude in pattern= py». 


It has therefore been established that the correlation between the amplitude and phase in a 
diffracted signal is zero at large distances, but at smaller distances is positive for a phase modu- 
lating screen and negative for an amplitude screen. This occurs at a distance from the screen 
given by a~l, or a~D*/\. This is closely related to the Rayleigh distance for an optical 
system of aperture equal to the structure size of the screen irregularities. 

This result has been derived by setting up the wave solution. It may be partially con- 
firmed by a physical-optical argument, as follows. Consider a one-dimensional phase-modu- 
lating screen, with a phase profile ¢(z). The signal at a point near the screen is determined 
only by the properties of the screen in its immediate neighborhood. In fact, any small portion 
of the screen centered about 2» acts as a positive lens of focal length 


2 
F=—2s) E . (36) 


dx? L=2p 


The plane wave incident on the screen can be imagined as coming from a point source of 
waves, very far away. The image formed by the effective “lens” of the phase screen appears 
to have angular magnification =(1—2/F)~' where z is the distance of the point of observation 
from the screen. The image of the source appears the same brightness as the source itself, so 
the signal power arriving at z is proportional to the angular magnification. So the signal am- 
plitude is proportional to (l—z/F’)~“?; if z is assumed much less than F, approximately 


signal amplitude « 1+ 2/2F=1—(hz/4z) [Ss 
, r= 


on substitution from (36). Any correlation between the amplitude and phase in the pattern 
must therefore arise through correlation between ¢(z) and (—d*¢/dz’). 

In fact, it can easily be shown that if ¢(z) is a Gaussianly distributed random function, 
there is a correlation between these two quantities of [p’’(0)/p’’ ’’(0)]'/2, where p(a) is the cor- 
relogram of ¢(z). For a Gaussian correlogram, which has been assumed throughout this 
section, this reduces to 1/3; a result which is identical with that given for the same case 
in (35). It should be emphasized, however, that geometric optics can only be used near the 
screen, where the diffraction effects have not yet developed appreciably. 
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5. Effect of a Point Source of Waves 


In all the previous sections it has been assumed that there is no systematic variation of 
phase across the random screen. However, if the screen is illuminated with radiation from a 
point source a distance z’ from the screen, instead of from a plane wave source, the phase of 
the wave before passing through the screen varies with distance. The signal emerging from 
the screen is now given by 


Ey (0,Yo) =A(to,Yo) Xp { 2mi(z’?+- a+ ys) /?/A}. 


If the cone angle of the diffracted waves is small, the important regions of the screen are for 
Lo, Yo<2’. So the approximation 


(2/2 + 22+ yg) 1/22’ + 23/22’ + yf/22’ (37) 
may be used. 

Exactly as in section 3, the vector signal E(z,y,z) at a distance z from the screen can be 
found and the phase components £,(z,y,z) and E;(z,y,z) evaluated by a relationship similar 
to (9). In that equation, the phase of the specular component is independent of x and y since 
the specular component is a plane wave. In the present case, however, the specular wave 
is spherical and so the analogous equation to (9) is 


E(2,y,2) =(Ei(2,y,2) + iEa(2,y,2)] exp (2nil(e+2')*+2°+y"})} 


and for r,y<z, 
[(ete!) + a?t-y?}ltteet 2! +2°/2(242')+y2/2(2+2"). (38) 


The corresponding equation to (10) can then be written down as 


E(x,y,2z) exp {—2ni (ete) tat yn = (un) | {ff A(2o,Yo) « exp {2mi [(x—a)8 


+ (y—yo) rt (2’?+ a8+- yg) 2+ 2(1—s?/2—1?/2) |/X} dapdygdsdr 
or 


E(2,y,2,)= amy { {| A(29,yo) exp { 2riG/A} dadydsdr 


where 
G= (x—40)8 + (y—yo)r + (2? +2h + y6) /?#+2(1 —8?/2—1?/2)—[(e+ 2’)? + 2?+y?]}?. 
By substitution from (37) and (38), approximately 
G=(x—20) 8+ (y—Yo)r tea’ +.23/22’ + 8/22’ +-2—287/2 —er?/2—z—2' — 27/2 (et 2’) —y?/2(2+2’) 


= (4— 20) 8+ (y—Yo) +23 /22’ + y2/22’ —er?/2—zs?/2—2x?/2(2+2’) —y?/2(z+2’). 


The phase components £,(2,y,z) and E,(2,y,z) are given by 


Ey (2,y,2)= BALE (a,y,2 = cen) ff ff A(%,Yo) {exp (21G/A) + exp (—27G/d) } daydydsdr 


(39) 
E,(2,y,2)=S { E(ay,))}= (1/264 | | {J A(a»,yo) {exp (2iG/A) —exp (—2iG/d) }dapdyodsdr 
and the spatial spectra of F,(z,y,z) and F2(z,y,z) are given, as before, by 
P, ox={ E,(2,y,2) « exp {—2ri(va+xy)/A}dady 
‘ (40) 


P.o»= || E,(2,y,2) «exp { —2ri(va+xy)/rA}dady. 
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On substituting from (39) into (40) and integrating with respect to s, 7, z, and y, an analogous 

expression to (30) is obtained: 

cos {x(v?+X?) z(z+2’)/Az’} - {| A(ao,Yo) exp {—2ri (vay t+Xxyo) (2+2’)/dz2’ } dapdyo. 
(41) 


The differences between (41) and (11) can be accounted for by supposing that there are ap- 
parent values of z and \ which, when substituted in (11) and in all the results derived from it, 
just allow for the fact that there is a point source of waves at a distance z’ from the screen. 

In comparing (41) and (11), the double integral will first be considered. Bothy and x in the 
exponent are multiplied by the factor (z+ 2’)/z; so the angular spectrum must be scaled down 
in width by the reciprocal of this factor, relative to its width in the plane wave case. 

The ‘generalized autocorrelation function” of Booker, Ratcliffe, and Shinn (1950) given by 
the Fourier transform of the angular power spectrum, is invariant with distance in the plane 
wave case. With point source illumination, its width increases proportionally to (1+ 2/2’), 
and finally increases proportionally to z. In the case of vertical incidence of waves on the 
ionosphere z=<’, and the apparent structure size of the screen irregularities is just twice the 


, 





P,(v,xX)= a! 


actual structure size. 

These results can be appreciated physically by considering the behavior of the diffraction 
pattern as A—>0, and an analogy can be drawn with geometrical optics. If an amplitude screen 
(see fig. 8) is illuminated with a point source of light at a distance 2’ from it, the scale of the 
amplitude variations a distance z away from the screen will just be (1+ 2/2’) times the scale of 
the amplitude variations at the screen. 

The other respect in which (41) and (11) differ is in the argument of the cosine. The quan- 
tity (v?+2x?)2(2+2’)/2’ replaces (v?+-x?)z. But the quantities y and x have already been found 
to be multiplied by (1+ 2/2’); so z in (30) must be replaced by 


meter). 2 gs a 


Ta (1+ 2/2’ a7) ett 





where Zer: is the effective value of z. In the case of vertical incidence on the ionosphere, 


e=2’, and Zen=2/2. 
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Figure 8. Diagram showing the increase in structure 
size of the amplitude variations with increasing 
distance from a point source of waves. 
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It is interesting to note that 





~ at 
“eff @ 


which is exactly comparable with the object and image distance relations for a zone plate, the 
theory of which is closely related to that given above. 


6. Effect of Oblique Incidence 


If a plane wave is incident at an angle 7 to a diffracting screen, the equations corresponding 
to those of section 3 can be deduced by an exactly similar method. If y axis is in the plane 
of the incident wave, we may write s=sin 7+8,, where (sin~'s—sin“'s;,) is the scattering angle. 
The phase of the specular signal at a distance z from the screen is now 

(2miz cos 12)/d 
and (9) must therefore be written 


E(2,y,2)=(Ei(2,y) + 71E2(2,y)] exp [(2riz cos 7)/d]. (42) 
The phase of the scattered wave at distance z is again given by 
Qwiz{ 1+ (s-2—1)7!+ (9? —1)“1}-¥2/a. 
Putting s=sin 7+, and expanding in powers of s; and 7, this becomes 
2riz(cos i—s, tan 1— } sj sec? i— 37? cos? 7)/d 


and the exponent in the resulting equation for (£,+7£;), analogous to (10), is 
; ae ae ee ee 
2ri4 (x—X—2 tan 2) + (y—Yo)r—5Fsi 2 sec® 1—sr*z cos® 7 /D- 


Comparing this with the exponent in (10), it is evident that 

(i) the origin of z has been effectively displaced a distance z tan 7 in the plane of propa- 
gation—namely, to the ‘mid-point of the path’’, in propagation terminology. 

(ii) the effective distance of the observing point from the screen has been multiplied by 
factors of sec* 7 and cos* 7 respectively for the x and y variations in the screen. 

To make this more easily comprehensible, consider the case of oblique incidence on a 
diffracting screen with an isotropic spatial correlation function, Gaussian in form. The 
structure size of the interference pattern is found by putting 


Za 
a,;=a sec? 7 : 
} k=sec’ i 


d,=a cos? 7 


in (28). Near the screen, from (29), 


a—D{ i+ (2/3) sec!? i+1 , 


sec*4 1+2 sec? 1+5 


i.-D{ sec?4 7+ (2/3) sec!? i+1 a 


5 sec*4 1+2 sec}? i+1 








Another important effect of oblique incidence, which does not enter explicitly into this 
analysis, is the following. If the random screen is composed of a “buckled specular reflector”, 
it will act as a phase screen for normally incident waves. However, at oblique incidence, the 
phase change it causes in the wave is scaled by a factor cos 7. This is related to the familiar 
phenomenon that an optical surface becomes more nearly a specular reflector at grazing inci- 
dence. The same effect is observed for radio waves. 
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7. Conclusion 


Detailed calculations have been given for the statistics of a radio wave diffracted by a 
random screen. Theresults, some of which have been used in a previous paper [Bowhill, 1957], 
now take account of anisotropy of the diffracting screen, and oblique incidence of the waves 
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The radio antenna is viewed as a space frequency filter with an output just equal to the 
convolution of its transfer function with the radio field illuminating its aperture. An equivalent 
uncertainty principle limits the accuracy with which the spatial distribution of the radio 
field may be determined. The radio field generated by an antenna is distorted in passing 
into the ionosphere, generating new space frequencies which is the information contained in 


the field at the receiving antenna. 


The energy diffracted into the different orders will appear 


to arrive from different directions, and the angle of arrival for a given order will be a func- 


tion of the radiofrequency. 
1. Introduction 


Radio fields are functions of space and time of 
the form 


flot+k-r), 


where w is 27 times the time frequency in cycles per 
unit time and 

k=tk,+jky+kk, (1) 
is 2r times the space frequency * in cycles per unit 
distance along the given direction. A monochro- 
matic plane wave propagating in a homogeneous 
isotropic medium will have space frequencies ranging 
from zero along directions in the plane of the wave 
front to a maximum equal to 27 divided by the wave 
length along the direction of propagation. For direc- 
tions making an angle 6 with the direction of propa- 
gation, the space frequency, or wave number becomes 


ko=— cos 6. (2) 


If the time frequency of the plane wave is not con- 
stant, then the space frequency will not be constant 
along directions other than #=90°, which cor- 
responds to zero space frequency. Thus, a radio 
wave modulated in the time domain has a correspond- 
ing modulation in the space domain. 

Even with a source monochromatic in the time 
domain, there is a spectrum of space frequencies 
which is generated by the radiating system, and 
this spectrum is a unique property of the current 
distribution associated with the radiating structure. 
In a homogeneous isotropic medium, the field every- 
where is described in terms of this spectrum. In 
nonhomogeneous media this spectrum is continually 
modified as the wave propagates along, the inhomo- 
geneities generate new space frequencies which feed 





‘Contribution from Smyth Research Associates, 3555 Aero Court, San 
Diego 11, Calif. 

2 Paper presented at Conference on Transmission Problems Related to High- 
Frequency Direction Finding, at UCLA, June 21-24, 1960. 

' The space frequencies are frequently called wave numbers and sometimes 
modes, but it is felt that the term space frequency is more appropriate in this 
discussion. 





on the energy of the original field. Other treat- 
ments of the effects of inhomogeneities, with special 
application to the ionosphere, have been reviewed 
by Ratcliffe [1956]. 


2. Space Signals Developed Over 
Ionospheric Links (1.5-15 Mc/s) 


Tropospheric effects at these frequencies will be 
macroscopic in nature and the gross decrease in 
refractive index with elevation will produce a down- 
ward bending which will amount to some 22 min at 
an elevation angle of 2°, and 3 min at 20°. The 
ionosphere on the other hand, will appear dense and 
turbulent to the radiation propagating in the forward 
direction. 

If, at each turning point of the radio ray, a reflect- 
ing plane parallel to the earth is assumed, and the 
intervening medium is considered homogeneous, 
then the radio field developed at the earth is just that 
associated with the image of the transmitter. Since 
the transmitter is at a large distance from the receiver, 
the waves arriving are almost plane and the phase 
wave length in a plane tangent to the earth’s surface 
is the ratio of the radio wave length to the cosine 
of the elevation angle. 

The medium between the reflecting plane and the 
earth isnonhomogeneous. In addition to the more or 
less regular change of electron density with elevation 
there are irregular patches of ionization and, in 
particular, filaments left in the wake of meteors at 
elevations between 80 and 120 km. In propagating 
through these irregularities, the radio waves are 
distorted and the image of the transmitting antenna 
in the ionospheric mirror becomes blurred. There 
are ionospheric conditions where the ray path may 
indicate directions up to 180° off the great circle 
path between the transmitter and the receiver. In 
any case, all that is available to the observer is the 
output voltage from the antenna which is generated 
by the radio field in its aperture. Certain deductions 
are made based upon the sampling of the radio field 
with the antenna which is nothing more than a band- 
pass space frequency filter. If the antenna is 
rotated, this is equivalent to shifting the pass-band 
across the space frequency domain in a given direc- 
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tion. If the radio field consists of a single plane 
wave, the output voltage from the antenna will be 
its transfer function for this specific cut in the space 
frequency domain. For an arbitrary radio field 
illumination, the output from an antenna is the 
convolution of the radio field with its transfer 
function. As the antenna rotates, there may be 
several peaks in the output, each corresponding to 
increased signal derived from certain of the space 
frequencies making up the space signal. These 
results are usually interpreted as multipath trans- 
mission, where the signal is assumed to be coming 
from several directions. 


2.1. Radio Field Distortion in 
Nonhomogeneous Medium 


Returning to the assumed case of a plane reflector 
at the ray turning point in the ionosphere, the radio 
field at the receiver can be considered as produced 
by transmission from the image transmitter through 
the inhomogeneous stratum, the ionosphere. Let 
it be assumed that the ionosphere is plane and 
horizontally stratified and the transmission problem 
is depicted in figure 1. If the separation of 7” 
the transmitter image and R, the receiver (along the 
2 direction) is large compared with the layer thick- 
ness 2/, the radiation incident on the layer may be 
assumed plane. In this case, the solution of the 
field as a function of 2 is all that is required. 

A plane wave polarized normal to the incident 
plane will produce a field in the inhomogeneous 
medium of the form 





E,(z,t)=A(z,t) sin [wt—k,z+¢(z,t)] (3) 
This field satisfies the inhomogeneous wave equation 
ok, yD =i (z.° 


If it is assumed that A and ¢ are slowly varying 
rag of z, so that to a first approximation f 


(z., DF os) may be neglected, then the unperturbed 
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Figure 1. Jonospheric link—downcoming waves generated by 
horizontally stratified inhomogeneous layer. 
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solution becomes 





E,=A sin @t—k,2z+¢), (5) 
and 
or = — Ak, cos (wt—ke+9). (6) 


Differentiating (3) with respect to z and substitut- 


z 


Oz 
A’ sin [wt—kz+¢]+Ag’ cos [wt—kz+¢]=0. (7) 


from (6) yields 





ing for 


Taking the derivative of (6) with respect to 2 and 
substituting into (4) yields 


ss : De he cos [wt—kz+g]—Ak oe sin [wt—kz+¢] 


=1(E. 3) ® 


Multiplying (7) by sin [ot—ke-+y] and (8) by (;) 
cos [wt—kz+¢] and 


adding gives 


=f (z., OE o) dhe. (9) 


In a similar manner 


oe —!1 
>= FA f (Ex, 9 Dt one) sin [wt—kz+g]. (10) 


It is assumed that the rates of change of the ampli- 
tude and phase vary little in a space period 2z/k. 


) Oe . : 
The mean values fo — and = in one period are 
~ 





given by: 
aoe 
= oo _ os )eoskzd(ke), (11) 
= , ok ) 
an _ I( B.S = sinkzd(kz). (12) 


It f(Z,, 2) i, is a function of E, only, the average 


change in the amplitude of the radio field over a wave 
length is zero. On the other hand, equation (12) does 
not vanish, and the phase change per wave length is 
the first order effect associated with the penetration 
of the field into the inhomogeneous medium. For the 
case where the dominant term in f is a function of 








oH, 
ot 

effect, as would be expected in cases where dissipation 
overrides linear and nonlinear restoring forces. 
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In the case of high frequency radio wave transmis- 
sion through the ionosphere, it is expected that the 
cases of interest are those for which the first order 
effect is a phase distortion. In any event, for direc- 
tion finding applications, conditions which produce 
attenuation will reduce the space signal at the receiv- 
ing antenna and limit the usefulness of this propa- 
gation mode. As the radio wave front progresses into 
the ionospheric layer, some parts of it are advanced 
in phase with respect to other parts as a result of the 
blobby, granular index of refraction structure of the 
medium. It will be assumed that the first order effect 
of the inhomogeneous layer on the emerging radio 
field is a phase “distortion ; which, of course, will pro- 
duce an amplitude variation in the far field. An 
optics system is a good example of this phenomenon: 
a plane light wave incident upon a lens is distorted in 
phase but not amplitude in a plane immediately 
behind the lens, and this aperture field distribution 
creates the usual variation in intensity in the focal 
plane. Thus, a distortion in the phase of the radio 
field by the ionosphere will create new space fre- 
quencies which modulate the space signal at the 
earth’s surface. The movements of the ionospheric 
irregularities will create a time and space variation 
of the phase distortion which in turn produces a time 
and space variation of the amplitude of the radio 

field at the surface of the earth. 

In order to describe the output of the receiving 
antenna, it is necessary to specify the space signal in 
its aperture plane. The phase distortion produced 
by the ionosphere may be described in terms of its 
Fourier components. For simplicity, let it be as- 
sumed that at z=/ the phase distortion has a domi- 
nant component along the y direction with ampli- 
tude, a, and space frequency or wavenumber k,. The 
amplitude of the radio field is assumed constant and 
the phase in equation (3) is given by 


$2=1(Y) =a sin kay. (13) 
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Figure 2. Energy diffracted into first three orders as a func- 
tion of the radio phase distortion. 





It can be shown [Smyth, 1959] that this phase dis- 
tortion over a distance y=2 6 will diffract energy 
into other radio wave space frequencies with ampli- 
tude given by 


ain — mk,)b 


— oS on 


A (ky) rm ye aia) — 


where J» (a) is the Bessel function of order m, and a, 
is the amplitude of the phase distortion.’ The radio 
space frequency in this case is its value along the 
y-axis, where, of course, the cos @ in equation (2) is 
replaced by sin 6. If the region over which the phase 
distortion takes place is large compared with the 
radio wave frequency, then the amplitude of the new 
space frequencies generated by diffraction will differ 
from zero only in the immediate vicinity of the 
spectral lines 

Rin wh, (15) 


since the function 


sin (k y— mkq)b b (16) 


a) 











is a maximum value for zero argument. The inten- 
sity of the radio field diffracted into each order is 
given by the square of the amplitude function, J», 
(a), which is just unity, when m=0 and a=0. This 
is the case of no distortion and for this case equation 
(14) reduces to a single line, k,=0, which is just the 
original plane wave. It is of interest to look at equa- 
tion (14) and observe the effect of radio frequency 
on the amplitude of the diffracted orders. Figure 
(2) gives J?, (a) for m=0, 1,2. It is clear that the 
energy diffracted out of the original plane wave never 
completely reappears in any one higher order. For 
small phase distortion we have 


J, (Fr), 


J; (F 1), (17) 





where a= (27L/\) is the amplitude of the radio phase 
variation from its mean value. Capital Z is the 
deviation of the radio phase from the mean. 


3. Antenna as Space Frequency Filter 


3.1. Antenna Response 


The output voltage of an antenna is determined 
by its filter characteristics and the electromagnetic 
field in which it is immersed. For a given orienta- 
tion of an antenna, the output voltage is a single- 
valued function ot the distribution of electromagnetic 
field over the effective aperture of the antenna, but 
unfortunately the converse is not true, there are an 
infinite number of field distributions that will yield 
the specified voltage output. 


295 











Let the aperture plane of the antenna be the zy 
plane z=0. If the antenna is fed radio energy at a 
constant time frequency, a resulting field distribu- 
tion E(z,y) will be maintained over the aperture 
plane. This aperture distribution may be repre- 
sented by a spectral distribution of space frequencies 
in much the same way that a time signal may be 
represented in the time frequency domain. Jn terms 
of plane wave functions, the space representation of 
the field is given in the space frequency domain bv 
the three dimensional Fourier transform: 


Alkeky k= f | | E(z,y, 2) 


ether tikwtti.? dadydz, (18) 


where A(k,,k,,k,) is the amplitude of the plane wave 
associated with the direction (k;,k,,k.), and, as usual, 
the space frequency components in the cartesian 
coordinate system are related as follows: 

B+B+R=-. (19) 
For the case where all the space frequency com- 
ponents are real, the fields are propagated waves 
and the amplitude function A(k) evaluated for the 
field distribution over any plane is just the polar 
diagram of this particular radio field distribution 
[Booker and Clemmow, 1950]. 

Fields in the zy plane, for which k?+k?>k? 
correspond to waves beyond cutoff in the z direction 
and these make up the so-called ‘induction field”’ of 
the aperture distribution. The distorted radio field 
emanating from the ionosphere may have components 
in the current distribution which are of this type, in 
which case the antenna will be in the induction field 
for these components. At the distance involved, 
only the spectral components of the induction field 
in the immediate region of cutoff will be felt at the 
ground, since the higher space frequencies will be 
attenuated rapidly with distance. 

Considering the field in the aperture plane z=0, 
the space frequency representation of the radio field 
and the field distribution in space are related by the 
two-dimensional Fourier transform pair 


A (ke, ky) ={ E (a, ye" *dady, 
E(x,y) =|| A (kz, ky) e~*"-*"dkedke,. (20) 


The output of an antenna placed in this field dis- 
tribution will simply be the superposition of the re- 
sponses to the various plane wave components in 
the field illuminating the antenna aperture. In 
terms of the space frequencies, the space signal is 
just 


A (kz, ky) = | | A (kz, ky) 6 (ka—kez ky — ky) dhezdky, (21) 
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where 6 is the impulse function in the space frequency 
domain; that is, a plane wave of unit strength arniv- 


ing from some direction k with the property 


5(k.—kzky—ky)=1 k= he 
b=k, 
=0 kk, 

ky Aky. (22) 


Let the impulse response of the antenna be G(k,,k,), 
then the output will be the sum of the responses 
created by the input signal and may be written thus: 


V(kek,) = | { G(kz—ki key —hej) A (ki ep dkeidk;, 
or 
k/)dkidke; 


(23) 


Vieoky)=| | Ghz ky) A(ke—ki ky 


since the output must be the same linear combina- 
tion of the displaced signal as the displaced response. 
Now, if the signal is the unit impulse function, the 
output will be just the space frequency filter response 
function. 


Viki ky) = [ [GG 4;)8(.— ke yk dks =G (key) 
(24) 


Since the aperture field input may be expressed as 
the superposition of impulse functions with appro- 
priate amplitude and phase, the output will be the 
convolution of the input space frequencies and the 
impulse response of the antenna, which is the voltage 
“oain-function” of the antenna. Thus, the antenna 
output is uniquely given by the input signal; on the 
other hand, the space signal is not a single-valued 
function of the output. 

It is clear that the output given in equation (24) 
will be the filter response function when the input 
is a plane wave function. When the input is of 
more complex nature, the output will not be simply 
related to the filter transfer function. For some un- 
known reason, many radio people have chosen to 
assume that the input signal is always of the form 
of a single plane wave, and that the output resulting 
from the convolution of the actual input signal and 
the filter response function of the antenna is a new 
response function for the antenna. This practice 
of absorbing the unknown aperture distribution of 
the radio field into the antenna filter function has 
led to considerable confusion, where some people 
have gone so far as to state that the antenna placed 
in a nonplane wave radio field “lost its gain.” 4 A 
more sophisticated way of putting the same thing 
is the more current statement that there has de- 
veloped an ‘“antenna-to-medium coupling loss.” 
Of course, the antenna is excited by the radio field 

4 This reminds us of the small boy riding his bicycle down a rough road, con- 


cluding that each time a wheel rolled over a depression that the wheel lost its 
diameter. 








nCy 
riv- 


Ses 
US? 


23) 


na- 
ise. 
the 
nse 


ky). 
24) 


as 
ro- 
the 
the 
age 
ina 
the 
ued 


24) 
put 
) OF 
ply 
un- 
to 
orm 
‘ing 
and 
1ew 
tice 
. of 
has 
yple 
ced 
A 
ing 
de- 
ss.” 


ield 


, cone 
ost its 





in its aperture and this radio field is influenced by | 


the medium; but in electromagnetic theory the 
differences in physical properties of media are ac- 
counted for by charge and current distributions. 
Thus, the only way an antenna can be coupled to a 
medium is through the fields. For example, it 
would be of questionable value to couple to a wave- 
guide capable of propagating many modes without 
regard to the mode being transmitted at the time. 
Even though the so-called medium remains the same 
independent of the mode being propagated, the 
fields are quite different, and the concept of antenna- 
to-medium coupling loss is of no real utility. 

In general, the field distribution across the aperture 
of the antenna will change with time and conse- 
quently the output from the antenna most likely 
will change with time. Both tropospheric and 
ionospheric transmission studies indicate that the 
functional form of the aperture field, E(x, y,), 
changes with time associated with a continual re- 
distribution of radio energy among the different space 
frequencies. In direction finding experiments this 
corresponds to a meandering and spread in the ap- 
parent direction and location of the transmitter. 
Even when the propagation effects are neglected, a 
time modulation of the radio field will produce 
an uncertainty in the determination of the angle of 
arrival. This effect has been discussed by Dausin 
et al., [1959], who have shown that the error can 
become appreciable when the time represented by 
the reciprocal of the bandwidth is less than the 
propagation time across the antenna aperture. 


3.2. Uncertainty Principle 


If the space signal could be sampled at a point by 
an isotropic antenna, and even if the radio field at 
the point could be measured with unlimited preci- 
sion, the signal variation in space or direction of 
arrival would be completely undetermined. To 
gain information on the space signal which contains 
information on the arrival direction requires the 
utilization of space. Again considering the zy 
plane z=0 as the antenna aperture plane, the range 
of uncertainty in the space frequencies making up 
the signal along a given direction in the plane is 
related to the available aperture in that direction 
by the relationship 


A\k|-5|r|> 6 (25) 





This is, in essence, the celebrated ‘Uncertainty 
Principle.”” The direction property of an antenna 
is a good example of this principle; if it is desired to 
radiate most of the energy in a specific direction then 
the effective aperture must be large. At first glance 
one might conclude that supergaining violates this 
relationship, but this is not the case as can be seen 
from the fact that the effective aperture of the 
antenna increases with increase in antenna gain 
according to the relationship (25). Diffraction of 
light by a slit clearly demonstrates this phenomenon. 
Any attempt to pass a plane wave through a confined 
region generates new space frequencies which spread 
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SPACE SIGNAL SPACE FREQUENCY SPECTRUM 
Constant space signal of finile extent and its space 
Jrequency spectrum. 


FIGURE 3. 


the energy out about the original direction & in 
inverse relationship to the width of the slit. 

Suppose that a plane wave is incident normally on 
the zy plane and its amplitude is measured over a 
strip of finite width ér. If it is assumed that the 
space signal is zero outside this range, then the space 
frequency spectrum is given by 
sin rk(dr) 


C(k) = (6r) —FEGr) (26) 





This space signal and its spectrum is shown in 
figure 3. On the other hand, any continuous space 
signal with continuous first derivatives will appear 
to be constant over a space interval chosen small 
enough. If, in this case, the space signal is assumed 
to be constant throughout space, the space frequency 
spectrum will reduce to a single discrete line. Clearly, 
then, there are an infinite number of spectral distri- 
butions which would yield the same space signal 
over a given range. To define a spectrum uniquely, 
the signal must be defined over the entire xy plane. 
Thus, given a region of space over which the signal 
will be sampled, (25) places a limit on the accuracy 
with which the space frequency may be specified. 

An example of some practical interest is the space 
signal distortion generated by the transmission of a 
plane wave through an inhomogeneous medium such 
as the earth’s atmosphere, including both the tropo- 
sphere and the ionosphere. This will deteriorate 
further the accuracy with which the space frequency 
may be specified and the relationship in (25) is indeed 
an inequality. The question asked concerns the 
possibility of filtering out the space frequency noise 
as a means of reducing the magnitude of the product 
on the left side of (25). If the characteristics of the 
space frequency noise are known, it should be pos- 
sible to design optimum space frequency filters to 
reduce the uncertainty in the direction of arrival of 
the desired signal. 
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The distribution function of envelope voltage for short samples of atmospheric radio 
noise as received by a communications receiver in the VLF range always shows a marked 
departure from that obtained for Gaussian noise. If it is considered that this departure is 
caused by strong noise pulses which do not overlap in time, the effect of changes in the 
receiver bandwidth on the observed distribution function can be deduced by consideration 
of the changes in the receiver impulse response. A transformation can be obtained which 
gives an excellent approximation to the change in a mathematical representation of the 
distribution function in the range of probabilities below 1 percent. Empirical relationships 
are suggested which give useful estimates of the change in the distribution function over 





the total range of probabilities. 


1. Introduction 


In the VLF range, atmospheric radio noise is one 
of the important factors in a system design. If all 
of the manmade interference is controlled by alloca- 
tion procedures and good engineering, the atmos- 
pheric noise provides the ultimate background dis- 
turbance from which the desired signal must be 
separated [CCIR Rpt. 65, ITU, 1957; Watt, Coon, 
Maxwell, and Plush, 1958]. Studies of the noise 
must be of a statistical nature, but there are problems 
which differ, for instance, from those of statistical 
studies of thermal noise, because the atmospheric 
noise is a nonstationary process. Since the statistics 
of the process change with time, the accuracy of 
measuring statistical parameters cannot be in- 
definitely increased by increasing the length of time 
of the measurements. Experience has shown that 
when measuring the amplitude probability distri- 
bution of the noise envelope, samples of noise 10 to 
20 min in length are short enough to avoid difficulties 
due to the nonstationary characteristics of the 
process, but are long enough to give useful informa- 
tion for system design purposes [Hoff and Johnson, 
1952; Watt and Maxwell, 1957]. 

One of the problems which occurs in system design 
is determining the percentage of time that the 
envelope of the noise will exceed the signal level, and 
how this varies as the receiver bandwidth is changed. 
In principle this can be calculated precisely if 
enough is known about the statistics of the noise. 
An applicable procedure, as described by Widrow 
[1957], is to consider a sampled version of the input 
noise for which the joint probability distribution of 
the samples can be determined, and to use this to 
calculate the probability distribution of a sampled 





‘Contribution from Central Radio Propagation Laboratory, National Bureau 
of Standards, Boulder, Colo. 

? An earlier version of this paper was presented at the Symposium on VLF 
Radio Waves held in Boulder, Colo., January 1957. 

‘Present address: Lockheed Aircraft Corporation, Palo Alto, Calif. 





version of the filter output. If the sampling rate at 
the input is high enough, the probability distribution 
of the output samples will accurately represent the 
probability distiibution of the continuous output. 
The difficulty with this procedure is that the interval 
between samples at the input must in general be 
short relative to the correlation time of the noise 
[Ragazzini and Franklin, 1958], so that a very high 
order joint probability density function must be 
calculated for the filter output. The purpose of this 
paper is to present a technique which requires much 
less computation, but gives an accuracy commen- 
surate with the statistical knowledge of the noise. 


2. Noise Characteristics 


Figure 1 shows two samples of measurements of 
atmospheric noise envelope distributions, one of rela- 
tively high dynamic range conditions and the other 
of noise of relatively low dynamic range; the dotted 





0.99 v 





P=Probability of Exceeding v ‘+ 
o—-o— High Dynamic Range Sample | 
lati Low Dynamic Range Sample | 
—|------ Distribution for Thermal Noisef 


0.98 ;}— 


0.85 -— 














100K 


a 
Ss 


v 
RELATIVE VOLTAGE 


Frieure 1. Typical measurements of atmospheric noise envelope 
distribution. 
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line represents the measurement which is obtained 
for thermal noise. At the low voltages which have 
a high probability of being exceeded, the atmospheric 
noise measurements parallel those for thermal noise ; 
this is characteristic of a phenomenon consisting of 
a large number of overlapping pulses, no one of which 
contributes a significant proportion of the total 
energy. At the high voltages, however, which have 
a low probability of being exceeded, the measure- 
ments depart radically from the shape of the thermal 
noise curve, in the form of a much higher probability 
of exceeding the voltage levels in the range of prob- 
abilities of 0.3 and below. Oscillograms of atmos- 
pheric noise show that the high voltage levels are 
attained by occasional noise pulses which greatly 
exceed the general level of the noise [Yuhara, Ishida, 
and Higashimura, 1956] and that the pulse shape is 
that of the impulse response of the receiver. These 
strong pulses occur infrequently and two such pulses 
very rarely overlap. 


3. Probabilities Associated With Impulses 


The effect of bandwidth changes on the probabili- 
ties associated with noise impulses can be calculated 
by considering the effect of bandwidth changes on 
the usual impulse response of the filter. The effect 
of increasing the bandwidth is to cause the impulse 
response to become shorter in time, and higher in 
amplitude [Guillemin, 1956]. More specifically, if 
the network characteristics are changed in such a 
way that the shape of the impulse response is pre- 
served while the bandwidth is multiplied by a factor 
B, the voltage scale of the impulse response will be 
multiplied by a factor B and the time scale by a factor 
1/B, as illustrated in figure 2 for B=2. This rela- 
tionship also holds for the envelope of the response 
of a band-pass circuit whenever the circuit possesses 
a low pass analog, as is normally the case with 
receiver circuits [Aigrain, Teare, and Williams, 1949; 

iuillemin, 1953]. Further, the probabilities asso- 
ciated with the impulse response are determined by 
the time scale; for example, the probability that a 
voltage level of one will be exceeded by the narrow 
band response is determined by the ratio of the time 
7t to the total observation time. Because of the 
relationship between changes in the voltage and time 
scales, the doubled bandwidth impulse response 
exceeds a voltage level of two for a time of 7/2, 
which represents precisely one-half of the probability 
corresponding to a time r. 

Thus the amplitude distribution of a short sample 
of noise which contained only one impulse would bk» 
determined by the time and voltage scales of that 
impulse response; if a particular voltage v is measured 
to have a probability P of being exceeded when a 
single impulse is observed through a network of band- 
width Aw;, and simultaneously the impulse is ob- 
served through a network of bandwidth Aw:= BAa,, 
the voltage Bv will be measured to have a probability 
(1/B) P of being exceeded. This relationship also 
holds for a noise sample which consists of a number 
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FigurRE 2. Envelope of band-pass impulse response for two 
bandwidths. 


of impulses, provided that the response of the net- 
work to one impulse does not overlap the response 
to another. 

These requirements are never exactly met in prac- 
tice because of the exponential decay of physical 
networks and because of the presence of thoroughly 
overlapping low level noise. However, the high 
voltage end of the distribution of VLF atmospheric 
noise appears to represent an excellent approximation 
to the requirements. Figure 3 shows the high volt- 
age portion of two representative sets of measure- 
ments at two bandwidths and also shows the results 
of translating the narrow-band measurements by the 
ratio of bandwidths to expected measurements at 
the wider bandwidth. It is evident that for both the 
low and high dynamic range noise conditions the 
points calculated from the narrow-band measure- 
ments agree with the measurements at the wider 
bandwidth for probabilities less than 0.01. 


4. Mathematical Representation of the Noise 


In dealing with a non-Gaussian random variable, 
it has proved profitable for some purposes to consider 
it as being generated as the output of a nonlinear 
resistance network which has a Gaussian input, even 
though the actual phenomenon may be very much 
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Ficure 3. Atmospheric noise distributions showing effect of translating measurements made at a narrow band- 
width to expected values at a wider bandwidth. 


different from this [Smith, 1959a and bj]. In the 
situation being considered here where the envelope 
of the noise is the quantity of interest, we can 
represent the envelope of the actual noise as being 
obtained from an envelope ‘(etector having a Gaus- 
sian input, followed by a ..colinear resistance net- 
work. The justification for “>ing this is simply that 
the resulting probability di-iribution can represent 
the probability distribution of the actual phenomenon 
quite accurately, and the model proves to be a con- 
venient one for the type of mathematical manipula- 
tions which are required. 

A nonlinear network which can produce the correct 
probability distribution has a voltage transfer char- 
acteristic which is a three term polynomial of the 
form 
b+1 

2 


+azy’ (1) 


where 4d, 2, a3, and 6b are chosen to match the 
atmospheric noise distribution under consideration; 
the probability that a voltage level v is exceeded at 
the output of the network is given by 


P=e-*, (2) 


The variable y may conveniently be considered just 
as a parametric variable, although it is the envelope 
voltage of a Gaussian variable having a variance of 
one-half. 


v=hhy +ay 





Each term of the polynomial, if considered sepa- 
rately, plots as a straight line of log-log of proba- 
bility versus logarithm of voltage, which are the 
coordinate scales used on the graphs of distributions 
in figures 1 and 3. The first term represents a dis- 
tribution having complete overlapping, such as 
thermal noise; the second and third terms represent 
departures from this distribution. The third term 
of the polynomial, a3y’, is the dominant term at the 
high voltage, low probability end of the curve, when 
the phenomena forming the distribution appear to 
be strong, nonoverlapping pulses. 


5. The Effect of Bandwidth Changes on the 
Mathematical Model 


Since the third term of the polynomial is consid- 
ered to represent nonoverlapping impulses, the effect 
of bandwidth changes can be obtained by the proce- 
dure described in section 3. That is, for an increase 
in bandwidth by a factor B, any point on the original 
line at voltage v and probability P is transformed 
to a voltage Bv at a probability P/B. After this is 
done, it is found that the points no longer lie on a 
straight line, but along a line having a slight curva- 
ture, as illustrated in figure 4 for an exponent of 8 
and a bandwidth increase of ten times. The curva- 
ture introduced in this way is so slight, however, 
that the curve may be replaced by a line tangent to 
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it at the point where P=e~*; as is apparent from 
figure 5, this is an excellent approximation. 

It is possible to derive a general formula for the 
effects of bandwidth changes on the exponent and 
the coefficient of the third term of the polynomial. 
For a bandwidth increase by a factor of B, the expo- 
nent 6, for the wider bandwidth is obtained from 
the exponent b, for the narrow bandwidth by the 
relationship 


8 


bu=3-T B 


b, (3) 


where L,B is the natural logarithm of B. The co- 


efficient is obtained from 


(<z3) bn / 
aae=009B L(8—L,B)8 re f (4) 
where dz, and dz, are the coefficients for the wide and 
narrow bandwidths respectively. The derivation of 
these relationships is given in the appendix. 

For the other terms of the probability function the 
effects of bandwidth changes are not as apparent as 
for the third term which describes the behavior at 
high voltages. The second term of the polynomial 
is important in the moderate voltage range where 
some overlapping of noise pulses occurs. For the 
values of exponents and coefficients which are re- 
quired to match the mathematical model to the 
observed atmospheric noise conditions, it happens 
that the second moment of the distribution is almost 
entirely determined by the second and third terms 
of the polynomial. Since for normal communica- 
tions bandwidths the total noise power received is 
directly proportional to the bandwidth, it is possible 
to adjust the second term to obtain the correct 
behavior for the second moment of the distribution. 





The procedure is to calculate the changes in the a, 
coefficient and exponent 6 as described above, and 
to insert these into the expression for the mean 
squared voltage of this form of distribution, which is 


(v*)=a?+a3r (*)+ar+1) +2007 on 


+2a,a,T (75?) +2000 (7) (5) 


where (v’) is the mean squared voltage. 

For a bandwidth increase by a factor of B, the mean 
squared voltage for the wide bandwidth is B times 
the mean squared voltage at the narrow bandwidth. 
With this and a value for a, equation (5) can then be 
solved for the coefficient dy. 

The first term of the polynomial, involving a, 
represents a situation where complete overlapping 
occurs. The slope of a distribution represented by 
only this term does not change with bandwidth, but 
only its amplitude, as long as the complete over- 
lapping is maintained. If this were the only term 
in the distribution, its amplitude would change 
precisely as the square root of the bandwidth. 
Since there are other terms present, however, there 
will be an interchange of energy between the various 
portions of the distribution, and the net effect will 
be that the amplitude of the first term will change 
more slowly than this; an increase according to the 
four-tenths power of the bandwidth appears to fit 
the available experimental data quite well. 

These procedures work quite well with noise 
samples of low and moderate dynamic range, but 
an error appears when used with noise of high dy- 
namic range, that is, values of 6 of about 11 or 
greater. This occurs because the mathematical 
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model being used is not sufficiently complex. With 
the very high dynamic range conditions, a sub- 
stantial contribution to the integral for the mean 
squared voltage occurs at probabilities below 10~°; 
in this region the tangent line and the curve repre- 
senting the distribution diverge significantly, causing 
considerable error in the evaluation of the mean 
squared voltage. 

Under the high dynamic range conditions it seems 
appropriate to make the same transformation on the 
second term as on the third. Modification of the 
coefficient by some factor to compensate for the 
partial overlapping might be considered, but also 
considering the nature of the experimental data, 
that is, short samples of a nonstationary process, 
and also the relative unimportance of the second 
term, any refinement in the transformation seems 
unjustified. 


6. Discussion 


These somewhat arbitrary operations on the first 
and second terms, together with the accurate trans- 
formation for the high voltage end of the distribution 
where the third term is dominant, provide a method 
of describing what changes would be observed in 
the polynomial representation of the distribution of 
atmospheric radio noise amplitude when observed 
through different bandwidths. The results of carry- 
ing out the complete transformation on representa- 
tive samples of high and low dynamic range noise 
are shown in figure 6 and figure 7. For the low 
dynamic range noise in figure 6, the a, coefficient 
was determined from the rms voltage relationship; 
for the high dynamic range noise in figure 7, it was 
determined by the transformation procedure for 
nonoverlapping pulses. It is apparent that these 
procedures provide good estimates of the distri- 
butions, and are useful and practical tools for use 
when more rigorous methods are not justified. 
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Figure 6. Low dynamic range atmospheric noise distribution 
measured at a wide bandwidth, and estimated distribution 
from simultaneous narrow bandwidth data. 


5§0999—61——_-8 











39 7 7 T T 
| | 
0.98 = } } robability of Exceeding v 
a | | nts Meosured at !I70 Cycies 
| — Estimated from Data at 160 Cycies 
ota token Oct 18,1956,2240 to 2310 
0.95 - ~ - t > soos 

















v,4v/m 


Ficure 7. High dynamic range atmospheric noise distribution 
measured at a wide bandwidth, and estimated distribution 
from simultaneous narrow bandwidth data. 
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7. Appendix: Variation of High Voltage Co- 
efficient and Exponent 


Consider a distribution having a form expressed by 
%1=ay 


P, =¢e- v (6) 
when it is observed through a particular bandwidth. 
If this distribution is composed of nonoverlapping 
impulses, then when observed through a bandwidth 
B times larger, the voltages and probabilities will be 
related by the ratio of bandwidths according to 


P, 

9 ees 

PL=3 
= Br, (7) 


where P; and 7, are obtained at the first bandwidth 
and P, and 7 at the second. From (6) and (7) we 


can write 
v1 
log ( =) 
=\a. 


2 / D> \ ; 
log (—LnBP,) =F log (73) vi 


log (—LInP;)= 


a! bo 


We wish to develop an approximation to the rela- 
tionship between v, and P; in the region of high volt- 
ages and low probabilities which will have the form 


v=aBCyy” 


P,=e"" (9) 
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(10) 


‘ o(_” “ 
log (—Lnl = 57; log ( BC, ) 


where C, and C, are constants to be determined. 
They may be evaluated by selecting two points on 
the P, versus v relationship where the approxima- 
tion is to be exact. Choosing the points Pq, v2.2, and 
Py, 028, we obtain from (8) 

of 

2 


Voa=.B(— In BP 0) 


2 
2 


0.3=aB(—LnBP.,)?. (11) 


Now in the approximation, 
. 2 Voa 
log (—LInl 2=) = 50, log (35; ) 


9 
log (— InP) = 57, log (7 BO 7) (12) 


C, may be eliminated from equations (12), obtaining 


> (LMP x ‘at Vra ? 
log (F258 )=c. Cp og (* *) (13) 


and from (11), 


Ya _(LnBP2.\? 
tng Fe (14) 


and from (13) and (14) we can obtain 


| pat) 
°8 \ Tn BPs 


oa) ‘ 
Ln Ps, 





(15) 


For the evaluation of C,, from (11) and (12) we can 
obtain 


b 
C (—InBP..)? sare i nBP)? 
im cs. oC2 2 


(—LnP2.)? (—LnP3,) ? 





(16) 


Going back to the polynomial term itself, we can 


change the third term to read a3,y’", where 


b= b,,C> 
and 


A3y=AgnBC, (17) 


a;, and b, being the values observed at the narrow 
bandwidth. 
A limiting value of CQ, is obtained as P,, approaches 
4g and this represents the approximation where 
the straight line is tangent to the curve. 
limiting 


This 


value may be found by application of 
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l’Hospital’s rule, and is 


Bane. (18) 
"Ln BP ii 
If we use the tangent approximation and evaluate for 


Poa=P3= e~ *=0.000335 
we obtain 


8 
b,=C.b,= = aB Pe (19) 
and 
8—LiB stat 
3 y= A3y, B ( —— ‘g (20 ) 


Figure 5 shows this tangent line for an assumed 5b, of 
8 and B of 10, and also shows some points from the 
curve in figure 4. The difference between the 
straight line and the points from the curve is zero, 
of course, at P=0.000335, and increases slowly to 
about 1 db at P=0.005 and P=10~°; it is evident 
that there is no significant difference between using 
the straight line and the curve, and the procedure 
using the tangent line provides a simple and accurate 
estimate of the changes in the high voltage end of 
the distribution curve. 
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Excitation of VLF and ELF Radio Waves 
by a Horizontal Magnetic Dipole ** 


Janis Galejs 


(August 29, 1960; revised November 10, 1960) 


The VLF and ELF modes excited by a horizontal magnetic dipole (vertical loop) in the 
spherical shell between a finitely conducting earth and an isotropic sharply bounded iono- 


sphere are shown to have a nearly transverse magnetic character. 
With the exception of the zero order mode, the propa- 


to those of a vertical electric dipole. 


gating modes excited by the magnetic dipole are of slightly higher 


The modes are similar 


amplitudes, provided 


that the far fields of the horizontal magnetic and vertical electric dipoles are equal over flat 


earth in the absence of ionosphere. 


The transient fields generated by a current step in the magnetic dipole are in the first 
approximation similar to the fields generated by a current impulse in a vertical electric 


dipole. 


1. Introduction 


The mode theory of VLF transmissions has been 
developed by Watson [1919], Budden [1953], Schu- 
mann [1954], and Wait [1957, 1960a, 1960b] for 
excitation by vertical and horizontal electric dipoles 
and also by vertical magnetic dipoles (horizontal 
loops). Additional references on this subject have 
been listed by Wait [1960a, 1960b]. 

The reciprocity theorem has been found to be use- 
ful for relating the fields of a horizontal electric dipole 
to fields of vertical electric and magnetic dipoles 
[Wait, 1960a]. The reciprocity relations may be 
also applied to relate the fields due to horizontal 
magnetic dipole excitation to known fields due to 
vertical electric and magnetic dipoles. 


2. Reciprocity Relations 


The reciprocity theorem states that the voltage V, 
induced in antenna 2 by current J, of antenna 1 is 
the same as the voltage V; induced in antenna 1 by 
an identical current J, flowing in antenna 2. 

A vertical electric (VE) dipole of length ds at 2,= 2, 
and a horizontal magnetic (HM) dipole parallel to 
the z-axis (or a vertical loop of area da in the y, 2 
plane) at 2=2, are shown in figure 1, where sub- 
scripts r and ¢ refer to receiver and transmitter 
coordinates respectively. The only nonzero H 
component generated by the VE dipole is Hy. The 
magnitude of the voltage induced in the HM dipole 
Ver is maximum if Hy is parallel to the z-axis, 
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Response of the zero order mode of the m 





magnetic dipole has been calculated. 


which occurs for o=~ and 5 For e'* variation of 


2 
the fields. 
VO" = — tool 5 (29= 21, Zo= 2n) Sin > da, (1) 


where the superscripts hm and ve refer to horizontal 
magnetic and vertical electric dipoles respectively. 
With the same current applied to the HM dipole 
only the vertical electric field #, will contribute to 
Vv” and 





V%= E2"(25= 2,,29= 2) 08. (2) 
z 
4 
HM 
VE ( VM) 
i 
rd 
x 
Zq* 2, 
ee 








x’ 


Fiaure 1. Coordinate systems for defining reciprocity relations. 
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Equating (1) and (2) results in 






eg ( Zo 24,29 Zr) 
—Wwyy sin GH (25= 2,,20= 2,)da/ds. (3) 


A vertical magnetic (VM) dipole (or a horizontal 
loop of area da’ parallel to the 2’, y’ plane) at 
2p=2, and a horizontal magnetic (7M) dipole par- 
allel to the z-axis (or a vertical loop of area da’” in 
the y, 2 plane) at 2=z, may be considered in the 
geometry shown in figure 1. The only nonzero H 
component generated by the V/ dipole is H,,.. The 
magnitude of the voltage induced in the HM dipole 
Vv" is maximum, if H, is parallel to the z-axis, 
which occurs for ¢=0 or x. Hence 


i] Thm __ 


= top? (25= 21, 20= 2) cos ¢da™ (4) 


where the superscript vm refers to the vertical mag- 
netic dipole. With the same current applied to the 
HM dipole only the vertical magnetic field H, will 
contribute to V™ and 

ee". (5) 


Vi = — toll?” (Zp= 21, 20> 


Equating (4) and (5) gives for da” =da’™ 
3*(2,.— 2.,%5= 2,)=—cos oH? (z= Z1y3<9— 2) . (6) 


For spherical coordinates the z component is replaced 
by r component and the p component by the 6 com- 
ponent of the fields. 

The reciprocity relations apply to dipoles in the 
presence of any linear media. The fields of dipoles 
above plane earth derived by solving the wave equa- 
tion should also satisfy (3) and (6). Thus (27) and 
(116) of Norton [1937], which have been derived for 
unit dipoles, check the validity of (3). Some algebra 
is involved in verifying (6) from (35.1) of Sommerfeld 
[1949] and from (2) to (6), (17) to (19), and (30) of 
Norton [1937]. 


3. Fields of the Horizontal Magnetic Dipole 


The dipole fields will be examined with the aid of 
reciprocity relations in an idealized spherical shell 
that is bounded by a homogeneous earth of radius a, 
conductivity o, and dielectric constant «, and by a 
homogeneous ionosphere of radius (a+A), conductiv- 
ity o; and dielectric constant ¢;. The above model 
neglects the effects of the earth’s magnetic field. 

The vertical electric field HE?" is computed by (3) 
applying (5.1), (6.16) of Wait [1960a] and using +e 
relation 


4 pi 


6 cos 6) 


=v{cos 6P,(—cos 6)+P,_1(—cos 6)]/sin 6. (7) 
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This results in 





I (da) ty sin os £edIAe) - v 
; a. if —t})Jn\*r 


2hr, n=0 sin yr sin 6 


[cos 6P,(—cos 6)+P,_1(—cos 6)]. (8) 


E**= 


The vertical magnetic field H?” is computed by 
(6) applying (9.1), (9.2) of Wait [1960a] and using 
the relation (7) as 


8m Sm (2, r) BR 


n=l sin ur sin A 


in, __[(da)wpy cos 


2khrir,n 


d 
[Yale +r ae fale | 
dz, 
[cos 6P,.(—cos 6) +P,_1(—cos 6)]. (9) 
Using the second order representation of the radial 


functions, applying (6.9), 
Wait [1960a] yields 





2f (2d=e™e[REUC,)] +e [REC))'? ~— (10) 
5 ai ee) (11) 
/q 


where q is an integer equal’to either n or m and where 





“= NyC,—C, ; 
R 2(C,) = n,C,+C, (12) 
m— MC, 
R2(Cn)=6, Cin, C, (13) 
C,=v1— 8}, (14) 
C)=v1— (S;/ny) (15) 

+0. : 
$=". (16) 


The plus sign of (11) should be used with g=n, 
while the minus sign is appropriate for g=m. In 
(16) p is a complex number, which is equal to » of 
(8) or equal to » of (9), k=w/c, c=velocity of light. 


The refractive index n of the boundary medium 
characterized by subscript p is 
n? = (o,+Wwe,)/€o (17) 


where co, and e, are the conductivity and permit- 
tivity of the boundary medium and where €, is the 
permittivity of free space (e¢=(367)~! 10~° farad/ 
meter). The boundary medium can be either ground 
(p=g) or ionosphere (p=7). The refractive index 
of ground is 


(18) 


n2, ~o,/(iweo) 


(6.17), (9.3), and (9.4) of 
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where og, is the ground conductivity. The refractive 


index of the ionosphere is 
ni=1-+-4;/(iweo). (19) 


The ionospheric conductivity o; is defined by 


= “te (20) 


where y is the electron collision frequency in the 
ionosphere and w, is the plasma frequency (w? =3180 
N, N=number of electrons per meter’). (C, and 
C’, are roots of the modal equations 


RY(C,)RUC,)=exp (+2ikhC,). (21) 


S, that is related to C, by (14) has a magnitude of 
approximately unity and a small imaginary part 
for propagating modes of low attenuation. With 
S, and S,, determined from (14), v of (8) and u of 
(9) follow from (16). This completes the formal 
specification of the fields EH’ and H™. 

The expressions for and H'™ may be simplified 
by introducing further ‘approximations. The con- 
sideration will be restricted to cases where 


ry =", =a, (22) 
kC.2r<1, (23) 
kQyz.<l, (24) 

Ng>|. (25) 


This gives 


R(C,) ~1 ™ ~~, 26) 

RV(Cn) a4 Mee —], (27) 
1 

ful) =lt+enca~), (28) 

fn(2) On| (mge,) = Cn (iM) , (29) 


Af m(2)/dz=tkngCfin(2) = tkngfin (2). (30) 


Approximating the Legendre functions by [Watson, 
1919; Bremmer, 1949] 


P,(—cos 6) = (2mp sin @)~°> exp | i(o+5)r—0 _ in/t | 
(31) 


results in 


pam —. hm dja yd/r ira; 
E?h"=—sin ok} 5 (dla) BX 


2) bnSt ., @@ 


n=0 








Qrd Tr 


+h ro o 2 . 
Eo” | dja yd)y ‘it 


Ah" —— 
acai Ng V sir sin (dja) h/d 


. 24d 
—i—S§S 


Save ° oe 


cs) 
‘in m| 


where the distance along the curved earth is 


d=aé, (34) 
and where 

. 27a 
Em 2rl | da)n cit 
0 rn% 
is the vertical electric field of the source at a distance 
d in the direction of maximum intensity on a per- 

rectly conducting plane earth. 


4. Comparison of Harmonically Excited 
Dipole Fields 


The fields of the horizontal magnetic dipole £?” 
and H?” will be compared first with the correspond- 
ing field components of the horizontal electric dipole. 
It follows from (9.37) and (9.42) of Wait [1960a] 
that for the n™ and the m™ mode 


Eh | qrhe 

=| =—n, tan d= | (36) 

ES” in E* vg in ; 

H™| H™| a 

Fam =n, cot d=, (37) 
“0 jm Bri. 


where the superscripts he and ve designate field com- 
ponents of horizontal electric and of vertical electric 
dipoles respectively. /}” is defined in (35), E}° is the 
corresponding expression of the vertical electric dipole 
defined by (6.25) or (9.38) of Wait [1960a]. The 
approximation C,~1 was used in the expressions 
for #”* and H” in order to make them consistent with 
the derivation leading to (32) and (33). C,=1 con- 
stitutes a better approximation than f,(z)+1 and 
the assumption f,(z)~1 [Wait, 1960a] should be 
followed by C,~1. 

Comparison of (32) and (33) shows that 7H} 
several orders of magnitude smaller than H?” bacedes 


hm j 


ng'<1 (38) 
and because 
1—S2,.<1 (39) 


for propagating modes of low attenuation. The 
Hh" fields of a given mode are proportional to 
ng'~ywlo,. The fields are decreased with decreas- 
ing frequency w or increasing ground conductivity 
g,. Similar comments apply to the transverse electric 
(TE) field components that may be derived from 
H,. The TE fields will not be considered in more 
detail. 

The component £2" that is associated with trans- 
verse magnetic (7M) fields may be compared with 
the corresponding component /%° of a vertical elec- 
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tric dipole using an expression of the E%* fields de- 


rived under similar assumptions [Wait, 1957]. Thus 

faut sin o£ 5” BE ikI’da in oE" (40) 
oe") = — ——- Et") =——— s > 
: n S,, oO é n S, i ds : n 


Provided that the fields of the two antennas are the 
same in the given direction over a perfectly conduct- 
ing plane earth (sin ¢£}"= E%*), (40) may be simpli- 
fied to 

E** | 


—S,,. (41) 


n 





| 





Substituting the reciprocity relation (3) in (40) fol- | 


lows that (/%*)/(nH3°) also satisfy (41). The latter 
ratio, obtained from (13), (15), and (16) of Wait 
[1960b], when used in conjunction with (3) provides 
an alternate way of obtaining (40) and subsequently 
(41). 

A more exact treatment of the horizontal magnetic 
dipole fields cannot be expected to lessen the im- 
portance of 7M relative to TE field components. 
A better approximation to the 7M fields of the hori- 
zontal magnetic dipole may be worked out from the 
exact ratio E’ ‘Ee derived in the appendix. As long 
as the source and observation point are at equal 
radii 7 the ratio involves only the azimuthal 6 func- 
tions explicitly. However, 
S, and of £%* depends on the approximations to the 
azimuthal and radial functions. 

Numerical values of S,, are available from mode 
calculations of the vertical electric dipole [Howe 








and Wait, 1957]. For o,=o the values are as 
follows: 
h/» 3.5 3.5 7 
av¢|(w>h) 0.1 0.01 0.1 
1 0. 9975-i 310-4 0. 997-i 10-8 0. 995-i 610-5 
2 . 98-i 3X10-3 .97-i 8X10°3 . 992-1 4X10°4 
3 | .93-i0.01 . 91-1 9X10-% . 98-i 10-8 
4 . 87-1 0. 02 . 85-1 0. 01 . 985-i 2X10-8 





The n=0 mode is severely attenuated for fre- 
quencies above a few kilocycles [Wait, 1957]. In 
the lower frequency range the first order perturba- 
tion solution of the modal equation may be used 
(Schumann, 1954], which gives 


’ 3-4 
he +i. A | 2 eee ne 24h Se) (42) 


The cutoff frequency of the nth mode is 


NTC 
nes 





(43) 


The perturbation method applies for n #0 if C,,~a,/w. 
Substituting this for Cz 


‘ in the right-hand side of 


the accuracy of the roots | 














(42) and expanding the square root expressions 


Unt = re. : 
« 24 (7 mtn 3 )] 
“a ies 


The latter expression may be also obtained by solving 
(42) as a quadratic and by subsequently ignoring the 


(53 





CO, = (44) 


second order perturbation. Applying (14), (42) 
simplifies for n=0 and |n,|, |n;|>>1 to 

c ue ; 

So= 4/1 +7 —-+— | 45) 

hiw Lats | _ 


For no>n,; and o,=1.2X10~*(w?/y=1.35X10°) Sp 


has the following values 


10? 
1.46 e-i6e 


10° 
1.14 e-i7° 


w| 10 
S)/1 .04 e—i2° 





The above tabulations show that |S,! or the ratio 
|Eee/E?"| are decreased with increasing mode index 
n and, with the exception of So, also with decreasing 
frequency. The phase angle between the propa- 
gating modes of the two dipoles is negligible. The 
largest difference between the fields of the two dipoles 
may be expected for the lowest propagating fre- 
quencies or near the mode cutoff of n#0 modes. 
The values of S, near the mode cutoff may be deter- 
mined from an approximate solution of the modal 
equation. The first order perturbation solution is 
applicable only to the lower order modes. A differ- 
ent approximate solution is obtained by observing 
that S, is small near the cutoff (S,=1/(2 ka) at the 
mode cutoff for perfectly conducting ground and 
ionosphere according to Schumann [1954]). Hence 
one may look for solutions of the modal equation of 
the form 
1—S?= 


Jail (46) 


(12), letting 


C,=v1—S; 


where |A,|<1. Substituting (46) 
p=7 and defining 


Ego OV 5 
Lan—-weSp (47) 
i p 


results after neglecting terms with A? and _ higher 
powers of A, in 


AnvL 


aa 
VL: 


substituting (46) in 


Rt=i(2L—i—2yLVL— ae L,+A,LZ2. 


(48) 


—} 


Considering 2>1, (12) and 


letting p=g gives 
?2=~1—2n; (49) 


1 git/4(1 +A.) =G,+A,G>. 
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The solution of the modal equation (21) becomes 


21(kh—mn) —InlGi 


Pca 
aikh-+ 424 Oe 


(50) 


It is simplified at the cutoff of the respective mode 
by applying (43) to 
4=— i InLiG; n 


dine 124 





(51) 


The calculated curves of |S,,/ and of its argument are 
depicted in figures 2 and 3 for infinite ground con- 
ductivity. |S,| is smallest for small values of 
I=rcv/ (hw) (or for large ionospheric conductivities) 
and for the higher modes (higher cutoff frequencies 
w,—nme/h). The argument of S, exhibits a rela- 
tively small variation with n. The increase of |S, 
and of arg S, with increasing J signifies the increase 
of mode attenuation with * decreasing ionospheric 
conductivity. The increase of S, from the ¢,= © 
to the o,=0. 07 /him curve (¢,=10~* mho/m for h=70 
km) is less than 20 percent for the n=1 mode in 
figure 4. This relative increase is even smaller for 
the higher n values. 

The: amplitude ratio between electric and magnetic 
dipole fields in (41) will be larger for propagating 
modes than near the mode cutoff. The cutoff ratio 
is larger than 0.25 for the n=1 mode and larger than 








0 n j 
0.0 0.1 10 
We Vv 
ae oP 
FicurE 2. Magnitude of the constant 8, near mode cutoff for 


perfectly conducting te os 











0.15 for the first 8 modes with 750.01 (or for 
¢,>1.2X10-* with h=70 km). It may be noted 
that this comparison is made for dipoles that exhibit 
equal far fields above a perfectly conducting plane 
ground in the absence of ionosphere. The knowledge 
of the above amplitude ratios would enter in a com- 
parison of antenna radiation efficiencies. However, 
no such calculations have been attempted in this 
paper. 











Figure 3. Cutoff arguments of the constant S, for perfectly 


conducting ground. 
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5. Transient Fields 


The transient fields of the horizontal magnetic 


dipole may be computed as the Fourier integral of 


(32) for a specified spectrum of the antenna current. 
For a step of current the far field approximation (32) 


gives 





__ @ ,od 
Epm(e) =Kyio D3 3,Sn5 ee (52) 
n=0 
with 
- sin a r 
— 7 $1 di - (53) 





2hey 2rac sin (d/a) 
and with 6:=0.5 and 6,..=1. 
The response of the n=0 mode may be reduced 


to tabulated Fourier integrals after approximating 
S, of (52) by a power series. With n,>n;>1 (45) 


reduces to 
; 1 
So ~~ y! - —_—————* 
hy WUC 7 


Expanding ¥S,) and S) in powers of (1/yiw) and 
ending the expansions with the quadratic terms 
results in 


(54) 


; > senia = 1 2 _iwd _ vied 
E(w) = Ke / (2cB?) \ i ———— — é c cB R 
28 8B°V Ww 
(55) 
where a 
B= 2hwy uo}. (56) 


Applying the transform pairs (806), (801), and 
(807) of Campbell and Foster [1948] gives 


Kel/u-0-25/2 1 ee r 
E*™(r) | tba EAC ——)— | (57) 
} 27° u ue 





2a! ire 
where 
2=T/a, (58) 
7=(t—d/c)>0, (59) 
ded? vod” 
a(5 ~ Aah? ath? a 
9 2 
yan 2 _ 80h [wo (61) 


d ca V 


The power series expansion of (54) is applicable to 
frequencies where w>wo=(h7uo;)~'. Hence (57) will 
be most accurate for r< <22/w, or for 


2h2o; [6 


LO <H=2r EE 


Considering h=90 km, d=10* km, o;=10~* mho;m 
as an example, 2)»=2.3 with w=2.4. The representa- 
tion of the transient peak in figure 5 will be accurate 
even for this pessimistically low o,; value at d=10* 


2 
=ru?/8. (62) 
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Figure 5. Response to a current step of the n=O mode of a 
horizontal magnetic dipole. 


km. Higher o; and smaller d result in larger a 
(and ~), which improves the accuracy of the transient 
tail representations. 

It follows from reciprocity relation (3) that for 
equal current waveforms /’”"(t) and J°¢(#) in the 
HIM and VE dipoles 

Er" (t) =—pwoH (t) sin ¢da/ds. (63) 
A current in the VE dipole that is proportional to 
the time derivative of the current in the HM dipole 


T(t) = TI" (t), (64) 
results in 
E(t) =—p sin ¢daHy (t)/(Tds). (65) 


The response /?”"(t) to a current step J,u(t) should 
be proportional to the response H%(t) to a current 
impulse J,76(t). Assuming that (64) applies and 
assuming further that 
So~1 (66) 
(this is obviously incorrect at the lower frequencies), 
it follows from (40) that 
E*(t) ~—sin oE}* (t)da/(eTds). (67) 
The impulse response of a vertical electric dipole 
has been calculated considering factors proportional 
to finite powers of S) equal to unity [Schumann, 
1952; Wait, 1960bj. The calculated impulse response 
of the VE dipole [Schumann, 1952; Wait, 1960b] 
will be compared with the step response of the HM 
dipole (57) by means of (65) and also by means of 
the approximate relation (67). The leading term 
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of the impulse response (42) of Wait [1960b] (eqs 
(38), (39), (40), and (42) of Wait [1960b] should be 
multiplied by z/2), when substituted in either (65) or 
(67) results in £?"(t) that is equal to (57) with the 
two terms proportional to u~! and u~? of the square 
brackets set equal to zero. There is no agreement 
between the higher response terms because of the 
S) approximation as discussed earlier. The impulse 
response of Schumann [1952], when substituted in 
(67) is the same as (57) with u= @, 

The first order perturbation solution of the modal 
equation becomes inaccurate for the higher fre- 
quencies involved in calculating the transient re- 
sponse of the n#0 modes. Even the simplest (and 
inaccurate) expressions of S, that may be obtained 
from (44) and (14) for |n,|>>|n;|>>1 involve 
integrals similar to those encountered in the transient 
analysis of lossy rectangular waveguides [Cerillo, 
1948]. The transient response of the n#0 modes 
has not been calculated in this paper. 


6. Appendix. Radial Electric Fields of 
Vertical Electric and Horizontal Magnetic 
Dipoles 


The fields of a vertical dipole may be derived from 
a single scalar function. Applying (29) of Schumann 
[1954] to (3) 


Os. Ou, 
Ee ds, y, (sin 4 >a) am 
io ; a , Ow x 
‘ fe “OS 1 98 6 om 
rk? cos ¢ daz sin 9 








where for uw. the coordinates of the source and of the 
observation points are interchanged relative to 4%. 
As long as the source and the observation point are 
approximately at the same radius (7,~r,), both u 
functions are the same. The ratio of the vertical 
electric fields of horizontal magnetic and vertical 
electric dipoles for equal currents on a perfectly con- 
ducting plane ground is 


tg = ids; 


Ek" kda; ta 





With 

U; ~U2~P, (—cos @) (70) 
(68) becomes 
Ee cc ” (v+]) sin 6 P, (—cos 8) 


Et Eb cos ork cos 6 P, (—cos 0)+P,_1 (—cos 8) 


(71) 


The earlier approximations will reduce (71) to (41). 
However, more accurate v values [Wait, 1960a, ch. 
12] and a better approximation to P, (—cos @) may 
be used in (71) to obtain a more exact E/E?” ratio. 
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The relation of h,,,,.F2 to M(3000) F2 and h,F2, J. W. Wright 
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409 (July 1960). 

Several easy methods of inferring directly from ionograms 
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M(3000)/2 or h,F2, are compared with values of hmazl’2 
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profile with the average of the profiles for the individual 
days. The differences are small compared with the disper- 
sion of the individual profiles about their mean. It is con- 
cluded that the method is valuable for obtaining a world- 
wide N(h) morphology with a minimum expenditure of effort. 


Spiral patterns in geophysics, V. Agy, J. Almospheric and 
Terrest. Phys. 19, No. 2, 136-140 (Oct. 1960). 

Recent analyses of magnetic and ionospheric data, notably 
by A. P. Nikolski, have resulted in spiral ‘precipitation’ 
patterns which lead the authors to claim support for St6rmer’s 
theory of the aurora. Although it may be possible to argue 
against these claims by attacking the methods of analysis 
and/or the Stérmer theory itself, an entirely different approach 
is used in the paper presented here: an examination of per- 
tinent points of St6rmer’s theory shows that the analytical 
spirals mentioned cannot be Stérmer spirals and this con- 
clusion holds regardless of the soundness of the analyses and, 
indeed, of the validity of Stérmer’s theory. 


Supplementary world maps of F2 critical frequencies and 
maximum usable frequency factors, Donald H. Zacharisen, 
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This report supplements NBS Tech. Note 2 (April 1959), 
and completes the basic data required for F2-layer maximum 
usable frequency predictions. Prediction charts are given 
for the months of February, April, May, August, October, 
and November. Auxiliary charts are included to aid in 
predicting F2-layer MUFs. 


The four parameters used for predicting MUFs are foF2 and 
the 4000 km MUF factor for a twelve-month running average 
Zurich sunspot number of 50, and the rates of change of foF2 
and 4000 km MUF factor with sunspot number. The first 
three parameters are presented in map form for each even 
hour of Greenwich Mean Time. The fourth parameter is 
presented on a chart of geomagnetic latitude and local time. 


The height of maximum luminosity in an auroral are, F. E. 
Roach, J. G. Moore, E. C. Bruner, Jr., H. Cronin, and 8. M. 
Silverman, J. Geophys. Research 65, No. 11, 3575-3580 
(Nov. 1960). 

The height of maximum luminosity of an auroral are is 
estimated from simultaneous observations at three stations 
in western United States during a night of general auroral 
activity (November 27-28, 1959). Photometrically this are 
is characterized by a selective enhancement of the [OI] 6300 A 
line. From twenty-four individual measurements the height 
is found to be 412 km with a standard deviation of +23 km 
for one observation and +5 km for the mean. The geo- 
graphical position of the are, its orientation, and its movement 
during the night are discussed. 


FM and SSB radiotelephone tests on a VHF ionospheric 
scatter link during multipath conditions, J. W. Koch, W. B. 
Harding, and R. J. Jansen, TRE Trans. Commun. Systems 
CS-8, No. 3, 183-186 (Sept. 1960). 

Experiments have been carried out on an ionospheric-scatter 
link to observe the effects of long-delayed multipath signals, 
caused by F; propagated back scatter, on the intelligibility 
of voice communication. Frequency modulation and single- 
sideband modulation equipments were used for the tests. 
During periods when the back-scatter signal levels approached 
the level of the normal ionospheric-scatter signals, the 
frequency-modulation voice transmissions were unintelligible ; 
however, under the same conditions, single-sideband voice 
communication intelligibility remained at almost 100 per 
cent although there was some loss in quality. 


Radio refractometry, Jack W. Herbstreit, NBS Tech. Note 66 
(PB 161567) (July 1960) 50 cents. 

The optical refractive index is known to be determined 
principally by the temperature and pressure of the atmos- 
phere, whereas the radio refractive index is, in addition, 
affected by the water content of the atmosphere, the relation- 
ship between these quantities being expressed in the following 
way: 





Mate) 10=(S) (P+[4810e/Tl) (1) 


where total air pressure P, and water vapor e, are in millibars, 
and the temperature 7’ is in degrees Kelvin [Smith, 1953]. 
The quantities P, e, and 7 have long been routinely measured 
at the surface of the earth by standard weather bureau sta- 
tions. For some time they have been measured from the 
surface up to great heights using balloon-borne radiosonde 
equipment at a large number of places over the earth’s 
surface. More recently, equipment has been developed to 
measure rapidly and directly the radio refractive index of the 
atmosphere using radio techniques. The measurement of the 
radio refractive index properties of the atmosphere and the 
application to radio propagation problems is the subject of this 
paper. 
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